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The infra-red spectrum of methylphenylacetylene between 500 and 2300 cm has been 
measured with the aid of a prism spectrometer. The results are compared with the Raman 
spectrum of this compound, previously measured by the authors. The results are inconclusive 
as regards the assignment of the fundamental whose overtone, according to Badger’s hypothe- 
sis, is responsible for the resonance splitting of the frequency near 2230 cm~, 


INTRODUCTION 


ADGER' has suggested that the appearance 
of two strong frequencies in the 2200-cm- 
region in the Raman spectra of nearly all of the 
disubstituted acetylenes is attributable to a 
Fermi resonance interaction between the funda- 
mental frequency expected in this region and an 
overtone of a lower frequency near 750 or 1100 
cm~, Since the latter frequency is inactive in the 
Raman spectrum, he suggested that an infra-red 
study would best decide which of the two fre- 
quencies more nearly satisfied the condition for 
resonance. 

At the time the present investigation was 
begun, the authors had just completed a study of 
the Raman spectra of a series of derivatives of 
phenylacetylene and had found that the second 
overtone of a frequency near 750 cm™ was 
considerably greater than the mean of the 
doublet near 2230 cm. It was therefore of 
interest to measure the infra-red spectrum of one 


of the compounds in order to find if the expected. 


frequency in the infra-red would better satisfy 
1R. M. Badger, J. Chem. Phys. 5, 178 (1937). 


the frequency condition for resonance. Methyl- 
phenylacetylene was chosen for the study because 
in this compound the two lines in the 2200-cm-! 
region were well separated and approximately 
equal in intensity. 


EXPERIMENTAL PROCEDURE 


The methylphenylacetylene was synthesized 
by the action of dimethylsulfate on sodium 
phenylacetylide in ether solution and was purified 
by distillation in an efficient heated column 
using a high reflux ratio. The boiling point of 
the liquid was 72°-73°C at a pressure of 15 mm of 
mercury. The sample was the same as that used 
in obtaining the Raman spectrum. 

The infra-red spectrum was obtained with the 
recording prism spectrometer described by 
Randall and Strong? with the addition of the 
periodic radiometer designed by Firestone* to 
eliminate extraneous effects. The slit width 
varied from 0.5 mm at 500 cm to 0.10 mm at 
2000 cm-!; near 1100 cm it was 0.2 mm. 

2H. M. Randall and John Strong, Rev. Sci. Inst. 2, 585 


(1931). 
3F, A. Firestone, Rev. Sci. Inst. 3, 163 (1932). 
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134 mM. 


Data were first secured with the empty cell 
alone in order to make possible corrections for the 
absorption of carbon dioxide and water vapor in 
the atmosphere. The thickness of the liquid in 
the sodium chloride absorption cell was 0.1 mm. 


TABLE I. Infra-red and Raman spectra of CGHs;C =CCHs3. 


INFRA-RED RAMAN 
APPROXI- DISPLACE- DEPOLARIZA- 
FREQUENCY MATE MENT ESTIMATED TION 
(cm™) INTENSITY (cm~) INTENSITY FACTOR 
149 
298 3 0.8 
388 4 0.7 
401 2 — 
510 1 _ 
526 8 528 3 0.9 
546 4 
594 6 
611 2 
624 4 626 3 —_— 
631 + 
649 6 
702 10 701 3 —_ 
754 10 758 2 = 
825 6 
839 4 
909 10 
963 6 953 1 —_ 
981 2 971 6 0.3 
999 4 1001 7 0.3 
1024 8 1030 2 — 
1064 8 1053 3 —- 
1 1 + 
1196 6 1181 i} <0.5 
1256 4 1245 1 — 
1262 8 0.3 
1300 Zz 1286 2 — 
1326 2 1328 1 — 
1359 6 1382 4b 0.9 
1418 8 1436 3 0.7 
1468 8 1450 3 ? 
1489 4 0.7 
1500 4 
1572 6 1601 9 0.7 
1667 6 
1821 4 
1890 4 
1957 4 
2137 6 
2214 9 0.6 
2254* 10 0.6 


* A frequency 2231(4) was observed on three different spectrograms 
and 2254 appeared resolved into the doublet 2246(10), 2262(6) when the 
camera lens was diaphragmed to 1.0 cm. 


MURRAY AND F. F. CLEVELAND 


The source of radiation was a Nernst glower 
heated by current from a 220-v, a.c. power line. 


RESULTS 


Comparison of the results with the Raman 
data obtained by the authors is made in Table I. 


DISCUSSION OF RESULTS 


The frequency in the 1100 cm region which 
most nearly satisfies the frequency condition for 
resonance is the weak one at 1095 cm~". The first 
overtone of this frequency, 2190 cm, is seen to 
have a value even less than that of the smaller 
frequency of the doublet, 2214, 2254. The second 
overtone of the 756 cm frequency, 2268 cm-, 
on the other hand, is greater than the larger 
frequency of the doublet. The results are thus 
inconclusive. 

Soon after these measurements were com- 
pleted, a paper by Crawford‘ appeared which 
reported the results of an infra-red and Raman 
study of dimethylacetylene. A weak frequency 
was observed at 1126 cm-. The first overtone of 
this frequency, 2252 cm, falls between the 
doublet frequencies, 2233, 2310, but is somewhat 
less than the mean, 2272 cm-!. The mean of the 
doublet, 693, 774 cm, is near the value 725 
cm which Crawford assigns as one of the 
fundamental frequencies of dimethylacetylene. 
The second overtone of this frequency, 2175 
cm-', is much lower than the mean, 2272 cm. 
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The Rotational Structure of the Hydrazoic Acid Bands in the Photographic Infra-Red 
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The third and fourth harmonics of the N—H stretching frequency in hydrazoic acid, at 
9547.31 and 12,412.19 cm™, respectively, have been studied photographically under high dis- 
persion. The rotational structure has been well resolved and completely analyzed, to yield the 
rotational constants in the ground and excited states. The effective moments of inertia in the 
ground state have been combined with other information to determine the structure of the 
molecule; the best values of the molecular parameters are: 


Ni-H =1.012A 
Ni—N2=1.241A 


N2—N3=1.128A 
AHN;N2=110° 52’. 


LECTRON and x-ray diffraction experi- 
ments have fully confirmed the predictions 

of the modern structural theory concerning the 
structure of the azide group in ionic crystals and 
in covalent compounds.! The infra-red spectra of 
covalent azides have not, however, been satis- 
factorily correlated with acceptable structures. 
It has been asserted that the rotational structure 
of the hydrazoic acid bands in the photographic 
infra-red demands the existence of two kinds of 
HN; molecules, one of which is structurally 
inacceptable, and that the envelopes of methyl 
azide bands in the same region indicate that this 
molecule has a linear C—N—N-—WN group, con- 
trary to the electron diffraction results.? It has 


consequently been thought useful to reinvestigate 


these spectra in order to determine whether or 
not these reported inconsistencies are real. The 
present paper presents the complete analysis of 
the rotational structure of the two strongest 
bands of hydrazoic acid in the photographic 
region, and demonstrates the fact that the 
resulting rotational constants determine a struc- 
ture which is in complete harmony with theo- 
retical predictions and with previous structure 
determinations in other covalent azides. 


EXPERIMENTAL 


The infra-red spectrum of hydrazoic acid was 
photographed by the author in the spectroscopic 


* National Research Fellow in Chemistry. 

'The considerations leading to a prediction of the 
structures of the azide group in covalent and ionic com- 
pounds are summarized on page 185, L. Pauling, The 
Nature of the Chemical Bond (Cornell University Press, 1939). 

* A discussion of these results, with references to the 
Original papers, is given in an interesting review article by 
Verleger, Physik. Zeits. 38, 83 (1937). 
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laboratory of the Gates and Crellin Laboratories 
of Chemistry of the California Institute of 
Technology. The acid was prepared from recrys- 
tallized sodium azide (Eimer and Amend) by the 
method of Dennis and Isham.* The gas was swept 
from the generator in a current of dry nitrogen, 
through a 50-cm calcium chloride drying tube, 
and into a trap cooled in liquid air. This trap, 
attached to the generator by a Pyrex standard 
taper, was then transferred to a similar fitting on 
the six-meter Pyrex absorption tube, which was 
at once evacuated and sealed off while the 
explosive acid was still at liquid-air temperatures. 
The gas pressure was subsequently regulated by 
adjusting the temperature of the liquid acid in 
this trap. The whole of the absorption tube was 
shielded from stray radiation with black cloth, 
and a filter which cut out all wave-lengths 
shorter than 7000A was interposed between the 
source and the tube to protect its contents from 
possible photochemical decomposition in the 
direct beam. 

With the trap at room temperature, explora- 
tion photographs were taken with a Bausch and 
Lomb glass Littrow spectrograph set to record © 
the region between 7500 and 12,000A on Eastman 
I-Z plates. Four bands were observed, and the 
two stronger ones were then rephotographed in 
the first order of a twenty-one-foot concave 
grating spectrograph. For the band at 10,477A, 
the trap was held at 0°C, and for the weaker 
8058A band it was held at 25°C. Wave-length 
measurements on the prism spectrograms were 


3L. M. Dennis and H. Isham, J. Am. Chem. Soc. 29, 
216 (1907). 
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_ Fic, 1. The hydrazoic acid band at 9547 cm, and its analysis into constituent sub-bands and branches. Calculated 
intensities of the K-structure lines are indicated by their heights, and the intensities of the J-structure lines are given as 


envelopes. 


made with respect to infra-red emission lines of 
copper, sodium, potassium, barium, and calcium, 
and on the grating spectrograms with respect to 
iron arc standards in the second and third 
orders. Eastman 144-R and I-Z plates, hyper- 
sensitized in ammonia, were used, and higher 
order spectra from the grating were removed by 
the filter already mentioned. 


ROTATIONAL STRUCTURE 


Since one may confidently expect the struc- 
tures of hydrazoic acid and its methyl derivative 


' to be analogous, the results of the electron 


diffraction study of methyl azide* should provide 
a basis for the calculation of approximate 
moments of inertia for the acid itself. Such 
approximate moments enable one to estimate 
as probable values for Badger’s rotational 
parameters S and p° the following: 


=—0.999; p=+61.4. 


in a a and L. O. Brockway, J. Am. Chem. Soc. 59, 
’R. M. Badger and L. R. Zumwalt, J. Chem. Phys. 6, 
713 (1938). 


These values constitute a numerical statement of 
the fact that the rotational levels of hydrazoic 
acid should approximate those of a symmetric 
rotator of an extreme spindle shape. Hence it 
should be possible to treat the hydrazoic acid 
bands as symmetric rotator bands, except that 
lines receiving a large portion of their intensity 
from levels of high J and low K may show a 
splitting or fuzziness as the result of asymmetry. 

The bands at 10,477A and 8058A undoubtedly 
correspond to the third and fourth harmonics of 
the N—H stretching frequency, and both the 
third and fourth harmonic states of this vibration 
belong to the irreducible representation A’ of the 
group C, of the planar hydrazoic acid. Molecular 
symmetry then does not require that these bands 
present rotational structures of pure types, and 
if we imagine the direction of electric moment 
change to be the N—H bond direction, parallel- 
perpendicular hybrid bands will in fact be 
anticipated. 


6 This molecule is so nearly symmetric that the type B 
and C bands may well be called perpendicular, and the 
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General features of observed bands 


These general expectations are well fulfilled by 
the observed bands. The third harmonic band is a 
hybrid of predominantly parallel structure, 


whose most striking feature is provided by the > 


fact that the various parallel sub-bands, which 
usually coalesce to form a single parallel band, 
are here clearly separated as a result of the strong 
dependence of the small moment of inertia upon 
vibrational state ; unusually marked convergence 
is correspondingly present in the perpendicular 
component. Fig. 1 shows a representation of the 
complete band, with its analysis into constituent 
branches. The coarse structure lines are sepa- 
rately represented and named, while only the 
envelopes of the closely spaced @Rx and °@Px 
branches are shown; though these latter are not 


specifically named, their naming is uniquely 
implied by the associated °Qx branches. The 
nomenclature of lines will be evident from the 
following: the large letter describes the behavior 
of J in the transition, the left-hand superscript 
describes the behavior of K, the right-hand 
subscript gives the value of K’’, and the value of 
J” follows the large letter in parentheses. 

The fourth harmonic is a pure parallel band, 
and hence much simpler in appearance than the 
band of Fig. 1. As might be anticipated, the °Q 
branch sweeps back even more rapidly than in the 
third harmonic. Since the appearance of this 
band can be readily visualized from a con- 
sideration of the parallel component of the third 
harmonic band, it has not been separately 
illustrated. The frequencies and identifications of 


TABLE I, The third N—H harmonic band of hydrazoic acid at 9547.31 cm. 


J | | | PRA) | PPA) | | K °0K “OK 
6 9539.95 0 — 9564.82 
7 39.09 1 9527.07 | 9544.60 | 597.44 
8 9474.37 38.35 2 485.10 | 537.50 | 622.02 
9 73.69 | 9497.53 9545.30 — 3 438.00 | 522.92 | 642.93 
10 72.89 96.92 46.02 36.71 4 384.25 | 504.74 | 657.30 
i . 71.88 95.96 | 9514.03 46.73 35.92 5 326.40 | 480.76 | 666.32 
12 71.14 95.14 13.24 47.51 | 9527.74 35.02 6 452.62 

13 70.32 94.22 12.50 26.89 34.15 7 420.20 

14 69.40 93.32 | —— 26.08 33.41 8 384.14 

15 68.60 92.46 — 25.23 32.50 

16 67.74 91.60 10.04 24.45 31.68 

17 66.99 90.73 09.20 23.64 | 9558.21 30.84 

18 66.04 89.91 08.27 — 58.96 30.02 

19 65.19 89.02 07.41 22.03 59.69 29.16 

20 64.34 88.11 06.56 21.45 60.40 28.35 

21 63.55 87.36 05.72 20.31 61.10 27.30 

22 62.64 86.39 19.39 61.78 26.46 

23 85.62 03.97 18.47 62.48 25.58 

24 84.76 03.14 17.57 63.17 24.74 

25 83.77 02.29 16.64 63.88 23.86 

26 82.99 01.36. 15.68 64.61 

27 82.11 00.40 14.79 65.25 

28 — 65.97 

29 80.08 66.65 

30 79.33 67.32 

31 78.44 68.06 

32 77.47 68.76 

33 76.59 69.44 

34 70.12 

35 70.91 

36 71.60 

37 72.29 

38 72.95 

39 73.60 

40 74.25 

41 74.92 

42 75.54 

43 76.21 

44 76.80 

type A bands parallel. For a discussion of pure and hybrid bands and their dependence on molecular symmetry see 
E. H. Eyster, J. Chem. Phys. 6, 577 (1938). 
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TABLE II. The fourth N—H harmonic band of hydrazoic acid at 12,412.19 cm=. 


J @Rx(J) RJ) K °0K 
6 12,405.40 1 12,410.60 
7 12,406.03 04.72 2 399.91 
8 06.71 03.99 3 382.21 
9 12,350.64 07.54 03.12 4 357.69 

10 49.69 08.39 02.18 5 326.59 

11 48.92 09.14 12,419.73 01.34 6 289.66 

12 48.02 09.87 | 12,390.03 20.33 00.67 

13 47.10 10.52 89.26 21.05 

14 46.26 11.34 88.45 21.77 

15 45.47 | 12,394.41 12.04 87.57 22.31 

16 44.62 95.14 | 12,369.13 12.80 86.69 23.06 

17 43.73 95.77 68.18 13.40 85.87 23.97 

18 42.86 96.50 67.19 14.19 84.91 24.71 

19 42.06 97.19 66.34 14.90 84.04 25.30 

20 41.15 97.93 65.52 15.54 83.17 25.97 393.58 

21 40.42 98.68 64.70 16.25 — 26.61 92.76 

22 39.38 63.83 16.93 81.45 27.31 91.82 

23 38.47 62.87 17.61 80.52 90.80 

24 37.56 61.94 18.29 79.72 

25 36.66 61.03 19.03 78.87 

26 35.77 60.10 77.96 

27 34.88 59.21 

28 34.01 58.34 

29 33.07 56.64 

30 31.99 55.64 

31 31.05 54.69 

32 30.21 53.83 

33 29.35 52.83 

34 28.37 

35 27.32 


the 250 lines that were measurable in these bands 
are given in Tables I and II. Missing lines and 
branches in these tables are ones which are 
completely obscured by stronger components of 
the bands; many of the reported lines are 
blends, and subject to unavoidable errors of 
measurement for that reason. This is especially 
true of lines in the third harmonic band, where 
there are several regions in which the overlapping 
is so multifold as to produce utter confusion in 
the fine structure. 


Rotational constants 

In the absence of strong perturbations by 
adjacent vibrational states, the rotational energy 
of a symmetric rotator may be conveniently 
written as follows: 


The primed summation indicates the omission of 
the constant term X99’, and v is to be considered 
as the totality of quantum numbers required to 
specify the vibrational state; in our case it will 
simply be 2;, the quantum number associated 
with the N—H stretching frequency. A formula 


of this type should be expected to describe the 
majority of the rotational levels of a molecule so 
nearly symmetric as hydrazoic acid; the leading 
coefficients will then be related to the effective 
moments of inertia as follows: . 


Xo" (2) 

where 

(3) 


It was found best to evaluate Xo,° and X 9° in 
the third harmonic band, with the aid of the 
following relations: 


AF” (K)/(2K+1) =X +2X o2°M+ 
M=(K+})'+3 


(4) 
(K)/4K =X +2X 
M=K?*+1. 
The term differences of (4) are given by: 
AF’ (K) =°Qx— ’Ox+1 =FOx—°Ox+1, 3) 


AoF”(K) ="Qx-1—?Qx+1. 


| 
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Eqs. (4) were tested with the experimental 
combinations, and found to represent the data 
within experimental error. The constants of the 
equations were then determined by the method 
of least squares to be: 


X 19.946  Xo2"= — 0.05616 


The dependence of Xo," and X02” upon v was 
established by the determination of Xou', 
Xo2*, and Xo! with the aid of the relation: 


(°Q:—°Qx)/(K?—1) —Xor' 
(6) 
M=K?+1. 


This relation was first tested with the observed 
combinations in both the third and fourth 
harmonic bands, and found to be well fulfilled; 
deviations from linearity in M were within the 
error of observation. The constants were then 
determined by the method of least squares, and 
it was found possible to represent the dependence 
of the terms Xo,” and Xo2” upon the vibrational 
quantum number », of the N—H stretching 
frequency as: 


X91"? = X 919 — = 0.905 cm—!, (7) 
X 92"! =X 92° — = — 0.000904 cm—!. 


It must be realized, of course, that the “‘lines”’ 
used in the above formulas are actually whole 
branches, and that the appropriate averages 
over their J structures should appear in Eqs. (4), 
(5), and (6) to render them precise. This effect 
will be small, however, since these ‘“‘lines” are 
surprisingly sharp. It must also be admitted that 
the constants Xo.” and Xo” which we have 
determined in this way may contain small contri- 
butions from cross product terms like Xi,” and 
X12", but we shall find in the following section 
that these cross product terms appear to be 
exceedingly small. 

Determination of the constants Xi9° and X 20° 
was made by the usual method, with the aid of 
the expressions : 


AoF x” (J) /2(2T +1) =X 10° + K2X 11° 
(8) 
M = {(J+4)?+ 3}. 


The combinations of (8) are the familiar ones: 
AoFR” (J) (9) 


When these combinations were evaluated from 
the three available values of K in the fourth 
harmonic, their dependence upon K was found 
to be too small to detect. Hence it has been 
concluded that the cross terms in the. energy 
expression are too small to be considered in this 
treatment. The constants and were 
then determined by fitting an equation of the 
form of (8), without the K dependence, to the 
observed combinations by the method of least 
squares. The resulting constants, evaluated from 
the data from the fourth harmonic where over- 
lapping was least serious, were: 


X 190° = 0.3996 Xa0°= —8.1XK10-§ 


The variation of Xj9” with vibrational state was 
best determined by the following relations: 


°Rx(J—1)+°P 


(10) 


These combinations were made in the two 
harmonics, and no systematic deviations from 
linearity in M could be detected. The average 
value of the constant a9! was 0.000610 cm-'. The 
change in the small quantity Xoo” with vibra- 
tional state was beyond the limit of experimental 
determination. 

The rotational constants may now be collected : 


19.946 —0.905v;— 

X 02” = — 0.05616+0.0009040;+ ---. 

X 19” = 0.3996 —0.0006100,;—---. 

Xo” = —0.0000081+---. 
Line intensities 

The relative intensities of the lines comprising 

the parallel component, and the relative in- 
tensities of the lines comprising the perpendicular 
component of the third harmonic band may now 
be calculated to the symmetric rotator approxi- 
mation with the aid of the matrix elements of the 
electric moment given by Dennison for symmetric 
rotator wave functions.’ In order to determine 
the intensity ratio of the two separate com- 


7D. M. Dennison, Rev. Mod. Phys. 3, 313 (1931). 
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ponents, however, one must appeal to the experi- 
mental intensities, which are very difficult to 
measure in this region. Fortunately this task was 
made easy by the fact that two sharp lines of 
moderate intensity, ?Q. and ”Q3, which are close 
together in the third harmonic and in a region 
quite free from other lines of appreciable in- 
tensity, are found to have equal intensities, both 
by visual examination of prints of the spectro- 
grams and also by inspection of microphotometer 
traces of the plates. This useful fact establishes 
the ratio of the parallel and perpendicular com- 
ponents in the third harmonic. Fig. 1 shows these 
symmetric rotator intensities for the individual 
coarse structure lines, and the envelopes of the 
fine structure intensities. It is interesting to note 
that the calculated intensities of the ?P, #P, ?R, 
and ¥R branches are too small to be represented 
on the scale of Fig. 1; no lines of these branches 
were found. Because of the inherent difficulties of 
intensity measurements in this region it is quite 
impossible to make any quantitative comparison 
of theoretical and observed intensities; suffice it 
to say that a summation of the intensities of 
Fig. 1 produces a band envelope whose similarity 
to a microphotometer trace of the third harmonic 
band is amazingly close, except for two ex- 
ceptional lines which will be mentioned shortly. 
A similar comparison in the fourth harmonic 
yields comparable results. 

This treatment also enables one to calculate 
the ratios of the vibrational matrix elements of 
the parallel and perpendicular components of the 
electric moment. One finds: 


(My) 


It is of interest to notice that there is no evident 
relationship between these observed ratios and 
those which one would calculate on the assump- 
tion that the change of electric moment takes 
place along the N—H bond, for such an assump- 
tion would predict a ratio of 6.7, nearly inde- 
pendent of 2.8 

Up to this point, no departures from the 


8 This assumption has been used by Badger and Bauer, 
J. Chem. Phys. 5, 854 (1937). It is probable that the 
appearance of the third harmonic of the N—H stretching 
frequency in ethylene imine is a result of the failure of 
this assumption, rather than an indication of an abnormally 
small angle between the N—H bond and the plane of the 
ring, as was first supposed ; see the reference of footnote 6. 
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symmetric rotator intensities have been re- 
marked ; such departures do, however, exist, and 
are found precisely where they should be ex- 
pected. The @P, and Ry branches, whose higher 
lines should feel the effects of asymmetry most 
markedly, are completely submerged by stronger 
branches. The first line of the °Q branch, °Q,, is 
already very sharp and in the correct position. 
The lines ?Qp and ’Q,, however, both of which 
obtain their greatest intensity from transitions 
involving levels with K=0 for J values around 
16, are noticeably weak and diffuse, the former 
merely giving anomalous intensity to two fine 
structure lines, and the latter being almost 
indistinguishable. The estimated positions of 
these lines, however, do not seem to deviate 
greatly from the symmetric rotator predictions, 
though they were felt unsuitable to use in 
determining the rotational constants. The next 
lines in these branches have become very sharp, 
and no further deviations from the behavior of a 
symmetric rotator were found in the third 
harmonic band. None were observed at all in the 
fourth harmonic where the perpendicular com- 
ponent is absent. 


VIBRATIONAL CONSTANTS 


The vibrational energy expression for hydrazoic 
acid will have the form: 


G(v) = Vi(oi td) ++ + Yur(oit3)? 


+---+-++cross terms. 


(11) 


Determination of the centers of the third and 
fourth harmonics of the N—H frequency of 
quantum number » will make possible the 
evaluation of the vibrational constants Y; and 
Y,;. The band centers are given by: 


Go" =G(v) —G(0) = + X01") 


+K4(X 02°—Xo2") +--+. (12) 


With the aid of the above expression and the 
values of the rotational constants already de- 
termined, values of Go” were calculated from 
each of the °Q lines in the third and fourth 
harmonic bands, giving as best values: 


Go = 9547.31 Got=12,412.19 cm. 


From these the vibrational constants are found 


“azoic 


(11) 
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Fic. 2. The assumed model of the hydrazoic acid mole- 
cule; the size parameters have the relative values given 
by this investigation. 


to be: 
Y,= 3497.32 —78.72 


This information permits the center of the 
fundamental band to be calculated as 3339.88 
cm~!. Davies gives 3336.2 cm as the apparent 
center of the parallel structure of the funda- 
mental,® and incomplete results in this laboratory 
give 3335.6 cm-. Since the apparent minimum of 
the parallel component will be found at a 
frequency somewhat lower than that of the 
actual band center, this agreement seems ade- 
quate, and constitutes the justification for the 
vibrational assignments made for the photo- 
graphic bands. 


STRUCTURE OF THE MOLECULE 


In order to interpret the rotational constants 
Xo:° and X49° in terms of effective moments of 
inertia of the ground state of the molecule, it is 
necessary not only to assume the planar struc- 
ture, but also to assume explicitly that I4°+J,° 
=I¢°; the latter assumption will of course be 
only approximately true since the I°’s are effective 
moments of inertia, but for the present purpose 
this is one of the least important errors. With 
this assumption and Eqs. (2) and (3), we may 
obtain: 


Ao=20.346 cm-!; Bo=0.4035 cm™!; 
Co=0.3957 cm-!. 
The effective moments of inertia, expressed in the 


convenient chemical units, i.e., mass on the 


chemical atomic weight scale and distance in A, 
are: 


I,” = 0.8286; Ip” =41.78; Ie” =42.61. . 
*M. M. Davies, Trans. Faraday Soc. 35, 1184 (1939). 


The assumed shape of the hydrazoic acid 
molecule is shown in Fig. 2. Having assumed the 
planarity of the molecule and the linearity of the 
azide group—assumptions in which one has 
great confidence—the molecule still has four 
structural parameters, which are conveniently 
chosen as: 7, the N—H distance; a, the H—N—N 
bond angle; d, the length of the azide group; and 
k=r2/d, which places the central nitrogen atom. 
The moments of inertia can define only two of 
these parameters, so that we must assume values 
for the other two. The following have been 
chosen as the most suitable assumptions: 

(1) The parameter ~ has been chosen as 
1.012A, on the basis of the 1.016A N —H distance 
in the ground state of ammonia. It is felt a 
decrease of 0.004A should be made to care for the 
absence of the H—H repulsions in HN; and for 
the fact that the nitrogen atom involved in N—H 
bonding in HN; is partially double-bonded in the 
azide group. 

(2) The parameter k has been chosen as 0.524. 
If one assumes the structures 


. 
‘N=N=N‘° “N=N=N 
H H 
(a) (b) 


to be the only ones contributing appreciably to 
the ground state of hydrazoic acid, then one can 
make an approximate calculation of the depend- 
ence of both d and k upon x, the ‘‘hybridization 
ratio,” which is a quantity ranging from zero to 
one which measures the fractional contribution of 
state (a) to the ground state.!° Fig. 3 shows the 
result of this calculation. Unfortunately we do 
not know the value of x which will prevail in the 
actual molecule, but one sees from Fig. 3 that the 
azide group can achieve greatest compactness, 
and hence presumably the greatest binding 
energy by choosing an x of about 0.63; the 
corresponding k is 0.524. 

One may now determine the values of a and d 
from the moments of inertia and the assumed 
values for k and r. One finds: 


a=110° 52’; d=2.369A. 


1° This method is described and applied to a somewhat 
different case on page 180, L. Pauling, The Nature of the 


Chemical Bond (Cornell University Press, 1939). 
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Though one might expect the value of d to be 
less accurately given by the curves of Fig. 3 than 
is k, it is interesting to note that the minimum 
value of d given by the simple theory is 2.371A. 
Such close agreement is of course accidental, but 
fortunately neither a nor d are very sensitive to 
reasonable changes of k. Table III shows a 
comparison of the parameters of the azide group 
as determined by this method, as calculated by 
the simple theory, and as determined in other 
compounds containing the azide group. The 
differences between the values obtained by these 
four methods are gratifyingly small, being within 
the limits of error of the diffraction and crystal 
structure determinations. It is consequently felt 
that the values 


Ni—H=1.012A Ne—N3=1.128A 
Ni—Ne=1.241A YH—Ni—N2=110° 52’. 


constitute a reasonably precise description of the 
mean positions of the atoms in the ground state 
of hydrazoic acid. The structure of the azide 
group isin complete agreement with the structural 
theory, and hence tells us nothing new in 
principle. The angle a is of greater interest, since 
it might seem somewhat smaller than one would 
expect on the basis of structures (a) and (b) and a 
hybridization ratio (a)/(b)=0.63. It is felt, 
however, that this is a real effect, and is another 
example of the tendency of the unshared electron 


TABLE III. Parameters of the covalent azide group. 


ELECTRON a. PRESENT SIMPLE 
DIFFRACTION] STRUCTURE Stupy THEORY 
PARAMETER CH3Ns CsNa(Na)3 —N3 
k 0.530 0.532 — 0.524 
d 2.34 2.37 2.369 2.371 
a 120° — 110° 52’ — 


EUGENE H. EYSTER 


“\ 
k 


2 6 init 10 


2.390 


Fic. 3. Dependence of the parameters d and k=r2/d upon 
the hybridization ratio x, 


pairs on the nitrogen atom to prefer the 2s 
orbital, and hence to give to the bonding orbitals 
somewhat more than the usual 2 character, a 
procedure which leads to smaller angles. 

Finally, this structure makes it easy to under- 
stand the large value of the constant ap,', which 
is responsible for the unusual nature of the 
photographic bands. In fact this constant can be 
roughly predicted, for with this molecular struc- 
ture, one sees that ao:' should be given approxi- 
mately by the value of a, for a diatomic molecule 
with B,=20.8, w.=3497, and w.x.=78.7. The 
Morse approximation, which usually gives too 
small values for a., then gives a-=0.705, which is 
certainly of the right order of magnitude. 

In conclusion, the author is pleased to acknowl- 
edge his special indebtedness to the National 
Research Council, and to express his gratitude to 
Professor Richard M. Badger of the California 
Institute of Technology, in whose laboratory the 
experimental part of this investigation was com- 
pleted, for his helpful advice and for his continued 
interest in this problem. 
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Theory indicates that the normal modes of vibration of 
chain molecules fall into two classes; end vibrations in 
which the amplitude falls off exponentially from the ends 
towards the middle, and chain vibrations in which it 
varies harmonically as in a standing wave. The transition 
from one class to the other takes place via chain vibrations 
with a small number of nodes of amplitude along the chain. 
For a long chain we expect only those chain vibrations to 
show up in the infra-red spectrum for which the phase of 
the electric moment is constant along the chain. ; 

These considerations have been applied to predict the 
vibration spectra of the straight chain paraffins, assumed 
to consist of a zigzag chain of carbon atoms each carrying 


two hydrogen atoms in a plane perpendicular to the chain 
and alternately above and below it. 

The infra-red spectrum of undecane has been measured 
from 1 to 15 on a prism spectrometer, the positions of 31 
bands being determined. Accepted values of the interatomic 
distances and angles were assumed and force constants 
were adjusted to fit four of the fundamental frequencies 
to observed bands. Bands were observed near all the other 
fundamental frequencies predicted certainly to occur (as 
chain vibrations) and near one out of three that should 
possibly occur (as end vibrations near chain vibrations 
with few nodes). All but three of the observed bands could 
be accounted for as fundamentals, harmonics, or simple 
combination frequencies. 


I. INTRODUCTION 


HE spectrum of undecane and those of 
similar compounds have been measured by 
several investigators. The first investigation was 
carried out by Coblentz.!:? Later Lambert and 
Lecomte’: measured the spectrum in the region 
between 6u and 16y and Liddell and Kaspar® the 
region from 0.5 to 1.84. More recently Barca- 
Galateanu®:7 has measured the spectrum of 
dodecane, the next member of the series, from 2 
to 15y, and the results are very similar to those 
for undecane. The work of Lambert and Lecomte 
on the other members of this series also show: 
this similarity. 
Various attempts have been made to explain 
the observed spectra. The first was by Ellis* and 
was purely empirical. He interpreted part of the 
bands as being harmonic members of a series 
whose fundamental at 3.34 was due to a C—H 
vibration. Other bands were attributed to a 


' Coblentz, Astrophys. J. 20, 207 (1903). 

* Coblentz, Phys. Rev. 16, 385 (1903). 

* Lecomte, Comptes rendus 183, 27 (1926). 

* Lambert and Lecomte, Ann. de physique 18, 329 (1932). 

5 Liddell and Kaspar, Nat. Bur. Stand. J. Research 11, 
599 (1933). 

° Barca-Galateanu, Bull. Soc. Roum. de Phys. 37-67, 
19 (1936). 

? Barca-Galateanu, Bull. Soc. Roum. de Phys. 38-70, 
104 (1938). 

‘Ellis, Phys. Rev. 27, 298 (1926). 


C—C vibration with an unobserved fundamental 
at 28u. Other observers have applied these series 
with fair success to various members of the 
methane series. However, without making use of 
improbable combination frequencies, they could 
not account for all of the observed bands, nor 
could they explain the intensities. Later workers 
have attempted to explain the observed spectra 
on the basis of a model of the molecule. Lennart 
Simons® assumed a model composed only of 
carbon atoms arranged in a straight line and 
obtained expressions for the normal frequencies 
of vibration. Kassel’® has obtained the normal 
frequencies of pentane taking account of the 
carbon atoms and assuming an angle of 110° 
between the C—C bonds. Kirkwood™ more 
recently has worked out the normal frequencies 
for a model of any length taking into account the 
angle between the C—C bonds. All of these 
workers have, however, neglected altogether the 
effect of the hydrogen atoms, which must cer- 
tainly be important. Accordingly it was thought 
worthwhile to investigate the whole problem 
taking into account all of these factors and thus 
to set up a fairly complete theory of the vibration 


spectra of these long chain compounds. 


* Lennart Simons, Soc. Scient. Fenn. Com. Phys. Math. 
8-16, 1 (1935). 

10 Kassel, J. Chem. Phys. 3, 326 (1935). 

" Kirkwood, J. Chem. Phys. 7, 506 (1939). 
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II. THe THEORY OF THE VIBRATIONS OF CHAINS 
OF SIMILAR DYNAMICAL SYSTEMS” 


Consider a chain of m links each described by 7 
coordinates, yur, #=1, 2, ---r, h=1, 2, 
completed at the ends before h=1 and after h=m 
by systems described by g coordinates, x, =1, 
2, and s coordinates, 2, w=1, 2, 
chosen so that the kinetic energy has the form 


The total number of degrees of freedom is thus 
qtmr-+s. 

Suppose that the potential energy contains no 
product terms involving x, and y., if h>1; yu, and 
Yor if >1; or yun and if h<m+1-—1; 
where / is a small integer usually 1 or 2. Suppose 
further that the coefficients of the terms in 
Yur¥on+k (which differ from zero only for |k| +2) 
are independent of h. 

This form does not represent a finite chain 
molecule exactly because of progressive small 
changes in the internal and relative equilibrium 
distances of the links as the ends are approached. 
On the other hand, the coordinates may be defined 
differently for these links so as partly to com- 
pensate this. (Indeed for C,Hen+2 with r=9 we 
may perhaps take m=n, g=s=3; x1, x2, and xs, 
being combinations of coordinates of the end 
CH; group, chosen with yu, yi2, so that 
these are coupled to ye: etc. by nearly the same 
potential energy terms as these with ys: etc. ; but 
it is safer to take whole CH; groups as end 
groups.) 

If the system is vibrating in a normal mode 
with frequency w we may write 


(2) 


We then regard the 0,,, as linear combinations of 
terms of the form 


b™. (3) 


The equations of motion for yu,, A=/+1, 1+2, 
---m—l, give a set of r homogeneous linear 
equations in },, linear in w*, and of degree / in 


® Routh, Advanced Rigid Dynamics (Macmillan, 1930), 
p. 287, where some references are given. The first consider- 
ation of = of this kind seems to be due to Lagrange; 
chains of gyrostats were considered by Kelvin; the general 
theory is not difficult to set up and presents analogies with 
the theory of electronic states of linear crystals. 


WHITCOMB, NIELSEN AND THOMAS 


cos\, the term in w not containing >. The 
equation 

D=0, (4) 
formed by equating to zero the determinant of 
the coefficients of }, in these, is an equation of 
degree 7 in w’, the coefficient of w** being in 
general of degree /(r—s) in cos \. This equation 
connects the simultaneously possible values of w 
and X for the oscillations of an infinite chain with 
phase difference \ from link to link. 

For a given value of w there are in general / 
values of cos \ and sets of ratios of b,. If we take 
for b., linear combinations of the 2/r terms 
b,e*""*, and substitute these, together with a, and 
Cw, in the end conditions, the equations of 
motion, that is, for x2; Yur, h=1, 2, m+1—1, 
m+2—l, ---m; and 3»: we obtain a set of 
q+2ir+s homogeneous linear equations in 
2ir+q+s variables. The determinant formed by 
the coefficients of these equations when equated 
to zero must give the condition that w is a normal 
frequency. This determinant contains only sym- 
metrical functions of the 2/r values of e® for 
given w and can therefore be expressed without 
actually solving the Eq. (4) for cos. So ex- 
pressed it must be equivalent to that obtained by 
eliminating all of ai, bua, and cy», from the 
(original) equations of motion. 

If the r values of w? satisfying D =0 are given as 
functions of \ from \=0 to \=2z7, it may happen 
that some of them vary little over this range and 
lie in bands separated from those covered by 
other roots. If any normal frequency lies in such 
a band, the corresponding normal mode must be 
a linear combination of at least one pair of sets of 
terms b,e*'* with \ real and at most /r—1 pairs 
of sets with \ complex, which must die off 
exponentially from the ends towards the center. 
The mode, for a sufficiently long chain, is then 
given by a sum over the sets with } real 


cos (hA—e), (5) 


modified near the ends to satisfy the end con- 
ditions. The modifications near the ends are 
usually small for vibrations that can be regarded 
as near the internal vibrations of the separate 
groups and of the order of 1 for those that can be 
regarded as relative motions of the groups, but 
may exceptionally be much larger. Such modes 
we shall term chain vibrations. (Fig. 1, A.) 
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' Fic. 1. These diagrams represent a chain molecule (e.g., undecane with bl 
nuclei, looked at perpendicular to its length and paraliel to the plane containing the C nucle 


chain vibration (A= m/4); (B) the last member of a set of m such modes 


; (D) a similar mode of which the frequency happens 
to lie inside a band of chain frequencies. For (A) and (B) the displacements with amplitude varyin 


exponentially along the chain are most important, for (C) and (D) these are small compared with 
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If a normal frequency does not lie in any of 
these bands, the corresponding mode must be 
compounded entirely of sets of terms with \ 
complex falling off exponentially from the ends 
towards the middle. Such modes will be termed 
end vibrations. (Fig. 1, C.) 

In the simpler situations all the modes including 
translations and rotations of the whole chains are 
accounted for as follows: one end vibration with 
frequency near that of each internal vibration of 
an end group by itself, and one for each vibration 
of an end group relative to the chain, g+s modes 
in all; and m modes for each of the 7 values of w* 
satisfying (4). Of the latter all, except, for some 
end conditions, a few, are chain vibrations, a few 
may be end vibrations with only complex values 
of \, the frequencies corresponding to which lie 
outside but near the bands of chain frequencies. 
(Fig. 1, B.) These modes have in the simplest 
cases just 0, 1, ---m—1, nodes of amplitude along 
the chain; and when each involves but one real 
value of \, these differ from one another, except 
near \=0 and \=7,, by approximately 7/m. 

This simple situation may be complicated by 
the overlapping of bands by bands, (for links 
with symmetry, overlapping of bands for modes 
of different symmetry may be immaterial) ; or by 
what would otherwise be frequencies of end 
vibrations lying within bands, giving rise to chain 
vibrations withlargeend modifications. (Fig.1,D.) 

The procedure suggested for a long finite 
chain, then, is first to solve the equation D=0 for 
r values of w? as functions of X, essentially the 
problem of an infinite chain; secondly to solve 
the equations for a semi-infinite chain from each 
end with no amplitudes increasing exponentially 
away from the end, /r+q or Ir+s homogeneous 
linear equations in as many amplitude coeffi- 
cients, these provide first approximations to the 
end vibrations (which may in case of resonance or 


Fic. 2. This diagram illustrates how for a general chain 
molecule the electric moments perpendicular to the axis 
are in phase when \= +4 (in this case 7/5). 


symmetry involve both ends), and, used as end 
conditions, fix more exactly the chain vibrations. 


III. THe Opticat Activity oF CHAIN 
VIBRATIONS 


The displacement from one member to the 

next in a chain of similar and similarly situated 
groups can always be regarded as made up of a 
translation along and a rotation of not more than 
two right angles about an axis with unique 
direction and unique position (unless the rotation 
vanishes). If this screw has sufficient pitch, 
chains of indefinite length should be possible. 
For a long chain molecule with many complete 
rotations from one end to the other, the axis of 
this screw may justly be called the axis of the 
molecule. 
. A chain vibration will give rise to bands 
observable in the infra-red only if the electric 
moments in any direction produced in the suc- 
cessive links are nearly in phase. For electric 
moment parallel to the axis this means \~0; for 
electric moment perpendicular to the axis \= +8, 
where @ is the angle of rotation from one link to 
the next. (Fig. 2.) In the special case of no 
rotation from one group to the next we must have 
<0 for the electric moment parallel to any 
direction in space. 

If chain vibrations are given roughly by 


burn= cos (hA—€), (6) 
x 


the electric moments for vibration parallel to the 
axis are proportional roughly to 


sin 


h=1 sin $d 
Xcos {3(m+1)rA—e}, (7) 


for vibration perpendicular to the axis at angle 
to the initial direction from which rotations are 
measured, to 


cos (AA— €) cos (hO— 


sin 3m(A+ 0) 
= cos [3(m+1)(A+8) 
sin 3(A+8) 
sin }m(A— 0) 
—(n+s)}+ 


Xcos (8) 


(6) 


the 


(7) 


le 7 
are 


(8) 
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Fic. 3. The structure assumed for the straight chain paraffins. 


where as } varies, the energy of vibration remains 
nearly the same, of order of magnitude 


(9) 


so long as there is at least one value of \ that is 
real and not too small. 

As the values of \ depart from nearly the above 
values (0 or +8), the intensity of emission or 
absorption must fall off rapidly according to the 
first factors of the above expressions (7) and (8) 
like that in the diffraction pattern due to a slit. 
If the molecule is symmetrical end to end, the 
intensity is exactly zero for modes with an odd 
number of nodes of electric moment (this applies 
to C,Hen42 whether is even or odd). If but one 
real \ is involved, changing in steps of +/m, the 
last factors have one value, cos y say, for modes 
with an odd number of nodes of amplitude, and 
then sin y for those with an even number of 
nodes. 

If a mode corresponding to a root of D=0 is an 
end vibration, however, the energy of vibration 
for given b,, may be much smaller than (9), 
while, for suitable end conditions, any component 
of electric moment may be appreciable, irre- 
spective of 6. 

To a first approximation the infra-red spec- 
trum of such a molecule can therefore be pre- 
dicted using only the solutions of D=0 for \=0 
for electric moment parallel to the axis and for 
\= +60 for electric moment perpendicular to the 


axis, to obtain the frequencies due to chain 
vibrations. In addition, for suitable end con- 
ditions, end vibrations may occur giving fre- 
quencies near the limits of the bands of solutions 
of D=0, most usually for \=0 or z. Only the 
qg+s modes properly belonging to the end groups 
are then left out of account. 


IV. THE VIBRATIONS OF LONG STRAIGHT CHAIN 
PARAFFINS 


The model assumed for the structure of these 
members of the methane series consists of m 
carbon atoms lying along two parallel lines in a 
plane so that there is an angle of 28 between the 
C—C bonds of length 6. The hydrogen atoms are 
assumed to be grouped in pairs in perpendicular 
planes through the carbon atoms, the angle 
between the C—H bonds of length a at each 
carbon atom being 2a. The end groups of the 
chain carry an additional hydrogen atom which 
is important only in the consideration of end 
conditions. The arrangement of five CHe groups 
near the middle of the chain is shown in Fig. 3. 

We take coordinates X;, Yr, Z,, for the dis- 
placements of the carbon atoms along one of the 
parallel lines, perpendicular to the plane con- 


taining the carbon atoms, and perpendicular to 


these away from the other parallel line, and 
Xhy Vay Zn, and yn’, 2x’, correspondingly for the 
hydrogen atoms of the hth. CH» groups. These 
coordinates are oriented as shown in Fig. 3 in 
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such a way that a rotation of the whole group 
through 180° about the parallel line and a 
translation along a C—C bond lead to the 
corresponding coordinates of an adjacent group. 
Thus the rotation angle @ of Section III is just 7. 

Taking the mass of a carbon atom as M and 
that of a hydrogen atom as m, the kinetic energy 
of the Ath. CHe group is given by 


(10) 
% ren 
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The potential energy is assumed to be that due 
to central forces; a linear combination of squares 
of the elongations of the distances between atoms. 
Gi, Y2, * denote elongations of distances with 
equilibrium values as indicated in parentheses: 


=a), gs(H;,—H,’ = 2a sin a), 
qa(Ci, Chai = b), qe(Ci Crise = 26 sin B), 
sin B), 
together with the analogous interactions of the 
other atoms. Here 
sin? B, 
d?=a? sin? a+ 0? sin? B+ (a cos a+b cos £)?, 
f?=4¢@ sin? a+ 40? sin? B. 


(11) 


Ve 
Ae 
% Aco 


Fic. 4. These diagrams show the nature of the displacements of the atomic nuclei in the chain vibrations of undecane 
(and of other straight chain paraffins) with \=0 or x. For motion perpendicular to the axis, e.g. v1, A=, a portion of the 
molecule is viewed along the axis; for motion parallel to the axis, e.g., v1, 4=0, the molecule is viewed nearly perpendicular 
to the plane containing the carbons, and three consecutive CH: groups are shown in projection. 
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The terms in the potential energy involving 
the coordinates of the Ath. CHe group are taken 
to be 


+ K + K +95") 
15 19 
> > qi? 
i=8 


i=16 


27 31 
+Kiw>d g?+Kiu qi? (12) 
i=20 i=28 

When the elongations are replaced by their 
values in terms of the coordinates, the equations 
of motion of the Ath. CHe group, where h may 
take any value from 3 to n—2, can be written 
down. 

We regard the motion in a normal mode with 
frequency w as a combination of terms in which 
successive groups oscillate with phase-difference 


—R+Aun, Aj», —1A};3, —Anu, tA is, 
An, —2r+A2, 0, tA 95, 
tA 51, tA 32, —R+A3;, Ax, As, 
tA ai, 0, A4s, —2r+Au, As, 

—1A 52, Ass, Asa, —2r+Ass, 
0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 


d and so substitute in these equations of motion 


We thus obtain a set of nine homogeneous linear 
equations in X, ---, 2’, involving w, A, and the 
angles, dimensions, and force constants. These 
will have nontrivial solutions if the determinant 
formed by the coefficients vanishes. This gives 
the Eq. (4) of the above general case. This ninth- 
order determinant may be simplified since 
X, x+x’, Z, 2+2’, and y—y’, form a set of 
coordinates unchanged, Y, y+y’, z—2’, and 
x—x’, a set changed in sign, by reflection in the 
plane containing the carbon atoms. Accordingly a 
new determinant is formed taking rows and 
columns the proper sums and differences of those 
of the former to correspond to these new coordi- 
nates in the above order, giving 


0, 0, 0, 0 
0, 0, 0 
0, 0, 0, 0 
0, 0, 0, 0 
0, 0, 0, 0 =0, (14) 
—R+Ac, Ag, Ags, — Aes 
Az, —2r+Az, Ax, tA zy 
Ase, As7, 0 
tA —1Agz, 0, —2r+A 


where R=}Mo’, r=}mw*, Aij=Aji, and with the abbreviations k=(b/d) sin 8, 1=(a/d) sin a, 
n=(a/d) cos a+(b/d) cos 8, p=(2b/c) sinB, q=(a/c) sina, s=(a/c) cosa, t=(2b/f) sin B, and 


u=(2a/f) sin a, we have 


A u=K, sin? B(1 — cos A) +2K sk? +2K p?+Ke(1 — cos 
2K,(1 — cos — cos 24) +2Ksk?+ 2K wp", 
A33=K, cos? a+K, cos? B(1+cos A) +2Kgn?+ 2K ys", 
Ass=A Ky cos? a+2Kgn?+2K ys’, 

A;;=Ki sin? a+2K3+2K,,u?(1+ cos 2d) + 2K gl? + 2K 109’, 
Ag=Ki sin? at+2Kl?+2K 

Az;=K, sin® 2h), 
Ayy=2K,(1 — cos 2d) 2d) + 2K sk?+ 2K 
Ay= 2K sk? cos \—2K 0p? cos 

A,3;3=Kz,sin cos B sin i, 

A4=2Kskn sin \+2Kyiops sin 24, 

sin \+2Kiopq sin 24, 

2Kskn sin sin 2x, 

sin 2X, 

—K, cos? a+2Kgn? cos 2K is? cos 2X, 
A¢s=A3;=K; sin a cos a—2Ksln cos A+2Ki0gs cos 2X, 
—K, sin @ COs a—2K,ln—2K 

Ag= —K, sin? \+2K;/? cos A— 2K cos 2x, 

Ago= 2K skl sin \—2Ki0pq sin 2x, 

sin X. 
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Fic. 5. The infra-red absorption spectrum of Ci:Hox. 


Thus the ninth-order determinant splits up for 
all values of \ into one fifth-order determinant 
and one fourth-order determinant. For the values 
4=0 and \=z7, which are of special interest for 
the infra-red spectrum, the determinants split up 
further into one first-order, one second-order, and 
two third-order determinants. The corresponding 
normal coordinates have then the forms indicated 
in the second and last columns of Table II (see 
also Fig. 4). Those containing x coordinates for 
<0 will correspond to optically active chain 
vibrations giving rise to parallel type bands, 
those containing y and z coordinates for \~ 7 to 
perpendicular type bands. 

There may also be end vibrations for fre- 
quencies near those found from (14) most 
probably for \~0, possibly for \~7, or since 
1=2, possibly even for other values of \, which 
could give rise either to parallel or to perpen- 
_ dicular type bands. The two zero frequencies for 
A= correspond to translations of the whole 
molecule along and rotation about the y and z 
directions; the two for \=0 to translation along 
and rotation about the axis (x direction). 


V. COMPARISON OF THE THEORY WITH THE 
EXPERIMENTAL RESULTS FOR 


In order to test the above theory the infra-red 
spectrum of undecane, C,,Ha, has been meas- 


ured in the region from 1 to 154. Undecane was 
selected because it has a chain long enough (nine 
CHe groups and two CHs groups) to make 
certain of its approximating an infinite chain, 
while it is still light enough to give a sufficient 
pressure of gas at fairly low temperatures. 

The liquid undecane was vaporized by means 
of an oven surrounding the absorption cell. The 
spectrum was then measured on a Wadsworth- 
Littrow NaCl prism spectrometer of the type 
described by Strong and Randall.” The results 
are shown in Fig. 5, where the percentage of 


TABLE I. The infra-red absorption spectrum of undecane. 


WAVE- WAVE- 
LENGTH |FREQUENCY LENGTH |FREQUENCY 
B I cm! I 

1.11 9000 2 5.815 1720 2 
1.183 8453 2 6.188 | 1616 1 
1.446 6915 1 6.447 1551 3 
1.780 5618 2 6.835 1463 9 
1.950 5128 1 7.297 1370 8 
2.305 4338 6 7.647 1308 5 
2.410 4149 1 7.901 1266 3 
2.668 3748 3 8.446 | 1184 5 
3.397 2944 2 8.806 1136 3 
3.600 2778 10 9.378 1066 8 
4.310 2320 6 10.508 951.4 4 
4.612 2168 2 11.023 907.2 8 
4.875 2051 3 12.425 805.5 3 
5.334 1875 3 13.310 751.3 6 
5.535 1807 5 13.85 722.0 8 
5.712 1751 3 


8 Strong and Randall, Rev. Sci. Inst. 2, 585 (1931). 
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absorption is plotted against the wave-length and 
frequency. The wave-lengths, frequencies, and 
estimated intensities of the vibration bands 
observed are given in Table I. 

An attempt was made to measure the spectrum 
at a higher dispersion using a prism-grating 
spectrometer. However the complexities intro- 
duced by the rotational fine structure and by the 
overlapping of bands prevented this work from 
being of much value in the determination of the 
envelope structure. After measuring one of the 
four regions this work was abandoned. 

In order to apply the theory to these results 
certain simplifications were made. It was as- 
sumed that the potential energy constants Ke, 
Ky, and Ky, are negligible and that Ks and Kyo 
are equal. The constants k, 1, n, p, q, s, t, and u, 
were found by taking @ and £6 each half the 
tetrahedral angle, i.e. a=8=55°. The C-H dis- 
tance, a, was taken to be 1.15A™ and the C-C 
distance, b, to be 1.56A." This gives the following 
values: 


c=2.8005A, d=2.2241A, f=3.1786A, 
k=0.57551, 1=0.42356, n=0.69957, p=0.91218, 
q=0.33567, s=0.23503, t=0.80538, u=0.59274. 


Approximate values for K,, K4, and K3, were 
found from the bands at 3000 cm!, 990 cm-', and 


1600 cm~', which are found in all compounds of 


4 Kistiakowsky, Phys. Rev. 45, 4 (1934). 
 Murison, Phil. Mag. 17, 201 (1934). 


this type. The constant K; is found by using the 
formula 


(16) 


where y; is the reduced mass of the vibrating 
group. These constants were used to make an 
assignment of the frequencies and then accurate 
values were found as follows, from (14) with A=0 
or A\=7, 


Ks and Ky from the parallel band at 1463 cm~', 

0.46557 K 10° dynes ‘cm, 
K, from the band at 2320 cm™, 

2.19551 K 10° dynes cm, 
K; from the band at 2778 cm™, 

1.04932 K 10° dynes ‘cm, 
K, from the band at 907 cm", 

2.87614 KX 10° dynes ‘cm. 


Using these constants the values of the fre- 
quency were obtained for the cases \=0°, 5°, 
45°, 90°, 135°, 175°, and 180°. These values are 
given in Table II and are plotted in Fig. 6. There 
is little variation of the frequency with the phase 
except for those types of which the frequencies 
vanish for phases 0 and 7. 

We expect chain vibrations with the following 
frequencies to be optically active and to give rise 
to fundamental bands: 


frequency incm 385 1106 1354 1463 2320 2778 )° 


TABLE II. The roots of D=0. 


0° 0° 5° 45° 90° 135° | 178° 180° 180° 

X+(x+x’) 0 75.72 | 456.80 | 423.10} 229.08 0 0 Z+(s+s') 

v4 — Z+54'(2+2') | va" — 
+es'(y—y’) 907.02 | 907.03 | 953.70 | 889.70 |1076.59 |1101.86 1105.65 | +e 


V5 
+es'(y—y’) |}1038.32 |1038.53 |1062.10 


1066.95 |1115.97 1112.17 1106.71 | as 


|| 1463.00 |1461.02 |1360.60 


1473.18 |1398.94 |1413.89 |1416.03 


¥3 


—€3'(y—y’) ||2751.94 |2751.53 |2756.48 


(s+s")— 142} 


2755.96 |2770.44 |2778.22 |2778.00 || 
—as'(s—2’) 0 19.32] 172.25 | 267.90} 151.61 0 0 ¥+(y+y") 
+o6'(z—2’) |} 604.20 | 603.13 | 585.65 | 513.50} 399.23} 385.30) 385.21 
V9 x—x’ 1354.46 |1354.57 |1356.06 |1359.67 (1371.66 |1354.88 |1354.46 || x— x’ 


+o7'(2—2") ||2272.83 |2272.78 |2272.20 


2276.30 |2302.52 |2319.75 | 2320.00)| 
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while end vibrations with the following fre- 
quencies may well exist, be optically active, and 
give fundamental bands: 


designation veloc) v7(oc) ve(oo) 
frequency incm=! 604 907 1038 2273 2752 J 


The number suffix in the above designations 
refers to the row of Table II and to Fig. 4, 
the first symbol in the parentheses gives the 
value of \, either 0 or 7, the second the type of 
band expected, either parallel, 7, or perpen- 
dicular, o. 

As stated above the frequencies at 907, 1463, 
2320, and 2778 cm™!, were used to fix the force 
constants; the observed bands at 1066 and 1370 
cm~! are assumed to be those predicted at 1038 or 
1106 and 1354 cm~. The envelope structure 
shows clearly that those at 1370 and 1463 cm 
are of parallel type. The best comparison of the 
theory with experiment is found, however, in the 


consideration of the harmonic and combination 
bands. These bands are given in Table III and it 
is seen that, with the exception of three weak 
bands, all the bands are accounted for within five 
percent. In some cases more than one combina- 
tion band seems to fit the observed frequency ; 
wherever possible on the basis of intensity con- 
siderations, the most important band is indicated 
by an asterisk. 


VI. CoNCLUSIONS 


The theory of the vibrations of long chain 
compounds as developed here seems capable of 
accounting completely for the infra-red spectrum 
of undecane, although the problem is incompletely 
solved because the end conditions have not been 


TABLE III. Observed and calculated frequencies. 


FREQUENCY, 


Oss. Can. Av DESIGNATION Oss. Cat. Av DESIGNATION 
9000 1875 1876 — 1 | 
8433 8334 +99 | 3y3(x0) *1848 +27 | ve(ro)+v2(or) 
6915 6960 —45 | 1807 *1814 — 7 | 2vs(oo) 

6910 + 5 | 1812 — 5 | 3ve(oc) 
5618 5556 +62 808 — 1 | 2v6(ro)+vs(o0) 


2v3(3a) 1 

5128 5098 +30 | |} 1751 17389 +12] 
4338 *4389 —51 | 3v2(or) 1720 1710 
4280 +58 | 2ve(or)+yr9(or) 1616 1642 —26 

4149 4171 —22 | 2vo(or)+ve(or) || 1551 1540 +11 | 
*4132 +417 1463 1491 —28 
3748 3816 —68 *1463 — vo(or 
*3783 —35 | || 1870 1874 — 4 | 

3685 +63 | || 1808 1354 —46 | (or) 
2944 *2926 +18 | 2ve(or) 1266 1292 —26 | 

2924 +20] 1184 1208 —24 | 2vé(oc) 

2778 2817 —39 | vo(or)+v2(or) 1155 —19 | 3v6(xo) 


2778 — 066 *1106 —40| vs(x0) 
2721 +57 | 1038 +28] 
2320 —50 | 951 989 —38 | 


907 97 — vs(oo) 

2314 + 6 | 805 
2168 2199 —31 | 751 770 —19 | 2ve(xo) 
*2144 +24 722 
2051 2076 —25 | 2v:(oc) 604 ve(oo) 
2067 —16| 385 
2013 +38 | 


considered explicitly. The theory may be applied 
just as easily to the other straight chain hydro- 
carbons of the methane series. It is possible that 
it may also be applicable to the derivatives of the 


‘series, with modifications to take account of the 


added groups. A similar theory should explain 
the spectra of the cyclic hydrocarbons. The 
agreement of the observed and calculated results 
for undecane verifies the theory and substanti- 
ates the choice of dimensions and constants. 
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Displacements, estimated intensities, and depolarization factors are listed for the main lines 
of diethyl ether, di-n-propyl ether, di-isopropyl ether, di-n-butyl ether, di-isobuty] ether, di-n- 
amyl ether, and di-isoamyl ether. All the frequencies in the 2900 cm region are polarized, 
except the one near 2970 cm™ which is, without exception, highly depolarized. 


INTRODUCTION 


Ww the exception of dimethyl and diethyl 
ethers, the Raman spectra of the simple 
ethers have not been subject to thorough investi- 
gation. High! reported data for the propyl and 
amyl ethers, but the data are obviously incom- 
plete, probably due to the fact that he was 
interested in measuring the main lines for many 
related compounds rather than in a detailed 
study of any particular one. Since High appar- 
ently used unfiltered mercury light for excitation 
of the spectra, there also existed the possibility 
of inaccuracy in the data due to incorrect 
assignment of the exciting frequency. More 
recently, Cleveland and Murray? have published 
results obtained for di-n-butyl ether with a 
spectrograph having a rather low dispersion 
(40A/mm at 4500A). At the time the present 
investigation was begun, depolarization data 
were available only for dimethyl and diethyl 
ethers. 

In view of this situation, it appeared worth 
while to undertake a systematic investigation of 
the simple ethyl, propyl, butyl and amyl ethers 
with the purpose of obtaining not only the 
displacements and estimated intensities, but also 
depolarization factors for the principal lines. 


EXPERIMENTAL PROCEDURE 


Light from a cylindrical Pyrex mercury arc* 
passes vertically upward through a filter solution 
contained in a glass tube (Wood’s arrangement). 


jm Now at Armour Institute of Technology, Chicago 
inois. 

1M. E. High, Phys. Rev. 38, 1837 (1931). 

*Forrest F. Cleveland and M. J. Murray, J. Chem. 
Phys. 5, 752 (1937). 

*Forrest F. Cleveland and M. J. Murray, Am. Phys. 
Teacher 5, 271 (1937). 


The filter tube acts as a cylindrical lens and 
concentrates the light in a nearly parallel beam 
along the axis of the horizontal Raman tube. 
A plane mirror placed above the tube sends the 
light back through the liquid, thus nearly 
doubling the intensity of illumination. The wave- 
lengths are measured by use of an iron com- 
parison spectrum. The method used in the 
depolarization measurements and the charac- 
teristics of the spectrograph have been described 
elsewhere.‘ 

The filters used were a solution of iodine in 
carbon tetrachloride when excitation was by 
Hg 4047A and an alcoholic solution of p-nitro- 
toluene and rhodamine 5 GDN Extra when 
excitation was by Hg 4358A. During the course 
of the investigation, it was found that the 
excellence of the filter used for isolating Hg 4358A 
made excitation by Hg 4047A superfluous. Conse- 
quently, excitation by Hg 4358A only was used 
for diethyl and di-isobutyl ethers. 


TABLE I. Raman spectrum of diethyl ether. 


PRESENT TABLES PRESENT TABLES 
RESULTS ANNUELLES* RESULTS ANNUELLES* 
Av Av Ay Ay 
| I | p | p p | I | p 
373 |} | 375 | — 1272 2b |0.8| 1275 | 2b/0.90 
439 |6 |04; 440 |5 |0.48| 1452 6b 1453 | 5 |0.89 
496 |}4 |—] 497 |1 | —] 1465 1 09}; — — 
[544] |—| — | 1480 —|— 
830 | 2 —j|— | 2696 —| 2694/1 
842 | 4 843 |8 |0.29| 2724 2 —|— 
927 |4 |0.9] 928 |36/0.61| 2774 1 
1035 | 1b | — | 1033 |3 | — | 2806 4 0.3) 2805 | 5 | — 
1074 |} 1 |0.6) 1080 |0 | — | 2866 | 106 |0.3| 2867 | 6b/|0.30 
1122 |}4 |08; — |—j| — 2898 3 — 
1152 | 46 |0.4/ 1150 |1 | — | 2935 | 106 |0.2| 2932 |10 |0.20 
2979 8 0.9} 2980 |10 | 0.89 


*M. Magat, Tables Annuelles de Constantes et Données Numériques, 
Effet (Hermann & Cie., 1937), Volume XII, Chapter 26, 
page 3/. 


‘Forrest F. Cleveland and M. J. Murray, J. Chem. 
Phys. 7, 396 (1939). 
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All the ethers studied were carefully frac- 
tionated in a heated column using a high reflux 
ratio. Details regarding the samples used are as 
follows: 


Diethyl ether. Dried over sodium. Boiling point, 
34.2°-34.3°. 

Di-n-propyl ether. Eastman product. Boiling 
point, 89.4°-89.6°. 

Di-isopropyl ether. Eastman product. Boiling 
point, 67.5°-67.7°. Distilled over sodium. 

Di-n-butyl ether. Eastman product. Boiling point 
142°. Distilled over sodium. 

Di-isobutyl ether. Synthesized from sodium iso- 
butylate and isobutyl bromide. Boiling point 
121°. (Product gave no action with sodium 
and showed a negative Beilstein test for 
halogen.) 


TABLE II. Raman spectra of propyl ethers. 


Di-NorMAL- Dr-Iso- 

PRESENT PRESENT HiGu's 

RESULTS VALUES RESULTS VALUES 
Av Av Av Av 
(cm™) I p J | (em™)| J p 
323 2 —!| — |—| 306; 3 10.5) — |— 
335 2 —!| — |— 
408 | 4 0.65 — |— 
437. 1 |— 
478 | 2 —| — |—| 492; 4 — |— 
852.| 2 — |}—] 794] 6 |0.2| 791] 6 
870 | 2 0.6 874} 2] 848] 1 | 0.8} 848 | 2 
890 | 3 — 
904 | 2 — |}—| 903] 4 |0.5} 903 | 5 
940 | 26 |o8| — |—| 930] 2 93011 
989 1 
1027 2b 10.9) — |— 
1055 1b — |— 
— |—] 1091*| 2 
1110 1b |—| — 1126] — | — 
1149 3 0.6} — — |— 
1243 2 — |—} 1188 — |— 
1282 3 0.8} 1282 | 3 | 1328] 1 10.7) — |— 
1298 1 —| — 1384] 1 10.7) — | — 

1460 | 6 }| 0.7; 1453 | 8& 
1462 | | 14461 4 | 2 
1482 4 
2737 26 |0.5}) — 1 — |— 
2794 | 3b — !— 
2839 2 —| 
2875 | 8} | 0.3) 2864 | 3 | 2874| 75 | 0.1] 2872 | 9 
2935 | 105 “| 2931 | 3 | 2937) 8 *"| 2929 | 9 
2968 | 6 0.8} — |—J| 2976;10 | 0.9) 2975 | 9 
— | — | —| 30734] 1 
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Di-n-amyl ether. Eastman product. Boiling point, 
68° at 12 mm pressure. 

Di-isoamyl ether. Eastman product. Boiling point, 
62° at 12 mm pressure. 


RESULTS 


The results obtained for diethyl ether are 
given in Table I. The displacements appear in 
column one, estimated intensities in column two, 
and depolarization factors in column three. 
Columns four, five and six give for comparison 
the values listed in Tables Annuelles. Broad lines 
are designated by the letter b and data in regard 
to which there is some question are enclosed in 
brackets. The wider slit used in the depolariza- 
tion measurements prevented resolution of some 
of the lines which were clearly resolved on the 
other plates. Such cases are indicated by joining 
the unresolved lines by a brace. Mercury lines 
(in angstroms) are referred to by letters in the 


TABLE III. Raman spectra of amyl ethers. 


Di-NorRMAL- Dr-Iso- 
PRESENT HiGu's PRESENT Hicn's 
RESULTS VALUES RESULTS VALUES 
v Av Av Ap 
(cm™) I p I |(m 1 p 1 
310 1 — |—] 421 1b}—|] — 
770 4b |—] 764*| 2 ped 3 |0.4 
848 — 1 1 1 
868 3} 0.6] _ |_| g49 
889 1 — |—] 913] 26 908 | 1 
907 1 — 953] 2 950) 3 
981} 1 — 
996; — |— 
1065 2b |0.7; — |—] 1014] 1 — 
nas | 23 (88 3 
114 2 — |— 
1170 | 3 
1298 | 2) |oal | | 1301 | 35 — 
129. 1 3 — 
1306 | 2f | 1292) 4 | 1335] 35 — | — 
1430 | 1448 | 5 
14 
1462 | 4/ | %8| 1447 | 7 | 14631 5} 10.7] 1449 | 9 
1483 3 1485 | 2 
2716 1 0.5) — |—] 2763] 1 |04) — |— 
1 
2871 10} 0.3) 2853 | 5 | 2904 | 9 4] 2906 | 7 
21 — |—|2 
3} — |—| 2953] 0.8) 2958 | 7 
2960 | 4 0.8) 2957 | 5 —_ — |—]3145t| 2 


* 1091(2)kei is probably 902%, 2858k and 1280k, respectively. 
t 3073(1)0g is probably 485k and 29780, respectively. 


* 764 may be due to the p e of di-i yl ether as an impurity. 
+ 3145ko is probably 2956i and 1444m, respectively. 
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Fic. 1, Graphical summary of the Raman spectra of simple ethers, R-O—R. 


tables, as follows: 


e=4358, f=4348, k=4047, i=4078, 
m= 3906, and q=3650. 


The results obtained for the propyl and amyl 
ethers are compared with those reported by 
High! in Tables II and III. Table IV lists the 
results for the butyl ethers. 

A graphical summary of the results for all the 
ethers concerned in this investigation is given in 
Fig. 1. The spectrum of dimethyl ether is in- 
cluded for sake of comparison. In this figure, 
the height of the line is proportional to its 
intensity and the depolarization factors are 
written near the lines. Decimal points have been 


TABLE IV. Raman spectra of butyl ethers. 


D1-NoRMAL- Dr-Iso- D1-NoRMAL- Dr-Iso- 
Av Av Av Av 
| I | p I | p | I | p | | p 
262 | 2 -- 263 | 4 |0.5}) 1120 3 1123 3 0.7 
276 | 2 . 1137 3) }0.8| 1137 3 
294 12 299 1 1156 3 1168 3 0.7 
319 | 2 _ 1181 3 > 
346 |} 1 |— 362 | 2 |— 1228 2 1249 4/108 
398 |} 1 418 | 3 | — 1300 5 |0.8| 1294 2 
450} 1 1339 | 4 |08 
467 1 1439 0.7 1448 7 
820 |; 2 |— 789 — 1459 6 1461 3) 108 
844 | 4 |0.4| 810 — || 1484 3 |0.7| 1478 2 
882 |} 1 |— 836 | 6 2732 2 |0.5| 2718) 3 |04 
979 3} 0.9) 958 | 4 1-7) 3864 | 10) | >| 2866 | 10 | 0.2 
1026 1 1001 1 2876 | 10 (2891); 1 |— 
068 | 4 |0.7 2912 8 03 2905 3 04 
2941 8 2952 1 ? 
2967 6 |0.9| 2963 4 08 
2975 3 


omitted from the depolarization factors in order 
to simplify the diagram. 


DISCUSSION OF RESULTS 


The most intense lines in the Raman spectra 
of the simple ethers are those corresponding to 
the C—H vibration frequencies near 2900 cm-". 
The C—H deformation frequencies near 1450 
cm~ are also relatively strong and are highly 
depolarized. It is interesting to note that no 
lines appear between 500 and 700 cm. This is 
true also for the aliphatic alcohols and for many 
aliphatic hydrocarbons. 

Assuming that, to a first approximation, the 
methyl group in the ethers may be regarded as 
having a pyramidal structure independent of the 
remainder of the molecule, one would expect 
this group to have four frequencies : Two of these 
will be deformation frequencies, in one of which 
(6,,) the change in electric moment is parallel to 
the symmetry axis, while in the other (6,) the 
change in electric moment is perpendicular to 
the symmetry axis. The other two frequencies 
will be vibration frequencies in which the change 
in electric moment is parallel (v,,) and per- 
pendicular (v,), respectively, to the symmetry 
axis. The former will be polarized, the latter 
depolarized. 

* », is believed to be the frequency near 2930 
cm which is polarized, while v, is believed to 
be the frequency near 2970 cm which is de- 
polarized. These assignments are in agreement 
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TABLE V. C—H vibration frequencies. 


—CH;3 —CH2-? 
ETHER 

Diethyl 2935(10b).2 | 2979(8).9 2866(105).3 
Di-n-propyl | 2935(10b).2 | 2968(6).8 | *2865(8).3 
Di-n-butyl *2927(8).3 2967(6).9 | *2870(10).2 
Di-n-amyl *2931(8).4 2960(4).8 | *2860(10).3 
Di-isopropyl | *2928(8).1 2976(10).9 | 2874(7b).1 
Di-isobutyl | *2929(2).4 *2969(3).8 2866(10).2 
Di-isoamyl 2932(5).4 *2959(5).8 2866(10).3 


* Mean of the doublet, which is presumably due to Fermi resonance 
interaction. 


with those made by Crawford for methy]l- 
acetylene® and dimethylacetylene® and by Craw- 
ford, Avery and Linnett’ for ethane. 

The —CH2— group may be considered, 
roughly, as a triangular molecule with inner 
vibrations as follows: A polarized deformation 
frequency (6,,) and two vibration frequencies, 
one strong and polarized (v,,), the other weak 
and depolarized (v,). The former (v,,) may be 
the strong polarized frequency observed at 2866 
cm~. On the other hand, this frequency appears 


> Bryce L. Crawford, Jr., J. Chem. Phys. 7, 140 (1939). 

5 Bryce L. Crawford, Jr., J. Chem. Phys. 7, 555 (1939). 

7 Bryce L. Crawford, Jr., William H. Avery and J. W. 
Linnett, J. Chem. Phys. 6, 682 (1938). 


in dimethyl ether,’ methylacetylene, and di- 
methylacetylene which contain no methylene 
groups. It is assigned by Crawford et al., in these 
cases, however, as an overtone of 1446 cm=. 
Only one depolarized frequency was observed, 
except for di-isobutyl and di-isoamyl ethers 
where the depolarized line appears double. It is 
conceivable that v, for the CHe group may 
coincide with v, for the CH3 group or may be 
close to some other line and, being weak, not be 
observable as a depolarized line. 

The values of the three frequencies, tentatively 
assigned above, are collected in Table V. The 
selection of doublets is somewhat arbitrary, 
especially in the case of the amyl ethers, but 
the constancy of the frequencies seems to justify 
the choice made. It should be noted that lines 
unresolved on the depolarization film and indi- 
cated by the braces in Tables I-IV are not 
necessarily those considered as resonance doublets 
in Table V. 
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The absorption coefficient of nitrogen dioxide in the visible spectrum has been measured at 
pressures up to 70 mm Hg with a 40A spectrophotometer slit. Beer’s law is obeyed in this 
region. The transmission as measured when using a tungsten source, is the same as that of- 
tained with either a high or low pressure mercury arc, except at the 4339-4358A group of 
mercury lines, where the absorption coefficient was found to be 15 percent greater with the 
mercury source. The absorption coefficient curve obtained with either a 40A or 15A slit is 
essentially sinusoidal in appearance, the absorption maxima and minima generally differing 
by about 10 percent in intensity. The maxima are separated by 700+100 cm~. They appear 
to be due, not to the ground state of the molecule, but to groups of excited levels which are 


uniformly separated. 


INTRODUCTION 


HE visible absorption spectrum of nitrogen 

dioxide has been studied in some detail in 
the past by a number of workers.’ Curry and 
Herzberg! have described briefly some features of 
the numerous bands which extend from 8900A to 
the ultraviolet and presented an interpretation of 
some of these bands. Daniels and Holmes? re- 
ported measurements on the absorption coeffi- 
cient at 4030 and 4359A in connection with 
similar measurements for several of the oxides of 
nitrogen. Since we have had occasion to measure 
the absorption coefficient at all wave-lengths in 
the visible spectrum with considerable accuracy, 
it has seemed desirable to record the values 
obtained and also to present several interesting 
features of the spectrum which hitherto do not 
appear to have been described. 


EXPERIMENTAL METHOD 


Nitrogen dioxide was generated by heating 
C.P. lead nitrate. The gas was dried by passing 
over phosphorous pentoxide and condensed in a 
suitable trap. The gas was freed of gaseous 
impurities by repeated freezing and evacuation. 
It was colorless on freezing. 

Equipment for carrying out the work was 
constructed throughout of Pyrex glass. The 
supply of nitrogen dioxide was connected by 


Curry and Herzberg, Nature 131, 842 (1933); L. 
Harris, Proc. Nat. Acad. Sci. 14, 693 (1928); Carwile, 
Astrophys. J. 67, 184 (1928). 
aa and Holmes, J. Am. Chem. Soc. 56, 630 


means of capillary tubing and stopcocks to the 
oil pump, a one-liter storage flask, a source of 
dry carbon dioxide free air, a mercury manometer, 
an oil manometer, and an all-glass photometric 
absorption cell. The manometers, which were 
constructed of capillary tubing, were evacuated 
on one side. 

The entire apparatus was mounted so that it 
could be moved as a unit, thus permitting placing 
of the glass absorption cell in position for 
measurement in a General Electric recording 
spectrophotometer.* The length of the cell was 
measured and found to be 5.01+0.03 cm. 
Photometric measurements were made versus a 
similar cell of the same length filled with air. It 
was possible to check the operation of the 
machine by calibration for percentage trans- 
mission versus a clear glass filter and a gray glass 
filter.* The wave-length calibration was made 
with a mercury arc. 

After generating and evacuating a supply of 
nitrogen dioxide in the trap, air was admitted to 
the storage flask to a definite pressure and 
nitrogen dioxide then allowed to enter from the 
trap to increase the pressure in the flask by some 
desired amount. The storage flask was shut off 
from the rest of the system, the latter evacuated 
and then filled to a measured pressure with the 
mixture. By leaving air in the capillary manome- 


ter, attack of the mercury and oil by nitrogen 


a — and Liebhafsky, Gen. Elec. Rev. 39, 445 
* These filters had been checked against the National 
Bureau of Standards’ filters, 
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Fic. 1. Absorption of light by NO». Cell length 5.01+0.03 
cm tungsten source, wave-length 4850A, slit 40A. 


dioxide was effectively eliminated. The tempera- 
tures of the absorption cell and storage flask were 
recorded. The entire visible range absorption 
curve of the gas in the absorption cell was then 
obtained. The storage flask and system were 
pumped down to a lower pressure and the same 
set of data was obtained again, thus permitting a 
number of measurements to be made at suc- 
cessively lower pressures on one mixture. In 
some of the later series of runs, no air was 
admitted to the system and nitrogen dioxide 
pressures were obtained from readings of the oil 
manometer alone. 

In some cases, the order of addition of the air 
and nitrogen dioxide was reversed to reveal any 
errors of procedure, but no effect was observed. 


‘CALCULATION OF RESULTS AND ERRORS 


Manometer readings were corrected to mm of 
mercury at 0°C.‘ 

The pressure of nitrogen dioxide was calculated 
from the measured pressure, which was equiva- 
lent to the sum NO:+N:0,, by using the 


‘Int. Crit. Tab., Vol. I, p. 69. 
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equilibrium constants given by Verhoek and 
Daniels.> In the dilute concentrations employed, 
N2O, was usually some 85 percent—95 percent 
dissociated, hence errors of +0.5°C in tempera- 
ture measurement were negligible. Although no 
extreme precautions were taken to obtain exact 
equilibrium, it was possible to check the vapor 
pressure measurements of Giauque and Kemp,® 
in the range from 20 to 40 mm with an accuracy 
of +10 percent. In this connection, it should be 
noted that these workers did not report measure- 
ments below 20 mm, but it has been found that 
one may safely extrapolate their curve to 
probably 5 mm, thus supplementing the rather 
uncertain values in the literature.’ 

In the blue and green portions of the visible 
spectrum nitrogen dioxide transmitted only a few 
percent of the light at 70 mm pressure and more 
than 80 percent at less than 4 mm. It was difficult 
to measure the absolute percentage absorption to 
better than 1 percent, thus making the error of 
the order of 1 part in 10 for 10 percent absorption 
and proportionately less at higher absorption 
values. In order to illustrate the randomness of 
the results, the percentage transmission at 4850A 
has been plotted in Fig. 1 on a semi-log scale 
against pressure of NO: (reduced to mm of Hg 
at 0°C) for measurements in six different series of 
runs. It is believed that the results of different 
series of runs are in satisfactory agreement since 
the standard deviation of the logarithm: of the 
transmission for the individual points from the 
line shown is only 0.035.8 

Validity of Beer’s law would require a straight 
line in Fig. 1. This type of curve appears to fit the 
data quite well since the computed estimate of 
standard deviation of the slope is 1.5 percent. 
The straight line may reasonably be extrapolated 
into the region from 90 percent to 100 percent 
transmission in order to establish the transmis- 
sion with more accuracy than obtained by 
drawing the best curve through the experimental 
data in this region.* 


— and Daniels, J. Am. Chem. Soc. 53, 1250 
ue and Kemet < 6, 40 (1938). 
og Crit. Tab., Vol. III, 

P. R. Rider, Statistical "Poinods (John Wiley & Sons, 

1939), Chap. IV. 

* Some points gp 10 mm pressure are omitted for 
purposes of clarity onl - All points have been used in 
estimating errors quoted above. 
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It was also easily shown that Beer’s law is valid 
at all wave-lengths as well as at 4850A. This was 
true with an error less than 5 percent in absorp- 
tion coefficient at all wave-lengths provided the 
absorption was greater than 10 peicent, below 
which figure photometric errors began to increase 
rapidly. 


EXPERIMENTAL RESULTS 


The slope of the transmission curve at 4850A 
(Fig. 1) was used to establish the absorption 
coefficient at this wave-length. The General 
Electric spectrophotometer establishes the trans- 
mission relative to that at any one particular 
wave-length with considerable accuracy, hence 
the transmission values at all wave-lengths were 
easily determined by averaging a large number of 
curves and then comparing the absorption at any 
wave-length for the average curve with the 
value at 4850A. In the relation (percent trans- 
mission /100) =¢= 10-*?! 


p=pressure of NO: corrected from the tempera- 
ture of measurement to 0°C in mm of Hg. 

1=cell length in cm. 

k=absorption coefficient, mm Hg-!Xcm.~1. 


a 


(Loe «) 


= 


ASSORPTION COEFFICIENT 


Lt | 
4000 


«6 
WAVELENGTH ANGSTROMS 


Fic. 2. Probable best value for absorption coefficient of 
NOs. Slit 40A, percent transmission /100=1=10-*?4, p in 
mm Hg: L in cm. k in 1/mm HgXcm. 
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Fic. 3. Slit width effect. " slit 15A, 


Values of logio & are plotted against wave-length 
in Fig. 2. 

The curve shown was obtained with a slit 
width equivalent to approximately 40A. When 
the slit width was decreased to approximately 
15A, the vertical spread between the maxima and 
minima was doubled. 

The appearance of the curve was not affected 
by fairly wide variations in the plotting speed or 
direction of plotting, nor was it changed by 
heating the gas to 85°C. It was also shown that 
the peaks and valleys were not due to fairly 
rapid time fluctuations since the transmission 
did not change with time. 

Since the spectrophotometer can be operated 
so that it plots the percentage transmission with 
fivefold greater scale, it seems desirable to 
present a trace of actual absorption curves 
obtained by this method at two different slit 
widths; namely, about 40 and 15A. Fig. 3 records 
these absorption curves which give the percentage 
transmission versus wave-length in angstroms, 
and illustrates qualitatively how the maxima and 
minima were spread when the slit width was 
decreased. With the narrower slit width used in 
Fig. 3; namely, 15A, further fine structure in the 
bands begins to be evident and with a finer slit a 
large number of absorption bands were observed 
which were in agreement with those evident in 
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TABLE I. Absorption coefficient of NO». 


k (DANIELS 


WAVE-LENGTH, A & Hotmes)* k (THIS WoRK)* 
4050 (mercury source) 0.00936 0.0093 

(tungsten source) _- 0.0093 
4359 (tungsten) _ 0.0085 
4359 (mercury arc) 0.00833 _— 
4339-4359A (mercury arc) _- 0.0098 


* (% transmission/ 100) = 10-*p! (see text). 


photographs published in the literature.® With a 
slit of 40A the absorption coefficient curve is 
that which is shown in Fig. 2 (full line). For 
wider slits the absorption curve would approach 
the average (dotted curve). It is likely that in 
most photochemical or analytical work, this 
curve will be most useful. 

Data available in the literature have been 
concerned chiefly with the absorption coefficients 
as measured with discrete line sources. In view of 
the known fine structure of the spectrum it is 
evident that agreement between different sources 
and slit widths should not be expected. However, 
it seemed of interest to establish the extent of this 

difference. A low pressure mercury source was 
first tried out in the spectrophotometer. Within 
the experimental accuracy of about +1 percent 
in the absorption coefficient, the absorption was 
the same as with the tungsten source at the wave- 
lengths of the 4050, 5460 and 5780A mercury 
lines. On the other hand at the 4339-4358A 
group, the absorption coefficient was 15 percent 
higher with the mercury arc than with the 
tungsten source. The same state of affairs existed 
on using a high pressure arc as source, for which 
source there is a broadening of the mercury lines 
(between which there is a strong continuum).?° 

These results are summarized in Table I. 

The absorption coefficients shown in Fig. 2 
may be utilized, with due regard to the effect of 
slit width, for a continuous source and for a 
mercury arc at its principal visible lines except 
for the group at 4339-4358A. In the latter case, 
it is suggested that a value of k=0.0098, as 
determined in this work, be employed. With the 
4359A mercury line as source, it seems desirable 


® Melvin and Wulf, J. Chem. Phys. 3, 755 (1935). 
10B. T. Barnes, J. Opt. Soc. Am. 27, 83 (1937). 
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to use k=0.0083, the value of Daniels. In cases 
where other sources with discrete spectra are 
employed it would be advisable to determine the 
absorption coefficient since otherwise there 
appears to be no certainty that the curve of Fig. 2 
is applicable. 

The absorption maxima differ on the average 
by 700+100 cm and those for the minima by 
the same amount. The average difference be- 
tween adjacent maxima and minima are 250+50 
cm~. In view of the complexity of each peak or 
valley, it is not desirable to stress intensity 
relationships, but it is notable that for the curve 
of Fig. 2, the absorption coefficient at a peak such 
as A in Fig. 2 is approximately 90 percent of that 
at the valley B over the entire curve. 

The ratio of absorption coefficients for any two 
adjacent absorption maxima in Fig. 2 is greater 
than 0.5. If the uniform difference of 700 cm™ 
between any two bands is due to transitions from 
two adjacent ground state levels to a common 
upper vibration level, then the intensity ratio for 
two adjacent peaks can be calculated with aid of 
Boltzmann’s law, assuming equal statistical 
weights, to be 0.03. Since the observed ratio is 
so much greater and, furthermore, because of the 
lack of any observable temperature coefficient, it 
is concluded that the constant differences which 
are observed are to be attributed to excited 
vibrational levels which tend to group themselves 
in such a way that the groups are separated by 
constant energy increments. Curry and Herzberg! 
stated that they were able to follow a series of 
bands in the visible, which were separated by 
740 cm. Very probably the band groups which 
have shown up in this work with a separation of 
700+100 cm can be identified with their 
frequency difference. Although the absorption 
maxima are found to be fairly complex, the 
envelope of the bands gives the simple picture 
shown in Fig. 2 in which the maxima are regularly 
spaced. It also appears to be significant that the 
progression of band groups is such as to give a 
minimum of absorption roughly 250 cm~! from 
the maximum and with an intensity 90 percent as 
great, for ail peaks. 

The advice and help of several members of the 
Physics Laboratory are gratefully acknowledged. 
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With a view to obtaining the internuclear distances and 
dissociation energies of individual bonds in polyatomic 
molecules purely from a knowledge of force constants the 
consequences of assuming that the potential energy of any 
di-atom may be expressed in the form V=—a/r™+8/r" 
have been further investigated. In addition to giving a 
relation between internuclear distance and force constant 
of the form suggested empirically by Badger, this assump- 
tion leads to a relation between force constant (.), inter- 
nuclear distance (r-) and dissociation energy (De) of the 
form D.=(kere)/mn. It is found that such a relation does 
exist between those quantities for series of diatomic mole- 
cules belonging to the same row in the periodic table, pro- 


vided one uses an “effective” internuclear distance re—dij. 
This was also used by Badger in his relation of force 
constant to internuclear distance and the values of d;; 
obtained by the two different methods are in reasonably 
good agreement. The values obtained for the product mn 
are 4, 6.2, 6.2 and 10.8 for the respective series LiH to FH, 
Liz to F2, Naz to Cl, and Kz to Bre. For the hydrides it 
would seem that —a/r+8/r‘ is the correct form with 6 
varying along the series. It is also shown that this function 
leads to a relation between the shift in frequency 6v and 
the change in bond energy 8D in the formation of a hydro- 
gen bond. 


T has been clear for some time that unless 
some very striking advance is made in the 
direction of obtaining higher resolving power in 
the infra-red, the future application of vibration 
spectra to problems of molecular structure will 
be largely concerned with the interpretation and 
use of force constants derived from purely vi- 
brational data. Indeed for complex molecules, 
even if the rotational fine structure of the bands 
could be obtained, the resulting information on 
the three principal moments of inertia would be 
of little value since the various internuclear dis- 
tances could not be obtained and it is the inter- 
nuclear distance which is now recognized as the 
most important characteristic of a chemical bond. 
For this reason the various attempts! which have 
recently been made to establish an empirical 
relation between the length of a bond and its 
force constant are of considerable interest and 
importance. The fact that a relation should exist 
at all between these two quantities for series of 
related di-atoms indicates that the potential 
energy curves must have certain common fea- 
tures. It was with this in view that the present 
author proposed? that the potential might be 


* Leverhulme Fellow. 

' Badger, J. Chem. Phys. 2, 128 (1934), 3, 710 (1935); 
Douglas Clark, Phil. Mag. 18, 459 (1934); D. Clark and 
J. L. Stoves, ibid. 22, 1137 (1936), 27, 389 (1939) ; Huggins, 
J. Chem. Phys. 3, 473 (1935), 4, 308 (1936). 

* Sutherland, Proc. Ind. Acad. Sci. 8, 341 (1938). 


represented by a simple function of the form 


a £6 
V=——+—. (1) 
g= 


From such an assumption it is easy to derive that 
the following relations should hold between the 
equilibrium internuclear distance, (r,) the force 
constant (k,) and the dissociation energy (D.): 


ma(n—m) nB(n—m) 


(2) 
kere” 
(3) 
mn 


The purpose of the present paper is to investigate 
a little more fully the consequences of such an 
assumption and the relation of the above equa- 
tions to the observed facts. 

We have essentially two different empirical 
proposals for the connection between internuclear 
distance and force constant, that of Badger viz. 


(4) 
(re—di;)* 


and that of Douglas Clark, and Allen and 
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Longair viz.* 
Ci;' 


k= (5) 


In order to correlate our theoretical relation (2) 


3The constant C;; in both formulae takes different 


values according to the rows in the periodic table to which 
the atoms belong and d;; in Badger’s formula varies in a 
similar way. In the earlier papers quoted it will be found 
that Clark’s relation is given in the form we=C;;’/./ Nr.’. 
Such a form is dynamically unsatisfactory and Dr. Clark 
informs the author that he has recently modified it to be 
essentially as given here. The symbol WN is defined by 
Clark as the number of electrons in molecular rather than 
atomic orbitals. It is omitted by Allen and Longair. 


with one of these findings we must decide that 
either m or n is constant for a series of molecules 
made up from atoms in fixed rows of the periodic 
table. Remembering that » must be greater than 
m the most likely alternatives are m=1 or n=4. 
The former would give agreement with Badger’s 
relation provided we write the potential function 
with an effective rather than a true internuclear 
distance. It would also imply that a and n are 
constant and that only 6 varies as we go from 
one member to another in the series of related 
molecules. The latter would give agreement with 


- Clark’s relation without the use of an effective 


Ye 


1.9 


2.0 3.0 


408 


Fic. 1. Graph showing linear relation between (D./k,)+ and 7, for diatomic molecules in same row 
of the periodic table. Intercept on r, axis gives dj; and slope of the line gives (mn)* in the relation 
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internuclear distance but since N varies the only 
conclusion to be drawn regarding m, a and Bis 
that a and either m or 8 or even both may vary 
in a related series. 

In order to get further information on this 
point we examine what happens for a small 
variation in 7,. According to Badger’s relation we 
should have 


ér. 1 6k, 
while Clark’s would give 
or. 1 dk, 
te 


Now the C—H bond is slightly different in 
methane and in acetylene. In the former‘ 
r-=1.093A and k,=4.9X105 dynes/cm while in 
the latter® we have r,=1.057 and k,=5.9X105 
dynes/cm. This gives 


bk, or. 


or, 
=17%, —=3.4%, 


e le Le ij 


=5%. 


It will be seen that these figures slightly favor 
Badger’s form of the relation. Unfortunately 
there are no other reliable cases in which this 
test might be applied since it is extremely sensi- 
tive to minute variations in either the distance or 
the force constant but as data accumulate it 
should be possible to use it as a discriminant. 
We consider next how the dissociation energy 
varies for a small variation in 7r,. Thus we have 
from 
a 6D. —ma(n—m) 
D,.= —V.=—(1—m/n), 
Ore 


But from (2) 
(m+2)ma(n—m) 


or. remts3 


and so finally 


(6) 


‘For distance, see Ginsburg and Barker, J. Chem. Phys. 
3, 668 (1935); for force constant, see Voge and Rosenthal, 
J. Chem. Phys. 4, 137 (1936). 

For distance, see Herzberg, Patat and Verleger, Zeits. 
f. Physik 102, 1 (1936); for force constant, see Colby, 
Phys. Rev. 47, 388 (1935). 


By substituting for D, and for k, this relation 
may be thrown into the form 


év 
(7) 


y(2—m) 


6D.=B 


2 
+2) ma(n — 2/(m+2) 
[ma(n—m) ] 


In these calculations we have assumed that a is 
constant and £ varies. It is fairly clear that the 
opposite assumption would have led to the 
relations 


6D. n 5k, év 
=——-—— and 6D,=B' 
D. nt+2k, 


(8) 


y(2—n) (2+n) 


If now we substitute the most probable value for 
each case viz. m=1 and n=4 we obtain the 
alternative forms 


év 
6D.=B— and 6D.=B'viév. (9) 


v 


There is as yet too little data to test these 
equations. Badger and Bauer® have proposed the 
empirical relation 


6D = B"— (10) 
v 
between the interaction energy and the fre- 
quency shift 6v when a hydrogen bond is formed. 
But until a relation has been established between 
the change in energy of the strong XH bond and 
the weak XH .. . Y bond has been established, 
we can make no use of Badger’s suggestion. 
We next consider the relation 


kr? 
D.= 


mn 


and whether it is fulfilled for those diatomic 
molecules of which the dissociation energies have 
been accurately determined. Unfortunately there 
are remarkably few diatomic molecules for which 
D, is accurately known and for polyatomic mole- 
cules the dissociation energies of particular bonds 
are still more uncertain. In Fig. 1 we have 


plotted (D./k,)! against the internuclear dis- 


6 Badger and Bauer, J. Chem. Phys. 5, 839 (1937). 


tance r. and we should expect the points for 
molecules in which the two atoms come from the 
same row in the periodic table to lie on a straight 
line the slope of which should be (mn)-}. If the 
relation is true as it stands then the straight lines 
should go through the origin but if one should 
use an effective internuclear distance then the 
intercept on the 7, axis will give that defect which 
has to be subtracted from the actual distance, 
viz. d;; in Badger’s notation. As far as the some- 
what meager data go it would appear that an 
effective rather than a true internuclear distance 
should be used and we have drawn the lines 
through the points to find the d;;’s and compare 
them with those used by Badger. The comparison 
is given in Table I which also records the values 
for the product mn deduced from the slopes of 
the lines. It should be added that the values of 
the dissociation energies of the alkali hydrides 
are not absolutely certain. The data employed 
have all been taken from Herzberg’s recent com- 
pilation.? There is still some doubt about the 
dissociation energy of Pz and so the points corre- 
sponding to the two alternative values have both 
been given. 

It will be noticed that the greatest deviation 
from any of the lines occurs for F2 and this results 
from the anomalously large internuclear distance 
which has been reported for this molecule from 
electron diffraction. Chlorine is next worst but 
if the lower value for the dissociation energy of 
P. turns out to be correct, then the line would 
have to be redrawn and the deviation would be 
considerably reduced while the resulting value 
for di; would also be in better agreement with 
Badger’s. It is interesting that the point for ICI 
falls on the Bre line. The values obtained for mn 
are all of a reasonable magnitude. In the case of 
the hydrides the excellent agreement with 


TABLE I. 
Liz Naz Ke LiH 
d;; (here) 0.57 1.02 1.13 0.34 
d;; (Badger) 0.68 1.18 1.35 0.335 
mn 6.2 6.2 10.8 4 


7 Molecular Spectra (Prentice Hall, 1939). 
8 Brockway, J. Am. Chem. Soc. 60, 1348 (1938). 
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40 20 30 40 a 


Fic. 2. Graph showing linear relation between (D./k.)4 

and r, for similar molecules. Line A refers to HF, HCl, 

HBr, HI; line B to Os, Se, Ses, Tez; and line C to Liv, 

Nae, Ke. A similar line can also be drawn for Cl2, Br, and 
I, but F2 does not lie on it. 


Badger’s d;; would indicate that m=1 and n=4 
gives a satisfactory representation. 

Until more experimental data are available it 
is not possible to be sure how safe it will be to 
extrapolate or interpolate from the above results. 
This method should at least prove useful in 
giving first estimates of dissociation energies. 
Thus it should be noticed that satisfactory 
straight lines may also be drawn connecting re- 
lated elements in the periodic table (Fig. 2) and 
also related hydrides. Although the reason for 
this is not yet clear on the above theory it can 
probably be used as an additional check in pre- 
dicting dissociation energies. The extension of 
this to polyatomic molecules where it will find its 
greatest application will be discussed in a sub- 
sequent paper. 

The author wishes to express his gratitude to 
the Trustees for a Leverhulme Fellowship during 
the tenure of which the above work was carried 
out. He also wishes to acknowledge the benefit 
of discussing these matters with Professors R. M. 
Badger and L. Pauling of the California Institute 
of Technology. 
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Exchange Reaction of Gaseous Bromine and Hydrogen Bromide 


LAURENCE C. LIBERATORE AND Epwin O. WuG 
Department of Chemistry, University of Rochester, Rochester, New York 


(Received August 10, 1939) 


The exchange reaction of bromine and radioactive hydrogen bromide in the gas phase has 
been found to fall just short of reaching equilibrium in two minutes. The very rapid rate of 
exchange would seem to be accounted for best by a chain reaction involving bromine atoms. 


HE exchange reactions of halogens with 
their corresponding halides in aqueous solu- 
tions have been studied by several investiga- 
tors.'~* In all cases exchange has been found to 
take place freely. Topley and Weiss‘ studied the 
exchange reaction of hydrogen bromide and bro- 
mine in carbon tetrachloride and found that 
equilibrium was reached in less than two minutes. 
The present gas phase measurements in the 
absence of water vapor were undertaken in the 
hope that the gaseous exchange might be some- 
what slower which would permit us to make a 
kinetic study of the thermal and photochemical 
bromine-hydrogen bromide exchange. This hope 
was not realized since our experiments indicate 
that equilibrium is almost, but not quite, reached 
in two minutes. Nevertheless, our results may be 
of interest. Brezneva, Roginsky, and Schilinsky® 
allowed radioactive bromine and hydrogen bro- 
mide to remain in contact for fifteen minutes and 
found that equilibrium was established according 
to the masses of bromine in the two gases. 


EXPERIMENTAL 


Radioactive bromine was obtained by bom- 
bardment of selenium with protons in the De- 
partment of Physics cyclotron. The bombarded 
selenium was finely ground and placed in a bulb 
attached to an all glass system using a minimum 
number of stopcocks and comprising a mercury 
diffusion pump (backed by an oil pump), ma- 
nometer, McLeod gage, reaction vessel, trains 


‘Grosse and Agrusse, J. Am. Chem. Soc. 57, 591 (1935). 

* Long and Olson, J. Am. Chem. Soc. 58, 2214 (1936). 
(19 sn Shiflett and Lind, J. Am. Chem. Soc. 58, 535 
as podeon and Fowler, J. Am. Chem. Soc. 61, 1215 

. Topley and Weiss, J. Chem. Soc. 912 (1936). 

Brezneva, Roginsky and Schilinsky, Acta Physico- 
chim., U.R.S.S. 5, 549 (1936). 


for preparing and purifying HBr and Bro, etc. 
Pure, dry HBr gas was admitted to the selenium 
and the mixture heated to bring about exchange. 
In the first six runs this HBr, containing radio- 
bromine, was distilled into the 600-cc spherical 
reaction bulb which already contained the de- 
sired pressure of bromine frozen down in a small 
side tube. After reading the HBr pressure on a 
phosphoric acid manometer the condensed bro- 
mine was quickly vaporized by placing alcohol 
at room temperature around the trap. To stop 
the reaction the two gases were condensed with 
liquid air. The time of reaction was taken from 
the moment bromine began to vaporize to the 
time when both gases were completely condensed. 
Obviously the actual time of reaction was some- 
what less than this. 

Following reaction, the pressure of the HBr 
was measured, keeping the Bre condensed, after 
which the HBr was completely absorbed in a 
tube containing aqueous NaOH. The Bry. pres- 
sure was also measured, followed by freezing it 
down in another tube containing aqueous NaOH. 
This tube was sealed off, the cooling mixture 
removed and the bromine allowed to react. The 
bromine in each solution was completely pre- 
cipitated as AgBr. The precipitate was filtered, 
dried, pressed into a disk about 15 mm in di- 
ameter, and the disk covered with Cellophane. 
The radioactivities of the two disks were meas- 
ured on an ionization chamber and compared. 

In runs 7-13 HBr at high pressure was quickly 
admitted to the bromine vapor in the reaction 
vessel. The gases were allowed to react for two 
minutes, frozen down and handled as before. 

Run 13 was made in the presence of oxygen. 
Since the gases were very difficult to condense 
after reaction, as much of them as possible was 
condensed and the rest completely pumped off. 
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The radioactivities were determined as before. 
In runs 10 to 12 the bulb B was filled with enough 
tubing to increase the surface exposed to the 
gases by a factor of seven. 

In certain runs, as shown in the table, the 
radioactivity was measured on y-rays, because it 
was thought that perhaps all the 6-rays did not 
emerge from the disk. However, negligible differ- 
ences were found in the ratios of the radio- 
activities of those samples which were measured 
both on 6-rays and on y-rays. 


RESULTS 


In the exchange reaction 
HBr*+Br-Br=HBr+Br-Br*, 


the probability that the radioactive Br atom will 
appear in the Br molecule in the forward reac- 
tion is twice the probability that it will appear in 
the HBr molecule in the reverse reaction. There- 
fore at equilibrium the radioactivity is dis- 
tributed between the HBr and Brz in proportion 
to the number of Br atoms present, and if the 
Br2 pressure is one-half that of the HBr, the 
radioactivities of the two resultant AgBr samples 
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are the same. If the gases are not present in the 
ratio 2HBr : 1Bre, then at equilibrium the frac- 
tion of the total radioactivity due to HBr should 
be numerically equal to the fraction of the total 
number of Br atoms present as HBr. 

The calculations in the last two columns of 
Table I were made as follows. The radioactivity 
of the HBr was calculated to what it would have 
been if the pressure of the HBr were twice that 
of the Bre. This value was then used in calculat- 
ing the total radioactivity and the fractions due 
to BrBr* and to HBr*. If reaction has reached 
equilibrium, the result is a value of 0.5 for the 
fraction of the total radioactivity due to either 
Bre or HBr. 

No change in the data could be detected in 
experiments in which the bromine and hydrogen 
bromide were allowed to remain frozen down, 
before analysis, for varying lengths of time up 
to an hour. This would seem to eliminate ex- 
change in the solid state. 

As can be seen from Table I, the thermal ex- 
change is not quite complete in two minutes. The 
pressures stated are those obtained after reaction 
had occurred. In the experiments in which the 
gases were left mixed for 15-20 minutes the 


TABLE I. Summary of experimental data. 


Rapio. oF Br2 
Bras Br: | _Rapio. oF Brs 
No. TYPE MIN. MM MM Br Rapio.* 2HBr : 

a 30 82 34 0.454 0.473 0.510 
1 Photochemical-Unpacked 20 95 48 0.503 0.507 0.505 
5 20 68 42.5 0.556 0.560 0.506 
65 111.5 42 0.429 0.465 0.535 
6 15 95.5 55 0.535 0.530 0.496 
3 2 79 52 0.468 0.490 0.423 
7 eileen 2 47 35 0.598 0.538 0.423 
8 2 66 48 0.593 0.517 0.424 
9 2 93 50 0.518 0.483 0.465 
13 Thermal-Unpacked +O, 2 90 48 0.516 0.450 
10 2 61 52 0.630 0.565 0.432 
11 Thermal-Packed 2 69 38.5 0.528 0.466 
12 2 69.5 38 0.523 0.458 


* Radioactivity measurements in 4, 5, 6, 8, 10, 11, 12, 13 made on gamma-rays. 
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radioactivity in the Bre is just that required for 
equilibrium. In the 65-minute experiment the 
Brz is somewhat more radioactive than expected ; 
this effect is apparently real and not due to ex- 
perimental error, but needs to be checked. The 
reaction appears to be homogeneous as there is 
no effect observed in packing the bulb with Pyrex 
tubing. The experiments indicated as photo- 
chemical were carried out in broad daylight; 
those marked “thermal” in a semi-darkened room 
with the reaction vessel completely covered. 

Since Br*® decomposes to give krypton, the 
radioactive hydrogen bromide would dissociate 
into a hydrogen atom plus krypton. Later con- 
siderations show that the effect of these hydrogen 
atoms is small and, as might be expected, any 
influence of the added oxygen in experiment 13 
would not be apparent. 

The total number of radioactive bromine 
atoms present can be estimated. A number of 
workers’: * have shown that radioactive Br®? 
having a half-life of 4.5 hours emits almost pure 
gamma-rays. These gamma-rays cause internal. 
conversion of electrons in the bromine atom so 
that an electron is removed, forming a radioac- 
tive bromide ion with a half-life of 18 minutes. 
This bromide ion changes to krypton with the 
loss of a beta-particle. Of the above particles only 
the beta-particle causes a current in the ioniza- 
tion chamber, along with the beta-particle 


emitted by radioactive Br® (half-life 32 hours). 


In calculating the total number of radioactive 
atoms present at the time of reaction, the amount 
of each isomer present must be calculated. A 
sample calculation follows. 


At time of measuring, activity Br®°=activity 
Br® (from half-lives) 

Ionization chamber sensitivity =100 per 
sec. per cm deflection (from standard U 
source) 

Sixty-five percent 6’s due to 
Br® get through sample 

Eighty-fi t B’s due to 
ghty-five percent B thick 
Br® get through sample . 

Or, 66.5 6’s per sec. per cm deflection due to 
Br® emitted in sample and 


| Estimated for 


" Segré, Halford and Seaborg, Phys. Rev. 55, 321 (1939). 
* DeVault and Libby, Phys. Rev. 55, 322 (1939). 


66.5 ’s per sec. per cm deflection due to Br® 
emitted in sample. 
Deflection = 263 cm (Total radioactivity) 


dN 1 66.5X3600X4.5 


dt X 
X 263 =4.0 x 108 


at time of measuring. This gives, by the 
decay law, 
5X10* radioactive Br®® atoms at time of 
reaction. 
Similarly, 
reaction. 
Total N(Br®°+Br®) =3.46X10° at time of 
reaction. 


N(Br*®) =2.96X10® at time of 


There are two ion pair producing processes 
occurring in the reaction vessel : (1) ionization by 
B-particles and (2) ionization by the internal con- 
version electron ejected by the y-ray. Because of 
its much smaller energy any effect of the recoil 
ion may be neglected. 

The number of ion pairs due to beta-particles 
was calculated from the ionization chamber de- 
flections as follows (the internal conversion elec- 
trons do not reach the ionization chamber). The 
energy of the beta-particles is rather low (2 Mev), 
so that the number of ion pairs formed per centi- 
meter in air is about 50° at normal temperature 
and pressure. The number of ion pairs formed 
was assumed to be proportional to the atomic 
number of the gas present. Thus the number of 
ion pairs formed in the reaction bulb is 

pressure temp. atomic No. Br 3 


50x x x xX- 
760 273 atomic No. air 4 


No. per cm 


deflection Xdeflection in cm Xradius of bulb 
Xtime of reaction in seconds. 


The three-fourths factor is used because three- 
fourths of the atoms are bromine, the rest being 
hydrogen. 
131 300 35 
N=50X—X— X— X133 X 263 
760 273 8 
X5X120=6.5 X 10° ion pairs 


® Rutherford, Chadwick and Ellis, Radiation from Radio- 
1930) Substances (Cambridge University Press, London, 
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formed by the beta-particles from Br*®* and Br® 
in the vessel in two minutes. . 

The number of ion pairs produced by the in- 
ternal conversion electrons may also be calcu- 
lated. According to measurements just completed 
by Valley'!® these internal conversion electrons 
have an average energy of about 40,000 ev and 
are about 75 percent emitted. This gives an 
average energy of 30,000 ev available and assum- 
ing about 30 ev per ion pair, we have an average 
of 1000 ion pairs per disintegration of Br®° (4.5 
hr.). From the concentration of Br®® and the 
decay law the number of disintegrations in two 
minutes is 2.5X10°, which gives 2.5X10°® ion 
pairs. 

The total number of ion pairs produced in the 
vessel in the two minutes the gases are mixed is 
thus 3.15 10°. 


DIscCUSSION 


The data may be summarized as follows: the 
thermal exchange is unaffected by increasing the 
surface sevenfold; equilibrium is not quite 
reached in two minutes; in an average experi- 
ment of this duration there are present ca. 
3.5 X 10° radioactive Br atoms; the total number 
of ion pairs formed in the vessel during this time 
is 3.15 10°. 

Two mechanisms have been suggested to ac- 
count for radiochemical reactions, viz., the 
cluster theory of Lind" and the primary forma- 
tion of atoms and ions followed by thermal 
reactions (Eyring, Hirschfelder and Taylor"). 

For the sake of simplicity let us assume a 
cluster of only three molecules (one of these being 
a positive or negative ion) as represented by 


HBr*+HBr+Brz 


* 
=( Bra) =BrBr*+2HBr. 
When such a cluster breaks up due to neutraliza- 
tion of the charge (by a cluster of opposite charge 
or by an ion or electron) the probability of the 


10 We are indebted to Dr. George E. Valley of the De- 
Smee of Physics for supplying us with these unpub- 
ished data. 

1S. C. Lind, The Chemical Effects of Alpha Particles and 
Electrons (Chemical Catalog Co., New York, 1928). 

2 Eyring, Hirschfelder and Taylor, J. Chem. Phys. 4 
479, 570 (1936). 
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cluster yielding BrBr* is equal to that of giving 
HBr*. If every ion pair forms such a cluster and 
each contains a radioactive atom, then 1.6X10° 
“tagged” atoms should appear in the Bre and 
actually at equilibrium 1.8 X10° “‘tagged”’ atoms 
do appear. However, the probability of finding 
a Br* atom in the cluster is given by the ratio 
of the number of Br* atoms present to the 
total number of HBr molecules present or 
(3.5 & 10°) /(3.2 X10”) = 110-”. The rate of the 
exchange reaction is thus about 10” times as 
great as such a simple cluster would predict. 
Obviously, no adjustment of the composition of 
the cluster can begin to account for this large 
difference. 

The rapidity of the exchange would seem to 
require a chain mechanism. Chain reactions 
would readily be initiated by the hydrogen and 
bromine atoms formed as a result of the neutral- 
ization of the ion pairs. These effects will be 
approximately the same as discussed in detail for 
a-particles by Eyring, Hirschfelder and Taylor.” 
The primary process in the case of a-particles 
consists of (1) excitation of molecules leading to 
a rupture into atoms and (2) a primary ionization. 
The latter is followed by secondary processes 
which lead finally to atoms. Thus an a@-particle, 
passing through Br and HBr forms by ionization 
processes, according to these authors, about 6 Br 
atoms and 2-3 H plus 3 Br atoms, respectively. 
To these must be added dissociation without 
ionization. Similarly, for 6-particles, we may 
expect from each ion pair several Br or H atoms 
or both. The modified cluster theory of Lind 
and Livingston," “ in which the cluster, after 
recombination of the positive and negative ions 
with the liberation of the energy of ionization, 
is thought of as a small isolated volume of gas 
momentarily at high temperature, might lead to 
a dissociation into atoms. The source of the 
atoms does not affect our consideration of the 
exchange reaction, except as it affects the number 
of atoms produced. 

The molecular species present are HBr, Bre, 
HBr*, and BrBr*, the ratio of the first two to 
the last two being of the order of 10". Conse- 
quently any “tagged” atoms produced as a direct 


13 Livingston and Lind, J. Am. Chem. Soc. 58, 612 


(1936). 
\ Earl F. Ogg, J. Phys. Chem, 43, 399 (1939), 
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result of the various dissociation processes can 
be neglected. Hydrogen atoms can play no role 
except to furnish Br atoms, since the only reac- 
tions which H can undergo are 


H+Br.—-HBr+Br 
H+HBr—H:+Br 
H+BrH—HBr+H, 


which would have no net effect. 
The only reactions in which Br atoms can 
participate are 


Br+BrBr’—BrBr+Br’ 
Br+HBr’—HBr+Br'’, 


the superscripts being used simply to differentiate 
between two atoms. These exchange reactions 
might well be expected to have only small or 
negligible energies of activation. The only effect 
then of a collision of a Br atom with ordinary 
Br2 or HBr is to pass a Br atom along. Collision 
of Br with a molecule containing a “tagged” 
atom may result in exchange as 


Br+HBr*—BrH+Br* 
and Br+BrBr*—BrBr+Br*. 


The reverse reactions would also occur. Gas 
phase recombination of atoms at the pressures 
used would be small, the atoms eventually 
disappearing at the walls. 

Equilibrium may be looked upon as being at- 
tained by setting free all the Br* atoms, after 
which the equilibrium distribution would be 
instantly attained by collision of Br* with HBr 
and Bre. This would require at least 3.5 10*! 
collisions of bromine atoms with HBr in two 
minutes, since the total number of HBr* is 
3.5 10°, and the ratio of HBr/HBr* is 110”. 
This corresponds to 5 X10"* collisions of Br and 
HBr per cc per second and assumes that exchange 
occurs at every collision. 

Aside from negligible thermal effects, there are 
four processes occurring which may result in the 
formation of atoms: (1) decomposition and 
ionization by the y-rays emitted, (2) liberation 
of an atom when radioactive disintegration with 
the ejection of a B-particle occurs, e.g., 


HBr*—H+Kr+8 
or BrBr*—Br+Kr-+48, 


(3) neutralization of the ion pairs produced by 
(a) the internal conversion electrons and (b) the 


and also 


B-particles and (4) excitation of molecules by the 
B-particle. The first two are negligible compared 
to (3) and (4), as is also any effect of the recoil ion. 

We shall assume an average of 6 Br atoms pro- 
duced per ion pair (Eyring, Hirschfelder and 
Taylor estimate 6 Br in Bre and 2-3 H+3 Br in 
HBr produced per ion pair from a-particles). 
Since about equal amounts of the §-particle 
energy go into ionization and excitation, an 
average total of 12 Br atoms should be formed 
per ion pair by both ionization and excitation 
processes. This corresponds to about 810° Br 
atoms formed in the 600-cc reaction vessel in the 
two-minute period by the 6-particles. Similarly, 
assuming 6 Br atoms per ion pair, the 2.510° 
ion pairs produced by the internal conversion 
electrons yield 1.510! Br atoms. All three 
processes thus yield 2.310" Br atoms in two 
minutes or an average concentration over this 
period of 1.1510! Br atoms in 600 cc. Assuming 
this concentration of Br atoms and a collision 
diameter of 3A, the number of collisions of Br 
and HBr per cc per sec. is 1.210!*. This is in 
surprisingly good agreement with the minimum 
number of collisions required, viz., 5X10!° 
ese.*. 

Two assumptions, which are certainly only 
approximations, have been made: that for the 
exchange reactions of Br with HBr and Bre the 
energy of activation is zero and that practically 
no Br atoms are lost by recombination. The 
inclusion of these factors would make the agree- 
ment worse, which may mean that the estimates 
of the number of Br atoms produced by the 
various processes is too small or that the collision 
diameter is too small, or both. The chain sug- 
gested might also involve ions instead of Br 
atoms. This contribution should be small, how- 
ever, since the number of ions produced is smaller 
than the number of Br atoms and the rate of 
recombination of the former much greater than 
that of the latter. 

A possible alternative explanation that sug- 
gests itself is that exchange occurs through the 
reaction 


BrBr+HBr*—BrBr*+HBr. 


However, calculation of the number of such 
collisions at the walls indicates that the hetero- 
geneous reaction could have an energy of activa- 
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tion not greater than 3 kcal. while for the homo- 
geneous reaction the energy of activation must 
be less than 15 kcal. to account for the observed 
rate. Both values seem far too low. The wall 
reaction also seems unlikely because no effect 
of increased surface was observed. 
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We wish to express our thanks to Professor L. 
A. DuBridge and his associates in the Depart- 
ment of Physics for supplying us with radioactive 
bromine, permitting us the use of their ionization 
chambers, and for many stimulating discussions 
and suggestions. 
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Proton-Attracting Properties of Liquids. II 
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(Received November 11, 1939) 


Forty-two different organic liquids are rated according 
to their proton-attracting power by comparing the strength 
of the hydrogen bonds which they form with heavy methyl 
alcohol. The classes of liquids studied are: aldehydes, 
ketones, esters, ethers and amines. The amount of per- 
turbation which the different solvents cause in the OD 
vibrational band of the CH;OD is taken as a measure of 
the strength of the deuterium bridges formed with the 
CH;OD. The study is a continuation of a previous work in 
which twenty-six similar liquids were studied, and the results 
agree with those of the former study in showing that ac- 
cording to proton-attracting power the different classes are 


INTRODUCTION 


HE present paper is a continuation of the 

spectroscopic comparison of the power of 
different proton acceptor solvents to form hydro- 
gen bonds with a common proton donor solute.! 
These studies should assist in the prediction of 
many important properties of solvents because 
the power of the solvents to form hydrogen 
bonds with a proton donor group is largely 
conditioned by its tendency to donate an 
electron pair. According to the modern concep- 
tion of a base as a proton acceptor, or electron 
pair donor,? the ability of a solvent to form 
hydrogen bonds with a proton donor would 
depend upon its strength as a base. Hence this 


* Department of Physics, Mary Hardin-Baylor College, 
Belton, Texas. 

+ Department of Chemistry, The College of Wooster, 
Wooster, Ohio. 

1W. Gordy, J. Chem. Phys. 7, 93 (1939). 

2K. W. Sherk, J. Chem. Educ. 13, 358 (1936); J. N. 
Bronsted, Chem, Rev. 5, 23 (1928), 


in the order: esters<aldehydes and ketones<ethers 
<amines. From the study it is possible to observe the 
effects of unsaturation, substitution, and other structural 
variations on the electron donor power of the O or N atom 
to which the OD group forms a bond. It is shown that the 
strength of the deuterium bonds formed with CH;OD 
bears a close relation to the basicities of the solvents, to 
their solubilities for other proton donor solutes, and to 
their reactivity with other compounds. There appears to 
be no correlation between the dipole moments of the 
solvents and the strengths of the bonds formed with the 
CH;OD. 


application of infra-red spectroscopy offers a new 
and direct method for comparing the intrinsic 
strength of bases. The recent comprehensive 
studies of Zellhoefer, Copley and Marvel’ on 
solubility of haloforms in electron donor solvents 
enable us to compare solvent power with the 
shift of the OD band. The remarkable agreement 
between the amount of perturbation brought 
about in the OD band and the solubility in the 
solvent suggest that abnormal solubilities depend 
to a marked extent on the electron donor power 
or the ability to form a complex with the 
solute. 

In the previous paper it was shown that 
according to proton-attracting power the different 
classes of liquids are in the order: nitro com- 
pounds <esters <aldehydes and ketones < ethers 
<amines. In the present work a more extensive 
study of the different classes is made. In general 


3M. J. Copley, G. F. Zellhoefer and C. S. Marvel, J. Am. 
Chem. Soc. 60, 1337, 2666, 2714 (1938). 
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the results bear out the conclusion previously 
reached, although the wider survey reveals many 
interesting variations depending upon the struc- 
ture of the individual members of each class. 


EXPERIMENTAL PROCEDURE 


The resolving instrument is the same as that 
described in the previous paper.! As in the earlier 
work on CH,OD solutions a 1 molar solution was 
found to be most satisfactory. For concentrations 
much more dilute than this the thickness of the 
cell required is such that the general absorption 
and harmonic bands of the solvent decrease the 
accuracy of the results. On the other hand, for 
a more concentrated solution a smaller per- 
centage of the alcohol molecules would be 
linked to the solvent. This, too, would tend to 
make the results uncertain. It was found that 
more reliable results could be had for most of 
the solutions by using a slightly thicker cell 
than was used before. Consequently a thickness 
of 0.008 cm was used instead of 0.006 cm. The 
curves given in the figures represent the ratio of 
transmission of the solution to that of an 
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373.9 41 43 3739 41 
WAVE-LENGTH 


3739 44 435 3.73.9 41 43 


Fic, 1. Effects of solvent on the OD band of CH;OD. 
Broken line, 0.1 M (t=0.32 cm) CH;OD in benzene. 
Solid line, 1 M (t=0.008 cm) CH;OD in: A. m-tolualde- 
hyde; B. cinnamic aldehyde; C. furfural; D. citral; 
E. heptaldehyde; F. 2 ethyl butraldehyde; G. methyl 
iso-butyl ketone; H. di-iso-butyl ketone; I. methyl amyl 
ketone; J. mesityl oxide; K. L. p-methyl 
acetophenone; M. cyclohexanone; N. fenchone. 
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PERCENT TRANSMISSION 


3.73941 45 3739 41 43° 3739 45 43 373.9 41 
WAVE -LENGTH 


Fic. 2. Effects of solvent on the OD band of CH;,OD. 
Broken line, 0.1 M (t=0.32 cm) CH;OD in benzene. 
Solid line, 1 M (t=0.008 cm) CH;OD in: A, iso-amyl 
formate ; B. ethyl chloroacetate ; C. y-bromopropy] acetate ; 
D. phenyl acetate ; E. furfuryl acetate ; F. methyl benzoate; 
G. benzyl benzoate; H. n-butyl pthalate; I. diethyl 
oxalate; J. ethyl carbonate; K. chute phosphate; L. 
ethyl n-butyl ether; M. 88’ dichloroethyl ether ; N. dibenzyl 
ether; O. phenetole; P. epichlorohydrin. 
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Fic, 3. Effects of solvent on the OD band of CH,OD. 
Broken line, 0.1 M (t=0.32 cm) CH,OD in benzene. 
Solid line, 1 M (¢=0.008 cm) CH;OD in: A. diamyl amine; 
B. triamyl amine; C. benzylamine; D. dibenzylamine; 
E, aniline; F. o-chloroaniline; G. m-chloroaniline; H. di- 
methyl aniline; I. o0-toluidine; J. ethyl o-toluidine; K. 
methylanthranilate; L. pyrrole. 


approximately equal thickness of the pure 
solvent. 

The bands are fairly sharp and the minima 
easily determined for all the compounds except 
the aliphatic amines and quinoline. Here the 
experimental difficulties are increased by over- 
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TABLE I. CH;OD solutions, 


SuirT FROM 
MOoNOoMOLECU- 
Postrion | LAR BAND IN 
or OD 
Banp SoLuTIoN 
SoLvENT Ap Av 
*Benzene CeHe 2681 
*Pure liqnid CH;OD 4.01 | 2494} 0.28 187 
*Nitrobenzene CeHsNO2 77 | 2653) 0.04 28 
Aldehydes 
*Benzaldehyde CeHs-CH : O 3.85 | 2597| 0.12 84 
m-tolualdehyde CH; 3.85 | 2597) 0.12 84 
Cinnamie aldehyde CcHsCH : CH-CH : 0 3.86 | 2591) 0.13 90 
Furfural CsH30-C 3.83 | 2611) 0.10 70 
Citral (CHC: CH: ‘CH: *C(CH3) : CH-CHO 3.88 | 2577| 0.15 | 104 
Heptaldehyde CH3(CH2)s -CHO 3.90 | 2564) 0.17 | 117 
2 ethyl butraldeh: = -CHO 3.88 | 2577] 0.15 | 104 
*n-butraldehyde C2Hs-CH2-CH 3.90 | 2564) 0.17 | 117 
Ketones 
*Methyl e ketone *CO-C2Hs 3.84 | 2604/ 0.11 77 
*Acetone 3.87 | 2584) 0.14 97 
Methyl 3.84 | 2604) 0.11 77 
Di-iso-butyl ketone 3.85 | 2597) 0.12 84 
Methyl amy! ketone C;Hn -CO-CHs; 3.84 | 2604) 0.11 77 
Mesityl oxide (CHs)2C : CH-CO-CHs 3.87 | 2584) 0.14 97 
CcHs -CO-CHs 3.84 | 2604) 0.11 77 
p-methyl acetophenone 3.84 | 2604) 0.11 77 
Cyclohexanone CH2(CH2)4-CO 3.90 | 2564) 0.17 | 117 
Fenchone CioH1cO 3.85 | 2597) 0.12 84 
Esters 
*Ethy! acetate CHsCO2-C2Hs 3.85 | 2597) 0.12 
*Amyl acetate boy 3.86 | 2591} 0.13 90 
Iso-amy] formate HCO2-C:H: 3.83 | 2611) 0.10 70 
Ethyl chloroacetate Cl-CH2- “CO:CoHs 3.83 | 2611) 0.10 70 
bromo-propyl acetate 3.85 | 2597) 0.12 84 
Phenyl acetate CHsCOsCcHs 3.84 | 2604 O11 | 7 
Furfuryl acetate CHsCO2-C.H30 3.82 | 2618} 0.09 63 
Methyl! benzoate CeHsCO2-CHs 3.82 | 2618] 0.09 63 
Benzyl! benzoate CsHsCO2-CH2-CeHs 3.82 | 2618! 0.09 63 
n-butyl pthalate CeH«(CO2CsHo)2 3.82 | 2618] 0.09 63 
Diethyl oxalate (COsC2Hs)2 | 3.82 | 2618} 0.09 63 
Ethyl carbonate CO(OC2Hs)2 3.80 | 2632} 0.07 49 
n-butyl phosphate PO(OC«Ho)s 3.90 | 2564) 0.17 | 117 
Ethers 
*Diethyl ether (C2Hs)20 3.92 | 2551] 0.19 | 130 
Ethyl n-butyl ether CoHs-O-CiHo 3.90 | 2564] 0.17 | 117 
*Iso-propy! ether ory -CH:0 3.91 | 2558] 0.18 | 123 
(C4Hs)20 3.89 | 2571} 0.16 | 110 
66’ dichloroethy! ether (Cl-CH2-CH2)20 3.85 | 2597] 0.12 84 
Dibenzy] ether (CsHsCH2)20 3.90 | 2564) 0.17 | 117 
Phenetole C2H;0-CcHs 3.81 | 2624] 0.08 57 
*Dioxan CH2OCH2CH20CH2 3.87 | 2584) 0.14 
Epichlorohydrin O-CH2-CH -CH?Cl 3.88 | 2577] 0.15 | 104 
Amines 
Diamy] amine yo me 4.02 | 2488} 0.29 | 193 
Triamy! amine (CsHi1)3N 4.00 | 2500} 0.27 | 181 
Benzyl amine CeHsCH2-NH2 4.07 | 2457] 0.34 | 224 
Dibenzylamine (CsH;CHz)2NA 4.05 | 2469) 0.32 | 212 
*Pyridine CsHsN 4.00 | 2500) 0.27 181 
*a-picoline CH3-CsHiN 4.02 | 2488} 0.29 | 193 
Piperidine CH2(CH2)sNH 4.10 | 2439) 0.37 | 242 
Aniline CeHs-N 4.00 | 2500} 0.27 181 
o-chloroaniline “CoH 3.96 | 2525) 0.23 | 154 
m-chloroaniline Cl -CeHs-NH2 3.97 | 2519} 0.24 | 162 
Dimethyl aniline (CH3)2N -CeHs 3.94 | 2538) 0.21 | 143 
o-toluidine CH3-CeHs-NH2 3.99 | 2506) 0.26 | 175 
Ethyl o-toluidine -NH-C:H; 3.95 | 2532] 0.22 | 149 
Methylanthranilate NH2-CeH«-CO2CHs 3.97 | 2519] 0.24 | 162 
Pyrrole (CH : CH)sNH 3.96 | 2525) 0.23 | 156 


* Repeated from previous paper, J. Chem. Phys. 7, 93 (1939). 
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lapping of the intense CH absorption of the 
solvent and by the broadness and shallowness of 
the band. 


DISCUSSION OF RESULTS 


The transmission curves for aldehydes and 
ketones are given in Fig. 1; for esters and ethers, 
in Fig. 2; for amines, in Fig. 3. With each of 
these curves is given a curve for CH;OD in an 
inert solvent, benzene (0.1 molar solution; cell 
thickness, 0.32 cm). The positions of the minima 
in these bands, with the shift from the position 
of the monomolecular band in the benzene 
solution are listed in Table I. Although there is 
considerable variation in the perturbation of the 
OD absorption made by the different members 


_ of each class, the classes themselves show notice- 


able unlikeness to each other. For example, the 
shifts produced by the esters are consistently 
lower than those produced by the aldehydes and 
ketones, which are about the same. For the most 
part, the ethers cause greater shift than the 
aldehydes and ketones; and the amines greater 
than the ethers. 


Effects of unsaturation 


The cases studied do not permit a direct 
comparison to ascertain the effect of an aliphatic 
double bond. Citral produces a slightly smaller 
shift in the OD band than does heptaldehyde, the 
saturated aldehyde most nearly like it of those 
studied. Cinnamic aldehyde has two carbon 
atoms between the carbonyl group and the 
benzene ring. For this reason it may be expected 
to behave more nearly like an aliphatic than an 
aromatic aldehyde. The fact that it produces a 
smaller shift than do the aliphatic aldehydes may 
be a result of its unsaturation. As these cases are 
not conclusive we expect to make further study 
of the aliphatic double bond. 

A comparison of the shift produced by pyridine 
with that produced by the similar but fully 
hydrogenated compound, piperidine, shows that 
the effect of unsaturation here is to decrease 
appreciably the proton-attracting power of the 
nitrogen. This result is in agreement with the 
lower basicity constant of pyridine as compared 
with that of piperidine. 

In general the effect of the attachment of a 
ring containing double bonds is to decrease the 
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perturbation produced on the OD band. As 
would be expected the effects are less pro- 
nounced the further the unsaturated ring is 
removed from the atom involved in hydrogen 
bond formation. Furfural, for example, causes 
less shift than any of the other aldehydes 
studied, while benzaldehyde and m-tolualdehyde 
produce slightly less shift than the aliphatic 
aldehydes. 


Effects of substitution 


The effect of the substitution of a chlorine for 
a hydrogen may be seen by comparing diethyl 
ether with dichloroethy] ether: 


H H H 


| | 


| | 
HH HH 


shift 0.19 


HH HH 


| | | | 


| J 
HH HH 


shift 


also by comparing ethyl acetate (shift 0.12) 
with ethyl chloroacetate (shift 0.10), and aniline 
(shift 0.274) with o-chloroaniline (shift 0.23) or 
m-chloroaniline (shift 0.244). The influence of 
the strongly electronegative chlorine is trans- 
mitted through the carbons to the oxygen or 
nitrogen and results in a decrease of the effective 
negative charge on the electron donor atom. 
In this connection it should be mentioned that 
the chlorine atoms in chloroform are believed to 
make the CH group sufficiently acidic to share 
its proton in hydrogen bonding.‘ 

The substitution of a phenyl group for an 
ethyl group causes a very marked decrease in 
the electron donor power of the ether oxygen. 
This may be seen by comparing diethyl ether 
with phenetole: 


C,H s—O—C:;H by 5 


shift 0.194 shift 0.084 


*S. Glasstone, Trans. Faraday Soc. 23, 200 (1936); W. 
Gordy, J. Chem, Phys. 7, 163 (1939). 


On the other hand, the influence of the phenyl 
group in decreasing the electron donor power of 
the oxygen is practically nullified when a carbon 
is placed between the phenyl group and the 
electron donor atom. Compare methyl ether with 
acetophenone : 


O O 
| LUA 
CH;—C—C.H,;, CH 


shift 0.114 shift 0.11y 


There is very little difference in the shifts caused 
by ethyl acetate and phenyl acetate: 


O O 
CH;—C—O—CG3H;, CH 


shift 0.124 shift 0.114 


a fact which probably indicates that the bonding 
is to the carbonyl oxygen of the acetate, which is 
somewhat removed from the phenyl group. 

Similarly, a benzyl when substituted for an 
ethyl group has only a slight effect as compared 
with a phenyl group. Compare diethyl ether 
(shift 0.194) with dibenzyl ether (shift 0.17,). 
That a phenyl group decreases the effective 
negative charge of the electron donor atom much 
more than does a benzyl group may be seen by a 
comparison of aniline (shift 0.274) with benzyl 
amine (shift 0.34). This result is in complete 
accord with other physical data. For example, 
phenol is much more acidic—less basic—than 
benzyl alcohol. 

To compare a furyl group with a methyl 
group, notice ethyl acetate and furfuryl acetate: 


CH 3—C—O—C;H 
shift 0.124 
O O 
| 
CH 
H—C 
shift 0.094 


to compare it with a phenyl group, notice furfural 
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(shift 0.104) and benzaldehyde (shift 0.12,). 
The effect of the furyl ring was mentioned above 
in the discussion of unsaturation. 

The shift produced by benzaldehyde is the 
same as that produced by m-tolualdehyde; and 
that produced by acetophenone as that by 
p-methyl acetophenone. These cases suggest that 
the substitution of a methyl group for a hydrogen 
has little if any effect. In these cases, however, 
the substitution is remote from the atom to 
which the OD is bonded. Naturally the most 
pronounced effect takes place when a hydrogen 
joined directly to the electron donor atom is 
replaced by a methyl group. Compare aniline 
(shift 0.274) with dimethyl aniline (shift 0.21,). 
From the latter example it would seem that the 
substitution of a methyl group for a hydrogen 
tends to decrease the electron donor power of the 
aniline nitrogen. Yet a-picoline, in which a 
methyl group is attached to a carbon adjacent 
to the nitrogen in the pyridine ring, causes a 
slightly greater shift (0.294) in the OD band 
than does pyridine (shift 0.274). The large 
difference in the effects of aniline and dimethyl 
aniline may be due to unlikenesses in the steric 
hindrance of the methyl group and the hydrogen 
rather than to unlikenesses in their effect on the 
electron donor power of the nitrogen. 

As a class the esters, with two oxygens linked 
to the same carbon, cause less effects upon the 
OD band than do the ketones, with one oxygen 
linked to the carbon. Diethyl carbonate, with its 
three oxygens linked to one carbon, produces a 
smaller shift than do the other esters. Phos- 


TABLE II. Comparison of dipole moment of solvent with 
shift in OD band. 


MOMENT SHIFT IN BAND 


SOLVENT X<10'8 e.s.u. 


Bromobenzene 
Nitrobenzene 
Benzaldehyde 
Acetone 
Acetophenone 
Ethyl acetate 
Amy] acetate 
Diethyl ether 
Dioxan 
Pyridine 
Aniline 
Dimethyl aniline 
Pyrrole 

Methyl cyanide 
Benzyl cyanide 


Wo 


DD 


00 in to 
— 


W. GORDY AND S. C. 


STANFORD 


phorous when joined to the electron donor atom 
appreciably increases its power to attract a 
proton. n-butyl phosphate causes a much greater 
shift (0.174) than any of the other esters. We 
may compare, for example, ethyl carbonate and 
n-butyl phosphate : 


CyH»—O 


C.H;—O 
\ 
Cc = O, C,H,—O—P =O. 
C,Hy—O 
shift 0.07 shift 0.174 


The pronounced variation here is probably not 
entirely due to the differences in the groups 
attached to the carbon and phosphorus but at 
least partly to differences in the electro-negativity 
of these two atoms. From other data’ carbon is 
known to be more electro-negative than phos- 
phorus. 


Other structural effects 


From a comparison of the shift caused by 
dioxan with that caused by ethyl, propyl and 
butyl ethers, it would seem that a cyclic ether is 
not as effective an electron donor solvent as is an 
open chain ether. This is borne out, but to a less 
degree, by epichlorohydrin. Except for cyclo- 
hexanone, 


there is not much variation among the different 
ketones. It produces a considerably greater shift 
than do the other ketones. Piperidine, the 
amine which is structurally similar to this com- 
pound, correspondingly produces a large shift of 
the OD band as compared with the other amines. 
Fenchone, with the carbon of its carbonyl group 
forming part of a ring structure, causes about the 
same shift as the ordinary ketones. 

The nitrogen in the pyrrole ring has an electron 
donor power equivalent to that of the least 


5L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, 1939), p. 60. 
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effective of the aromatic amines. Compare, for 
instance, the shift produced by pyrrole (0.23,) 
with that produced by o-chloroaniline (0.23). 


COMPARISON WITH OTHER PHYSICAL PROPERTIES 


Dipole moment of solvent 


In Table II we have arranged a comparison of 
the dipole moments of several typical solvents 
with the shifts of the OD band which they pro- 
duce. There seems to be little if any relation 
between the two. Nitrobenzene, which has the 
largest dipole moment of any of the solvents 
studied, produces the smallest shift of any 
except bromobenzene, while dioxan, with the 
smallest dipole moment of any of those listed, 
causes a shift of 0.144. Both pyridine (dipole 
moment 2.2) and aniline (dipole moment 1.55) 
produce a shift of 0.274. Considered as a class 
ethers have much lower dipole moments than 
ketones and aldehydes, yet they produce greater 
shifts of the OD band. Similarly, the amines, 
with lower dipole moments than the ketones, in 
general cause greater shifts than any of the 
other compounds studied. The nitriles and nitro 
compounds, with the highest dipole moments, 
produce relatively small shifts. This comparison 
indicates that the complex formation here is not 
due simply to dipole attraction, and that the 
tendency of a solvent to acquire protons from a 
hydroxyl group bears little relation to the dipole 
moment of the solvent. 


Basicity constants 


Few of the aromatic amines absorb appreciably 
in the region of 4y, so that it is possible to make 
reasonably accurate measurements on the OD 
band for CH;,OD solutions in these solvents. 
Also for most of them the band is sufficiently 
sharp and symmetrical for its center to be located 
without difficulty. In Fig. 4 we have plotted 
basicity constants against shifts of the OD band 
for ten aromatic amines and one aliphatic amine. 
The other amines for which basicity constants 
are available in the literature have bands so 
complex and shallow, or so affected by intense 
absorption of the solvent, that determination of 
the centers is uncertain. Because the basicity 
constants used are for water solutions they 
represent essentially the measure of the ability 


103 
1o* 


BASICITY CONSTANTS 


= 


20 25 30 35 
WAVELENGTH SHIFT INJJ 


Fic. 4. Relation of basicity constant of solvent to shift 
in OD band. © dimethyl aniline; © o-chloroaniline; 
m-chloroaniline ; 0-toluidine; A aniline; A pyridine; 
O a-picoline; © tributylamine; a benzylamine; @ 
piperidine. 


of a base to capture a proton from water. 
Except in certain cases where difference in steric 
hindrance may cause slight variations—for ex- 
ample, repulsion of the methyl groups of dimethyl 
aniline for the CHs;OD may explain why the 
point for this solvent lies off the curve—we see no 
reason why the tendency of different solvents to 
capture protons from water, as indicated by their 
basicity constants, should not vary in the same 
way as their tendency to capture protons from an 
hydroxyl group of alcohol, as in either case it 
would be conditioned by the electron donor 
power of the solvent. The shift of the OD band 
made by these different solvents is the measure of 
the weakening of the binding force of the 
hydroxyl group and is, therefore a measure of 
the tendency of the base to capture the hydroxyl 
proton. Assuming the linear relationships sug- 
gested in Fig. 4, we may obtain by extrapolation 
basicity constants for the oxygenated solvents 
recorded in Table I. This method gives the 
following values for Ky: nitrobenzene, 4 10-*; 
benzaldehyde, 8 X 10-*°; acetone, 2 aceto- 
phenone, 1X10-*°; ethyl acetate, 8X10-*°; iso- 
propyl ether, 9X10-"*. 


Solubility 


In their studies of the solubility of monofluoro- 
dichloromethane (CHCI.F) in electron donor 
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TABLE III. Relation of solubility to shift in OD band. 


SOLUBILITY OF 
CHChF in OD BAND 
SOLVENT MOLE/FR. Ap 


Nitrobenzene 0.360 0.04 
Benzaldehyde 0.436 0.12 
Furfural 0.367 0.10 
Heptaldehyde 0.519 0.17 
Cyclohexanone 0.548 0.17 
Furfuryl acetate 0.438 0.09 
Diethyl oxalate 0.556 0.09 
68’ dichloroethy] ether 0.451 0.12 
Dioxan 0.479 0.14 
Phenetole 0.425 0.08 
Aniline 0.258 0.27 
Dimethy] aniline 0.425 0.21 


solvents, Zellhoefer, Copley and Marvel,** in- 
vestigated many of the same solvents which are 
used in this study. They interpreted the solu- 


. bilities which they found to be in excess of 


Raoult’s law as being due largely to hydrogen 
bonding between the CHCI.F and the solvents. 
In Table III we have arranged a comparison of 
their data on solubilities of CHCI,F in 12 
different solvents with our data on the shifts 
which these same solvents produce on the OD 
band of CH;OD. For most of these solvents the 
correlation is remarkable, and is strong evidence 
that the interpretation of Zellhoefer, Copley and 
Marvel is correct. For diethyl oxalate the 
solubility obtained (0.556) is larger than that 
which the small shift of 0.094 would seem to 
indicate. In this molecule, however, there are two 
carbonyl groups with the possibility of two 
molecules of CHCI2F tying to one molecule of the 
solvent. This possibility could, we believe, ex- 
plain the lack of correlation in this case. The 
solubility of CHCI:F in aniline (mole fr. 0.258) 
is less than the ideal solubility predicted by 
Raoult’s law (mole fr. 0.381), whereas this 
solvent brings about a rather large shift (0.27) 
in the OD band. Here the NH—N hydrogen 
bonds between the aniline molecules® are proba- 
bly stronger than the CH—N bonds which would 
be formed with the CHCI.F. On the contrary, 
OD-N bonds formed between CH;OD and 
aniline would probably be stronger than the 
hydrogen bridges linking the aniline molecules. 
Thus the marked difference in data on solubility 
and shift is additional evidence that the proton 
donor power of the groups studied is in the order: 


6 W. Gordy, J. Chem. Phys. 7, 167 (1939). 


CH <NH <OH.’ As pointed out in the discussion 
of the relation of shifts to basicity constants, 
CH;0D does not form as strong a bond with 
dimethyl! aniline as one would predict from its 
basicity constant. It was suggested that this may 
be due to differences in steric hindrance of this 
solvent for water and methyl alcohol. The 
solubility of CHCl.F in dimethyl aniline, lower 
than would be predicted from the shift which the 
solvent produces on the OD band of CH,OD, 
may be caused by a difference in steric hindrance 
of the dimethyl aniline for the CHCl.F and the 
CH,OD. 

Although they did not use the identical 
phosphate, they found, in agreement with our 
data, that a phosphate forms much stronger 
bonds than other esters. Other points of agree- 
ment are: that the chlorines in 8’ dichloroethy! 
ether decrease the electron donor power of the 
ether oxygen; that a cyclic ether has a weaker 
electron donor power than an ordinary ether; 
that a phenyl group attached to an oxygen or 
nitrogen atom markedly decreases its electron 
donor power ; that the effect of unsaturation is to 
decrease the electron donor power slightly ; and 
that in their power to donate electrons the classes 
of oxygenated solvents studied are in the order: 
nitro compounds < esters < aldehydes and ketones 
<ethers.* There seems to be one major point of 
disagreement: these observers state that the 
basicity of the amine solutions in water ap- 
parently is not related to the donor properties of 
the nitrogen in forming hydrogen bridges, 
whereas in our measurements there is a close 
correlation with the donor properties of the 
nitrogen in forming hydrogen bridges. 


Chemical reactivity 


Conant and Bartlett® have compared the rates 
of semi-carbazone formation for some aldehydes 
and ketones. The general reaction which they 


studied may be illustrated by the reaction with 


acetone: 


<(CHs3)2C = NNCONH2+H:0. 


™N. V. Sidgwick, The Covalent Link in Chemistry 
(Cornell University Press, 1933), p. 166. 

8 See also, M. J. Copley and C. E. Holley, Jr., J. Am. 
Chem. Soc. 61, 1599 (1939). 

*Conant and Bartlett, J. Am. Chem, Soc. 54, 2881 
(1932); Conant, Ind. Eng. Chem. 24, 466 (1932). 
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Only four of the eight carbonyl compounds 
studied by these observers were included in the 
present study. Fortunately they are ones which 
show rather wide variations. The results of 
Conant and Bartlett are compared with ours in 
Table IV. In Fig. 5 we have plotted OD shift 
against log ke. Although there are too few cases 
for generalization, the comparison suggests an 
interesting relation. Also, in agreement with our 
data, these observers found furfural to be least 
reactive of all the aldehydes studied, and 
cyclohexanone to be most reactive of all the 
ketones studied. 

Though steric effects undoubtedly cause con- 
siderable variation in certain cases, the most 
important factor in determining the relative 
strengths of the deuterium bonds which these 
different solvents form with CH;OD is probably 
the effective negative charge on the electron 
donor atom, or its readiness to donate its 
electrons. Thus the close relation existing be- 
tween amount of shift of the OD band produced 
by the various solvents and their basicities, 
solubilities and reactivities with certain other 
compounds is not surprising. Indeed all organic 
molecules may be broadly classified as electron 
acceptors or electron donors—i.e., all their 
chemical properties are in some way affected by 
their ability to accept or donate electrons. The 
TaBLE IV. Comparison of reaction rate of semi-carbazone 


formation of aldehydes and ketones with 
shift in OD band. 


VELocITY CONSTANTS OF 
SEMI-CARBAZONE 
FORMATIONS$H yDROLYSIS SHIFT IN 
OD Banp 
ke ki X105 Ap 
Furfural 0.73 =| 0.55 0.10 
Benzaldehyde 2.05 0.62 0.12 
Acetone 6.02 1800.00 0.14 
Cyclohexanone 36.00 7600.00 0.17 


Or 


1 
10 A2 14 16 p 
OD SHIFT 
Fic. 5. Plot of OD shift against log ke of semi-carbazone 


formation: @ furfural; A benzaldehyde; 0 acetone; 
O cyclohexanone. 


method employed here in estimating this im- 
portant property is more direct and far less 
laborious than many other types of measurement 
used for this purpose. 
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Adsorption Isotherms for Mobile Monolayers 
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(Received October 3, 1939) 


A three-dimensional gas is pictured in contact with a liquid, forming at the interface a 
mobile monolayer. The density of the monolayer is maintained in equilibrium with the pressure 
of the gas. An equation is deduced giving as a function of temperature, the gas pressure which 
is just sufficient to initiate condensation in the monolayer. The only physical parameters in- 
volved in the relationship are the mass of the molecule and the heat of evaporation per molecule 
from the completely clustered monolayer. The fatty acids which exhibit monolayers with 
negligible external gas pressures, must have heats of evaporation greater than about 100 X 10-4 
erg. Molecules as small as 20 atomic mass units and heats of evaporation less than about 
75 X 10-4 erg should require pressures of the order of 10-? atmosphere. An adsorption isotherm 
equation is found for monolayer concentrations less than the saturation value. Its form is 
similar but not identical with Langmuir’s isotherm for immobile monolayers. 


INTRODUCTION 


OWLER! has succeeded in deducing Lang- 
muir’s isotherm equation for monolayers in 
which the molecules are adsorbed at definite 
locations on a solid surface. Using the results 


recently found? for mobile monolayers on liquid 


surfaces, it is now possible to apply Fowler’s 
general arguments to this case. The results are: 
first, an equation giving as a function of tem- 
perature, the minimum pressure of gas required 
to maintain the monolayer in its condensed 
state; and secondly, the pressure of the gas as a 
function of the corresponding concentration of 
the unsaturated monolayer, at any temperature 
—in other words the adsorption isotherm for the 
mobile monolayer. 

There do not yet appear to be any data re- 
quiring detailed interpretation by such formulae; 
actual mobile monolayers being produced with- 
out the aid of any three-dimensional gas in 
contact with the liquid. This means merely that 
the equilibrium gas pressure is unmeasurably 
small with existing technique; it does not mean 
that in future mobile monolayers may not be 
produced which require an appreciable pressure 


of gas for their maintenance. 


1Fowler, Statistical Mechanics (Cambridge University, 
1936), second edition, §21.53. 

? Band, “Statistical Mechanics of Condensation Phe- 
rere in Mobile Monolayers,” J. Chem. Phys. 8, 116 
1940). 


GENERAL THEORY 


Let the complete assembly contain X mole- 
cules and possess total energy E. Consider an 
example of the assembly in which M of the 
molecules are attached in the monolayer, and 
N=X-—M are in the gas. The analysis must 
find the example, specified by M, for which the 
number of complexions is a maximum.* 

Write f.(z) for the partition function of a 
single cluster of type s in the monolayer, and 
g(z) for the partition function of a molecule in 
the gas. Write x; for the variable corresponding 
to the selection of a state in the monolayer, and 
x2 for the selector variable corresponding to the 
gaseous state. 

Then the number of complexions correspond- 
ing to a particular assignment of energy E, in 
the gas and E,, in the monolayer will be, in the 
example considered, the product of the coefficient 
of in exp {> .x1‘f.(z)} and the coefficient 
of x2”z¥g in exp {x2g(z)}. To include all possible 
assignments of energy, will require the coefficient 
of 28, E=E,,+£,, in the product of the two 
expressions. This leads to the integral expression 
for the total number: 


f f dx 
x 


Xexp {xog(z)}-exp (1) 


3 Fowler, reference 1, §5.4, uses a similar argument in 
the study of equilibrium between a crystal and its vapor. 
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ADSORPTION 


Using, without rigorous proof, the theorem of 
Darwin and Fowler, reference 1, §2.42; it then 
follows that: 


In Cy = — M In 
(2) 


where Ai, Ae, ¢ are the unique roots of the equa- 
tions: 


In Cy = 
—E/¢=0, (3) 


In Cy =g(¢) —N/d2=0. (5) 


The example with largest Cy will be that in 
which N is adjusted to give 


In Cy =0. (6) 


Since M=X—VN, and X is not varied, this 
means that 
In Ae=0, 


or Ai =A2=A say. (7) 
Eqs. (3)—(5) then become: 
(8) 
(9) 
N=hg(f). (10) 


From these equations the law of “mass action”’ 
follows by elimination of \ between (9) and (10): 


From this equation either concentration, mono- 
layer or gas, can be found in terms of the other 
under equilibrium conditions. The usual argu- 
ment, from (8), gives the familiar relation be- 
tween the parameter ¢ and the absolute tem- 
perature: 

¢=exp (—1/kT). 


MintmuM GAs PRESSURE FOR CONDENSED 
MONOLAYER 


Making use of the arguments employed before,* 
the partition function for the monolayer cluster 
type s is: 


f(T) =s*A(1—B,)/p exp (.—¢1—sW)/kT. (12) 


‘Reference 2, Eqs. (1), (2), and (8)-(11). 
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For the meaning of p see reference 2, Eq. (19). 
The energy of a cluster of type s is sW—4,, the 
term ¢; occurs in (12) to cancel a similar term 
in p. Classically 


g(T) =(2amkT)iVh-, (13) 


where V is the volume of the gas. The internal 
vibratory states of the molecules have been 
omitted from (12) and (13) because, if they are 
assumed unaffected by adsorption, they will 
cancel from Eq. (11). This last equation now 
becomes : 


Xexp (¢.—¢1— W)/kT 
X {(N/V)h exp (14) 


The critical value of this expansion occurs when 


(N/V) 
=(2rmkT)'h- exp (W/kT). (15) 


The critical sum corresponding to substitution of 
(15) into (14) gives the result: 


(M/A). exp ($.—1)/RT. (16) 


This equation gives information concerning the 
monolayer condensation, and should agree with 
the result of the previous work.? To prove this, 
note that (16) gives 


(M/A) <p/(1—B1) = —B,)/(1— 
Xexp (¢:—91)/kT (17) 


exp (¢.— 1) /RT 
= 1+0*, (18) 


in the notation of reference 2, Eq. (23). But (18) 
is identical with Eqs. (20) and (22) of the paper 
quoted, showing that the results are the same. 

Equation (15) above leads to the desired 
relation between the critical gas pressure and 
the temperature at which the monolayer is to 
be maintained in the saturated state. Assuming 
that the gas pressure is so low that it obeys the 
perfect gas law, it will follow at once from (15) 
that 


exp (W/kT). (19) 


Again from reference 2, Eq. (9), it will be seen 
that the quantity W means the negative heat of 
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evaporation from the completely clustered mono- 
layer; it is remarkable that this is the only 
physical parameter, except for the mass m, upon 
which the minimum gas pressure depends as a 
function of 7. The equation is of course valid 
only for temperatures below the critical point for 
the phase change in the monolayer. 


THE ADSORPTION ISOTHERM 


It is usual to define the concentration of an 
immobile monolayer as the ratio between the 
actual number of sites occupied by adatoms, and 
the total number of sites. Here the relative 
concentration of a mobile monolayer at any 
given temperature will be defined as the ratio 
of the actual number of adsorbed molecules to 
the number required to produce saturation at 
that temperature. Note that this is not the ratio 
to the number in the condensed state; this 
part of the theory does not cover conditions in 
the monolayer after condensation has set in. 

Write for the relative condensation 


0=(M/A)/(M/A)., (20) 
and use (18) to deduce that 
(21) 
Then (14) gives 


0= {()*/(1+0*)} 
Xexp (¢.—¢1)/RT, (22) 


where v=(N/V)/(N/V).. (23) 


This series converges for ¥ <1. It can always be 
made to converge as rapidly as desired by taking 
y small enough; which means physically taking 
conditions far enough from saturation. Assume 
conditions such that it is enough to take only 
the first two terms of the series, and write 
approximately : 


exp (¢2—¢1)/RT. (24) 


Since far enough from saturation the 8, can all be 
neglected, these terms have been omitted from 
(24), and o written for o*, see reference 2, 


Eq. (27). 


In actual cases (¢2—¢1)/kT is about negative 


3, so Eq. (24) can be written approximately 


6(1+0)=y/(1—ay) (25) 
or (26) 
where a=8 exp (¢2—¢1)/kT. (27) 


To express these in physical terms, note the 
definition (23) and the value (15) of (N/V)., then 
writing N/V=P/kT, (26) becomes: 


(28) 
where A=(2rm)*h-*(kT)5 exp (W/kT) (29) 


has the same value as the constant in Langmuir’s 
isotherm, namely 


P=A0/(1—-6). (30) 


Eq. (28) is in the form of an adsorption isotherm 
also, but differs slightly from Langmuir’s, (30), 
not only in mathematical form, but physically 
because of the different meaning of @. 


NUMERICAL PREDICTIONS 


In practice the fatty acids form mobile mono- 
layers under conditions where gas pressures are 
negligibly small. From (19) this means that the 
(negative) value of W must be relatively large. 
A simple calculation shows that for molecules 
with mass greater than about 200 units, values 
of —W greater than about 100X10-" erg will 
suffice to make the minimum pressure P less 
than 10-‘ atmosphere, when TJ is around 300°K. 
For pressures to reach as high as 10-? atmosphere 
at room temperatures presumably smaller mole- 
cules will be required, say around 20 mass units, 
and in this case —W will have to be smaller 
than about 75 X10-" erg. 

Even in this case a condensed monolayer will 
be possible only if the critical point temperature 
for the phase change is above the range con- 
sidered, and this can hardly be determined in 
advance. 

The dependence of P as a function of W is 
thus quite sharp. Should suitable substances be 
found and the required technique be developed 
for producing monolayers and dealing with the 
equilibrium here discussed, Eq. (19) should pro- 
vide a very useful method of analysis. 
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The Entropy of Long Chain Compounds in the Gaseous State* 


Maurice L. HuGoins 
Kodak Research Laboratories, Rochester, New York 


(Received October 28, 1939) 


Equations are derived for the entropy of long chain compounds in the gaseous state at 
25°C. (1) The translational entropy, computed classically, varies with chain length according 
to the relation Ss=a+6 In n. (2) The entropy of rotation of the molecule as a whole, calculated 
statistically, assuming either a ball-like molecule or absolutely random kinking, is also linear 
with respect to In n. (3) The entropy of internal vibrations, assumed (following Pitzer) to be an 
additive function of contributions associated with the various types of bonds and bond angles 
present, obeys the equation S;,=a+bn. (4) The entropy of internal randomness, for randomly- 
kinked molecules with no rotation, with hindered rotation, or with completely free rotation, 
is likewise linear in n. The equation deduced for Jong chain paraffins, containing but one adjust - 
able constant (connected with the degree of hindrance of the rotation about the C—C bonds) 
agrees quite well with values deduced from specific heat data for short chain paraffins. 


INTRODUCTION 


ANY of the properties of high polymers 
and other substances consisting of very 
long molecular chains depend to a large extent 
on the entropy. A knowledge of the dependence 
of the entropy on the state, composition and 
chain length is, therefore, of prime importance. 
As a step in this direction, this paper deals with 
the dependence on chain length and composition 
of the entropy in the gaseous state of long chain 
compounds which may be assumed to be kinked 
in a random or nearly random manner. 

The total entropy may be considered to be 
equal to the sum of (1) the translational entropy, 
S1, (2) the entropy of internal vibrations, S,, 
(3) the entropy of external rotation—that is, of 
rotation of the molecule as a whole—S,,, and 
(4) the entropy of internal randomness, S;,, due 
to alternative orientations and more-or-less free 


rotation and oscillation about the interatomic 
bonds. 


THE TRANSLATIONAL ENTROPY 


The molal translational entropy of a dilute 
gas is given’: * by the equation 


*Contribution No. 692 from the Kodak Research 
Laboratories. Much of the material in this article was 
Loseees in a paper before the Division of Physical and 
norganic Chemistry at the Milwaukee Meeting of the 
American Chemical Society, September 7, 1938. 

1G, N. Lewis and J. E. Mayer, Proc. Nat. Acad. Sci. 15, 
208 (1929). 


*W. F. Giauque, J. Am. Chem. Soc. 52, 4808 (1930). 


V \ 
S,=RIn ) 
N 


in which R is the gas constant per mole, in 
calories per degree, V is the molal volume, 
m is the mass of one molecule, T is the absolute 
temperature, k is the Boltzmann constant, h is 
Planck’s constant, and e is the base of natural 
logarithms. 
Into this equation one substitutes for the 
volume 
R'T (2) 
P ’ 


where R’ is the gas constant per mole in cc— 
atmos. per degree and P is the pressure in 
atmospheres, and for the mass 


M. 0” (3) 
m=mon = 


where mo is the mags of one of the  “sub- 
molecules” in each molecule and Mo is the sum 
of the atomic weights of the atoms in the sub- 
molecule. (For a long paraffin chain, (CHe),, Mo 
is 14.03.) 

One thus obtains the equation 


[3/2 


S.=Rin| 


At 298.1°K and 1 atmosphere, this reduces to 
S:=26.0+3R In Mo+8R In n. (5) 
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For a long paraffin chain, 
S:=33.8+3R In n. (6) 


THE VIBRATIONAL ENTROPY 


If the vibrational frequencies are known, the 
vibrational entropy can be calculated in the 
usual way?’ * as the sum of Einstein functions. 
In the absence of the direct experimental fre- 
quencies, one may assume a uniform frequency 
for each C—H bond, another for each C—C 
bond, another for each H—C—H bond angle, 
etc. For a long paraffin chain, tentatively 
adopting Pitzer’s* frequency assignments, the 
entropy contributions given in the last column 
of Table I are calculated. Unfortunately, the 
lowest frequencies, which are most important for 
this purpose, are least well known. The calculated 
S, may therefore be considerably in error. 

For long chain molecules in general, one may 
put 

=k,’ +k n, (7) 


computing rough values of k, and k,’ as in Table 
I, if the requisite fundamental data are ob- 
tainable. 


THE ENTROPY OF ROTATION OF THE MOLECULE 
AS A WHOLE 


The entropy associated with rotation of a mole 
of rigid molecules, all alike, is given by the 
following equation :5~7 


Ser=R In | — 


TO 


TABLE I, Rough calculation of the vibrational entropy for 
normal paraffin hydrocarbons. 


1 
( (8) 


TYPE NUMBER FREQUENCY Sy (298°K) 
CH 2n+2 3000 0.000 

HCH n+4 1440 .064+0.016n 
CC n—1 1000 — 092+ .092n 
HCC 3n—2 950 — .226+ .339n 
ccc n—2 320 —2.628+1.314n 


Total S,= —2.88+1.76n 


*L. S. Kassel, Chem. Rev..18, 277 (1936). 

‘K. S. Pitzer, J. Chem. Phys. 5, 473 (1937). 

°J. E. Mayer, S. Brunauer and M. G. Mayer, J. Am. 
Chem. Soc. 55, 37 (1933). 

6 ys Halford, J. Chem. Phys. 2, 694 (1934). 
93s) L. Eidinoff and J. G. Aston, J. Chem. Phys. 3, 379 
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Here o is the symmetry number, which is unity 
for kinked long chain molecules, and Jx, Jy and 
Iz are the moments of inertia of the molecule 
about mutually perpendicular X, Y and Z axes 
through the center of gravity. 

Randomly-kinked molecules will, in general, 
not have the same moments of inertia. Hence, 
one must average, as follows: 


1/8n*keT\3] R 
-Sepe=R In ) (IxlIvyIz)) (9) 


h? 


Tv 
R 
=284.5-+—(In (IxIyIz))m at 298°K. (10) 


Although the individual moments, Jx, etc., of 
the molecules in a given sample of gas doubtless 
vary considerably from the average moments, 
Ix, etc., the moment products, Ix-Iy-Iz, would 
be expected to be much more closely bunched 
around the average moment product. Like- 
wise, the logarithms of the moment products, 
In (Ix-Iy-Iz), would be expected to be closely 
bunched around the average, (In (Ix:Iy-Iz))w. 
The following approximation, therefore, intro- 
duces but little error. 


(In (Ix:Iy-Iz))w. (11) 


Another justification for the use of this 
approximation is the fact that, on the average, 
randomly-kinked molecules are not very elon- 
gated. Kuhn® has shown that the root-mean- 
square length of the straight line between the 
first and the last atoms is only 1/6 times the 
root-mean-square distance of the middle atom 
from this straight line. (See Eqs. (16) and (19) 
below.) 

Moreover, it seems likely, as Langmuir® has 
pointed out, that the van der Waals attractions 
between different parts of the molecular chain 
cause the molecule, in the gaseous state, to coil 
up into a roughly spherical ball. If so, the three 
moments of inertia of a given molecule are 
nearly equal to each other and to the moments 
of all. the other molecules. Eq. (11) is then a 
very good approximation. 


8 W. Kuhn, Kolloid Zeits. 68, 2 (1934). 

9]. Langmuir, in Colloid Chemistry, edited by Jerome 
Alexander (Chemical Catalogue Co., New York, 1926), 
Vol. I, p. 525. 
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ENTROPY OF LONG 


Two calculations of S,, will be made: one for 
spherical molecules having a uniform density 
equal to the density the substances would have 
in the liquid state at that temperature ; the other 
for symmetrical ‘‘average’’ molecules, all alike, 
having dimensions equal to the root-mean- 
square dimensions calculated on the assumption 
of absolutely random kinking. The actual value 
of S., would be expected to be between these two 
results. 

For spherical molecules of uniform density p 
and radius r, 


Ix =Ty=12=3mr (12) 
and 
(13) 
(4/3) 


Eliminating r from these equations, using Eq. 
(3), and substituting into Eq. (11), leads to the 


relation 
M 
) (14) 


For a long paraffin chain, giving p the value 
0.78 g/cc, 


(In —265.7+5 In 


(In (IxyIz))w= —5.62-+1n ( 


(15) 


For large, randomly-kinked chain molecules, 
the mean square of the distance between the 
first atom and the ith atom is given, except for 
very small values of 7, by the relation®: 1% ™ 


(Ri*)w = Bali, (16) 


where / is the distance (1.54-10-® cm for a 
paraffin chain) between neighboring atoms in the 
chain and B, is related to the angle a between 
each bond in the chain and the extension of 
either of the two adjacent bonds by the re- 
lation 


1+ cos a 
(17) 
1—cos a 
For ‘‘tetrahedral”’ angles, 
B,=2. (18) 


The mean square of the distance between the 
ith atom and the straight line joining the end 
10H, Eyring, Phys. Rev. 39, 746 (1932). 


"E, Guth and H. Mark, Monatsh. f. Chem. 65, 94 
(1934), 
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atoms of a randomly-kinked chain is given 
(see Appendix) by the equation 


except when 7 is small. 

The calculation of moments of inertia is now 
made for a symmetrical molecule, having R? and 
r values equal to the averages in Eqs. (16) and 
(19), relative to mutually perpendicular X’, Y’, 
and Z axes through one end atom, the Z axis 
passing also through the other end atom. These 
moments of inertia are given by the equations 


(19) 


Tx =mo(Ly?+ L327), (20) 
Ty (21) 
and 
Iz=m(0x?+ Ly). (22) 
Now, 
1 
23) 
n 
and 
n 
For large i, 
n 
Li=— (25) 
1 2 
and 
n® 
Le=—. (26) 
1 3 
Substituting, one obtains 
4 
Ix: =Ty, (27) 
and 
I (28) 


The center of gravity of this hypothetical 
symmetrical molecule is on the Z axis, midway 
between the two end atoms. The moments of 
inertia about X and Y axes through this point 
are obtained by subtraction from Ix, and Iy: of 
the product of the mass of the molecule by the 
square of the distance between the X and X’ 
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(or Y and Y’) axes: 


1 7 
Ix = ly = Ix: (29) 
4 36 


It is reasonable to assume that moments of 
inertia calculated according to Eqs. (28) and (29) 
are fairly close approximations to the average 
moments. for randomly-kinked molecules. Sub- 
stituting in Eq. (11) gives 


(In UxlIyIz))w= —5.47+3 In (——). (30) 


For a long paraffin chain, 
(In (IxIyIz))w= — 267.6+6 In n. (31) 


In general, the center of gravity of a randomly- 
kinked molecule is, of course, not exactly at the 
midpoint of the straight line joining the end 
atoms. The average moments of inertia about 
axes through the actual centers of gravity are, 
therefore, different from those just deduced. To 
correct for this, a statistical calculation has been 
made of the average squares of the coordinates 
of the center of gravity, relative to the origin 
at the midpoint of the line joining the end atoms, 
with the following results: 


= yo? ~ 0.034B,.7? (32) 
20° (33) 


Since the decrease in (In (IxIyIz))s due to this 
correction is very small (less than 0.1), the de- 
tails of this calculation are not reproduced. 
Possibly more important is the neglect in the 
above treatment of any allowance for the finite 
sizes of the atoms. It was assumed that all 
orientations (consistent with a constant bond 
angle) of each bond relative to an adjacent bond 
are possible. Actually, some are blocked, be- 
cause of the large repulsions between atoms when 
their centers approach each other too closely. 
As a result, the molecules are, on the average, 
somewhat more extended—less closely packed— 
than if truly randomly kinked.” For a first 
approximation, one may assume a slightly larger 
“effective” angle between adjacent bonds (a 
smaller ‘effective’ a) and hence a _ larger 


#2 L, Laskowski and R. E. Burk, J. Chem. Phys. 7, 465 
(1939). 
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“effective” B, and a larger coefficient of In m in 
Eq. (31). 

This correction may be important for mole- 
cules in the liquid state, but in the gaseous state 
it is probably more than balanced by the correc- 
tion due to the van der Waals attractions within 
the molecule. As a comparison of Eqs. (15) and 
(31) shows, these attractions tend to decrease 
the coefficient of In n. 

The difference between Eq. (15) and Eq. (31) 
is relatively small and unimportant for the 
present purpose. The former is probably more 
nearly correct. 

Substituting Eqs. (14) and (30) into Eq. (10), 
one now obtains at 298°K 


5 5 
on the assumption of spherical molecules, and 
3 
Sa= In Inn, (35) 


assuming random-kinking. Using Eqs. (15) and 
(31), for large paraffin chain molecules, these 
equations reduce to 


5 
Ser=21 +S Inn (36) 


and 
Ser=234+3R In n. (37) 


ENTROPY OF INTERNAL RANDOMNESS 


One must still take account of the entropy due 
to alternative orientations, rotations and oscilla- 
tions around the bonds within the molecules. 
This entropy contribution would be expected to 
be a linear function of n: 


Sir= +Risn. (38) 


The constants must depend on the flexibility of 
the chain and so on the size and shape of the 
energy humps resisting free rotation. 

For the hypothetical case of rigid molecules, 
with three equally probable orientations around 
each bond, 


Sie=R In = —642.18n. (39) 


For the opposite extreme—the hypothetical 
case of absolutely free rotation about each bond 
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—a theoretical calculation has been made, 
yielding 


3 
In (B..Mol”) 


R 
+[364+— In (40) 


For normal paraffin chains at 298°K, (if there 
were free rotation) 


Sir= —26+9n. (41) 


Since the method of derivation is quite complex 
and involves certain somewhat questionable 
approximations, it will not be described here. 
The writer believes, however, that the results 
obtained are not far from correct. 

For actual chains, in which the rotations are 
known to be considerably restricted; one would 
expect constants intermediate between those of 
Eq. (39) and those of Eq. (40)—or Eq. (41) for 
paraffins. 

According to Eqs. (39) and (41), the simple 
relation 

ki,’ = —3kir (42) 


holds approximately for hypothetical paraffin- 
like chains, in which either free rotation or no 
rotation is assumed. Making the reasonable 
assumption that it also holds for actual paraffins, 
in which the rotation is hindered, one obtains 


—3kirtkin. (43) 


This should be a sufficiently good approximation, 
for very long chains at least. (Although, from 
theoretical considerations, deviations from Eqs. 
(39), (41) and (43) of the order of magnitude of, 
say, 15 units would not be surprising; such 
deviations would be of importance only for 
short chains.) 


THE ToTAL ENTROPY 


Adding together Eqs. (5), (7), (34) and (43), 
the total entropy at 298°K and 1 atmos. of a 
ball-like molecule is 


Me 
S=33+k,/ +R In (—) 
p 


+4R In n+(kotkir)n. (44) 
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For a randomly-kinked molecule, using Eq. (35) 
instead of Eq. (34): 


In 'l) 
9 
(45) 


These sums reduce, for long chain paraffins, to 


S=52—3ki-+8 In n+(1.8+k;,)n (46) 
and 
S=54—3ki- +9 In n+(1.8+k;,)n, 


respectively. 

Although these equations have been derived 
for chain molecules in which m is large, and 
would be expected to be inaccurate for short 
chains, a comparison with the values which have 
previously been deduced for the simplest normal 
paraffins is of interest. Such a comparison is 
shown in Fig. 1, in which the solid line repre- 
sents Eq. (46), with k;, given the arbitrary value 
6.2. The values for the individual points, listed 
in Table II, were obtained largely from specific 
heat data. 

The values given for the last four compounds 
were computed from the S29s° values given by 
Parks,'® using estimates of the heats of vaporiza- 
tion kindly furnished by Rossini,” and vapor 
pressure data from Landolt-Bérnstein’s Tabellen, 
interpolated and extrapolated to 25°C. 

For comparison, curves for rigid chains and 
for absolutely free rotation, calculated on the 
assumption of spherical molecules rather than 


(47) 


TABLE II. Entropies of gaseous normal paraffin hydro- 
carbons at 25°C, 


a K. Witt and J. D. Kemp, J. Am. Chem. Soc. 59, 273 
1937). 
( 4 J. D. Kemp and C. J. Egan, J. Am. Chem. Soc. 60, 
1521 (1938). 

 G. S. Parks, C. H. Shomate, W. D. Kennedy and B. L. 
Crawford, Jr., J. Chem. Phys. 5, 359 (1937). This value 
“‘ probably errs on the side of being a trifle high.”’ 

16G, S. Parks, Chem. Rev. 18, 325 (1936). 

17F, D. Rossini, private communication. 
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for completely random kinking, are also shown in 
Fig. 1. 

From the agreement shown in the figure one 
would seem to be justified in using equations of 
the form of Eq. (46), with appropriate changes of 
the constants, to deduce approximate entropy 
values for very long chain compounds in the 
gaseous state. 


APPENDIX 


Derivation of equations for the average distance 
_ and average square of the distance between 
the ith atom in a randomly-kinked chain 
and the straight line joining the end atoms 


Kuhn has derived the equations for the middle 
atom of the chain. His procedure, suitably 
generalized, is followed here. 

For any value of the bond angle other than 
90° or 180°, the directions of any two bonds are 
independent of each other, provided they are 
separated by a sufficient number of intervening 
bonds. Let the whole chain be divided into N 
segments, each large enough so that the straight 
lines connecting the end atoms of each segment 


MAURICE L. 


2Qi+1. QN-i+1 


HUGGINS 


are independent of each other as regards direc- 
tion. Let the average length of the projection of 
one of these N straight lines on one of these 
mutually perpendicular axes X, Y and Z be 
called b. On the average, each segment adds or 
subtracts a length b to the projection of the chain 
on each of these axes. 

Taking the origin of coordinates at one end 
of the chain, let the Z axis pass through the 
other end also. Consider the x coordinates of the 
atom at the end of the jth segment. Let P, be 
the number of positive values of b and P_ the 
number of negative values of b between the origin 
and this atom. Then 


P,+P_=j. (48) 
Defining a number p by the relation 

P,—P_=2p, (49) 
one deduces 

P,=j/2+p (350) 
and 

P_=j/2—p. (51) 


For the part of the chain between the atom 
being considered and the end (j= WN), similarly: 


P,'—-P_'= —2p, (52) 


P,’+P'=N-j, (53) 
P,'=N/2—j/2—p, (54) 
P_'=N/2—j/2+p. (55) 


The probability of having P, additions and 
P_ subtractions in the first part of the chain and 
P.,' additions and P_’ subtractions in the second 


part is 


(56) 


Substituting Eqs. (50), (51), (54) and (55) and 
using Stirling’s approximation for the factorials, 
one obtains 


(57) 


— 1 +) ( ( 
j j N-j N-j 


ENTROPY OF LONG CHAIN COMPOUNDS 


Taking logarithms, expanding In (1+2p/}), etc., 
dropping terms which are negligible for large j or 
large (N—j), and taking the antilogarithm of the 
result, the following equation is obtained: 


exp (- ). (58) 
j(N-j) 


From the relation 


f Wr,p,'dp=1, 


one deduces 


Substituting this and also 


p=x/2b 
and 
dp=dx/2b 


into Eq. (58), there results: 


1 2N 

26\xj(N—J) 

x ( (63) 
ex xX. 


Similarly, 


1 2N 


2b\1rj(N—j) 


“Ny? 
exp (-— 


dy. (64 


For a given pair of x and y values, 


N(x? +?) 


65 


dxdy. 


The distance r of the atom being considered 
from the Z axis is given by the equation 


(66) 


Hence, 


Nr? 
— }rdr. 


xp ( 


W,dr 


Averaging, 


N 


One now introduces the proportionality 

j/N=i/n (70) 
and the relation 

PN=B.I'n/3, (71) 
obtained by Kuhn* from a comparison of his 
statistical equation for R® with the equation 
derived for this quantity by Eyring.” Substitu- 
tion into Eqs. (68) and (69) gives 


3 n 


This last equation is identical with Eq. (19) 
above. 


(72) 


(73) 
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The problem of directed valence is treated from a group theory point of view. A method 
is developed by which the possibility of formation of covalent bonds in any spatial arrange- 
ment from a given electron configuration can be tested. The same method also determines 
the possibilities of double and triple bond formation. Previous results in the field of directed 
valence are extended to cover all possible configurations from two to eight s, p, or d electrons, 
and the possibilities of double bond formation in each case. A number of examples are discussed. 


INTRODUCTION 


ROBLEMS of directed valence, like most 

problems of molecular structure, can be 
attacked by either of two methods: the method 
of localized electron pairs (Heitler-London) or 
the method of molecular orbitals (Hund-Mul- 
liken). While it is now realized! that these 
methods are but different starting approxima- 
tions to the same final solution, each has its 
advantages in obtaining qualitative results. 
Theories of directed valence based on the 
methods of localized pairs have been developed 
by Slater? and Pauling? and extended by Hult- 
gren.4 The method of molecular orbitals has been 
developed principally by Hund® and Mulliken.® 
These methods have been compared extensively 
by Van Vleck and Sherman.' 

In the previous papers no attempt has been 
made to discover all the possible stable electron 
groups which lead to directed valence bonds, nor 
have the possibilities of double bond formation 
been completely explored. In the present paper 
both of these deficiencies in the theory have been 
removed. 

METHOD 


Pauling’s method consists of finding linear 
combinations of s, p, and d orbitals which differ 
from each other only in direction. Thus he has 


* Presented at the Boston meeting of the American 
Chemical Society, September 15, 1939. Publication assisted 
by the Ernest Kempton Adams Fund for Physical Re- 
search of Columbia University. 

1 Van Vleck and Sherman, Rev. Mod. Phys. 7, 167 (1935). 

2 Slater, Phys. Rev. 37, 841 (1931). 

3 Pauling, J. Am. Chem. Soc. 53, 1367, 3225 (1931). 

4 a. Phys. Rev. 40, 891 (1932). 

5 Hund, Zeits. f. Physik 73, 1 (1931); 73, 565 (1931); 
74, 429 (1932). - 

' 6Mulliken, Phys. Rev. 40, 55 (1932); 41, 49 (1932); 
41, 751 (1932); 43, 279 (1933). 


found combinations which are directed toward 
the corners of a tetrahedron, others directed 
toward the corners of a square, and others di- 
rected toward the corners of an octohedron. 
Electrons occupying these new orbitals can then 
resonate with unpaired electrons occupying or- 
bitals of other atoms lying in the directions of 
these new orbitals and so form covalent bonds 
with these atoms. The further these orbitals 
project in the direction of the surrounding atoms, 
the stronger should be the resulting bonds. 

In order to construct such sets of orbitals, it 
is most convenient to make use of group theory. 
Each set of equivalent directed valence orbitals 
has a characteristic symmetry group. If the 
operations of this group are performed on the 
orbitals, a representation, which is usually re- 
ducible, is generated. By means of the character 
table of the group’ this representation, which we 
shall call the o representation, can be reduced to 
its component irreducible representations. The 
s, p, and d orbitals of the atom also form repre- 
sentations of the group, and can also be divided 
into sets which form irreducible representations.* 

Let us refer to the set of equivalent valence 
orbitals as the set S. If the transformation which 


reduces this set is 7, then the set 7S can be 


broken up into subsets, each of which forms a 
basis for one of the irreducible representations of 
the symmetry group of ©. 


TS=LaSi, (1) 


7 For the group theoretical methods used here see Wig- 
ner, Gruppentheorie (Vieweg, Braunschweig, 1931); Weyl, 
Theory of Groups and Quantum Mechanics (tr. Robertson) 
(Methuen, London, 1931); Van der Waerden, Gruppen- 
theoretische Methoden in der Quantenmechanik (Springer, 
Berlin, 1932). 

8 Bethe, Ann. d. Physik (5) 3, 133 (1929). 
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where ©; belongs to the ith irreducible repre- 
sentation. If we let R be the set of available 
orbitals of the atom, these may always be 
chosen so that they fall into subsets each of 
which is a basis for one of the irreducible 
representations : 


(2) 


If each of the coefficients b; in (2) is equal to or 
greater than the corresponding a; in (1), we may 
form the new subset §’, given by 


R’ =T aN, (3) 


in which each orbital transforms in exactly the 
same way as the corresponding orbital in (1). 
If we now apply the inverse transformation T—! 
to this set we obtain a set of orbitals 7—'R’ which 
must have exactly the symmetry properties of 
the desired set S. 

To illustrate this process, consider the set 
G=i1+02+3 consisting of three valence or- 
bitals lying in a plane and making equal angles 
of 120° with each other. If the atom in question 
has available s, p, and d orbitals, the set ® 
consists of the nine orbitals s, pz, py, pz, dz, des, 
dyz, dzy, dzryy. The symmetry group of © is 
and the set © may be reduced by the trans- 
formation 


=1//3(o1+02+¢3) 
=1/4/6(201—02—<3) (3) 
o3 

Of these orbitals, o,’ belongs to the representa- 


tion® A,’ while and belong to E’. Hence 
we find for (1) 


S=A/'+E£’. (4) 


TABLE I. Character table for trigonal orbitals. 


Dsh E oh 2C3 2S3 3C2 3en 
A,’ 1 1 1 1 1 1 
A,! 1 1 1 1 -1 -1 
A,” 1 -1 1 -1 1 -1 
A," 1 -1 1 -1 -1 1 
E’ 2 2 -1 -1 0 0 
Dg 2 —2 -1 1 0 0 
s 1 1 1 1 1 1 
p 3 1 0 —2 -1 1 
d 5 1 -1 1 1 1 
o 3 3 0 0 1 1 
© 6 0 0 0 —2 0 


* The notation used here is that of Mulliken, Phys. Rev. 
43, 279 (1933). 


The set ® is already reduced, for s and d, 
belong to A,’, p, to Ae”, the pairs pz, p, and 
—d,z4, to E’, and the pair d,., dy. to E”. 
Hence (2) becomes 


R=2A +A +E’ +E". (5) 


Since the coefficients in (5) are each not less than 
the corresponding coefficient in (4), the desired 
directed orbitals are possible. In fact, since we 
have a choice between s and d, for the orbital 


TABLE II. Reduction table for trigonal orbitals. 


Ar Ai” A? Ad’ E’ E” 
5 1 0 0 0 0 0 
p 0 0 0 1 1 0 
d 1 0 0 0 1 1 
o 1 0 0 0 1 0 
x 0 0 1 1 1 1 


belonging to A,’, we may construct two different 
sets of directed orbitals. If we choose for SW’ 
the set 


R’ =stp.tp, (6) 
and apply the inverse transformation 
1 2 
(7) 
J/3 Yo v2 
v3 V6 V2 
we obtain Pauling’s trigonal orbitals 
1 2 
v3 
1 1 1 
ae (8) 
1 1 1 


On the other hand, this theory shows that we can 
everywhere replace s by d,, or in fact by any 
linear combination as+ $d, provided a?+6?=1. 
Similarly p, and p, can be replaced by d., 


and 
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It is not necessary, however, to carry out the 
actual determination of T and the reduction of 
S and § to decide whether or not a given set of 
directed orbitals is obtainable from given atomic 
orbitals, for the coefficients a; and }b; can be 
obtained very simply from the character tables 
of S and Rt. The transformation matrices of the 
o representation can be written down and the 
trace of each gives the character for the corre- 
sponding element of the group. The characters 
for R are easily found by the method of Bethe.*® 
We can thus construct a character table for the 
representations based on © and Rt. We also can 
enter in the same table the characters for the 
irreducible representations, as given for example 
by Mulliken.® Such a character table for the 
plane trigonal case we have been considering is 
shown in Table I. By means of the ortho- 
gonality theorems for group characters, the 
coefficients a; and 0b; can easily be found. These 
are entered in a second table, to which I shall 
refer as the reduction table. The reduction table 
for our trigonal case is shown in Table II. 
Comparing the coefficients a; (given in the row 
marked o) with the coefficients b; (given in the 
rows marked s, p, and d) we see immediately 
that the trigonal orbitals require one orbital 
which may be either s or d, and two p orbitals 
(remembering that the representation E’ is of 
degree two) or two d orbitals. Hence the possible 
valence configurations are sd*, and d’. 

It is interesting to note that the method of 
molecular orbitals leads to identical results, but 
by a rather different route. In this method we 
consider first the set of orbitals on the atoms 
surrounding the central atom. If this set con- 
sists of orbitals symmetrical about the line 
joining each external atom to the central atom, 
then these external orbitals form a basis for a 
representation of the symmetry group which is 
identical with the o representation. The reduction 
of this representation then corresponds to the 
resonance of these external orbitals among them- 
selves. The formation of molecular orbitals then 
takes place by the interaction between these 
reduced external orbitals and the orbitals of the 
central atom. This interaction can only take 
place, however, between orbitals belonging to 
the same representation. Hence, to obtain a set 
of molecular orbitals equal in number to the 


number of external atoms, it is necessary that 
each of the reduced external orbitals be matched 
with an orbital from the set which belongs to 
the same representation. The condition for this 
is again that b;2a;, so that the same result is 
reached as before. 

The possibilities of double or triple bond 
formation are most easily discussed in terms of 
molecular orbitals. The principal type of multiple 
bond consists of two parts: first, a pair of 
electrons in an orbital symmetrical about the 
axis of the bond; and second, one or more pairs 
of electrons in orbitals which are not sym- 
metrical about the axis. The first pair of electrons 
form a bond which differs in no way from the 
ordinary single, or ¢, bond. The other pairs are 
ordinarily in orbitals which are antisymmetric 
with respect to a plane passed through the axis. 
They may be regarded as formed by the inter- 
action of two p orbitals, one on each atom, with 
axes parallel to each other and perpendicular to 
the axis of the bond. We shall refer to orbitals 
of this type as 7 orbitals. 

In a polyatomic molecule, consisting of a 
central atom and a number of external atoms 
bound to it, bonds of this type may also be 
formed. As far as the external atoms are con- 
cerned, the condition for the formation of bonds 
is the presence of » orbitals at right angles to 
the bond axes. These p orbitals, however, will 
resonate among themselves to form new orbitals 
which are bases of irreducible representations of 
the symmetry group of the molecule. This 
reduction can be carried out in the same way as 
the reduction of the o representation. We first 
determine the representation generated by the p 
orbitals of the external atoms. Since there are 
two such p orbitals per external atom, this 
representation, which will be referred to as the 7 
representation, will have a degree twice that of 
the o representation. The characters of this 
representation are then computed, and entered 
in the character table. This has been done for 
the case of a plane trigonal molecule in Table I. 
The component irreducible representations are 
found as before, and entered in the reduction 
table. 

The condition for the formation of a 7 bond is 
now that there be an orbital of the central atom 
belonging to the same representation as one of 
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TABLE III. Resolution table for linear bonds. 


DIRECTED VALENCE 


Doon Zut Ny Ag Au 
Ss 1 0 0 0 0 0 
p 0 1 0 1 0 0 
d 1 0 1 0 1 0 
o 1 1 0 0 0 0 
cy 0 0 1 1 0 0 


TABLE IV. Resolution table for angular bonds. 


Cop A A: By 


1 0 
Pp 1 1 
d 2 1 1 
o 1 0 0 
T 1 1 1 


the irreducible components of the z representa- 
tion. Since this molecule is already supposed to 
be held together by o bonds, it is not necessary 
that all of the irreducible components of the z 
representation be matched by orbitals of the 
central atom. However, unless at least half of 
the irreducible components of the 7 representa- 
tion are so matched, it will be impossible to 
localize the z bonds, and the resulting molecule 
will be of the type ordinarily written with reso- 
nating double bonds. 

Thus in the plane trigonal case, Table II 
shows that 7 bond formation is possible through 
the p orbital belonging to A2’’(p.), the d orbitals 
belonging to E’’(d,.,dy:) or through the two 
orbitals belonging to E’ which are not used in 
forming the original « bonds. Since, however, 
the « bonds are probably formed by a mixture 
of both the p and the d E’ orbitals, these last two 
7 bonds are probably weaker than the others. In 
general we shall divide the 7 bonds into two 
classes, calling them ‘‘strong”’ if they belong to 
representations not used in ¢ bond formation, and 
“weak” if they belong to those representations 
already used in o bond formation. 

Because of the “resonating’”’ character of = 
bonds, it is usually difficult to form a mental 
picture of them. In this trigonal case which we 
have been discussing, we may imagine the p, 
orbital of the central atom to interact in turn 
with the p, orbitals of the external atoms. If the 
orbitals d,, and d,, are available, these may be 
combined with p., in the same way that p, and p, 
can be combined with s, to form three directed 
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m orbitals. No such simple picture, however, 
seems to be available for the t bonds formed by 
dz, and dzr+y. 

It should be noted that this method does not 
predict directly the type of bond arrangement 
formed by any given electron configuration. 
Instead it merely tells whether or not a given 
arrangement is possible. In many cases it is 
found that several arrangements are possible for 
a single configuration of electrons. In these 
cases the relative stability of the various arrange- 
ments must be decided by other methods, such 
as Pauling’s “strength” criterion, or considera- 
tion of the repulsions between nonbonded atoms. 


RESULTS 


The results of these calculations are most con- 
veniently arranged according to the coordination 
number of the central atom. For the sake of com- 
pleteness all of the results, including those pre- 
viously obtained by Pauling, Hultgren and 
others, are contained in the following summary. 


Coordination number 2 


If the central atom forms bonds with two 
external atoms only two arrangements of the 
bonds are possible: a linear arrangement (group 
Dz.) and an angular one (group C>2,). The resolu- 
tion tables for these are given in Tables III and 
IV. The configurations sp and dp can lead to 
either arrangement. The linear arrangement, 
however, is favored by both the repulsive forces 
and the possibilities for double bond formation 
and is therefore the stable arrangement for these 
configurations. The configurations ds, d? and p* 
on the other hand must be angular. 

In the linear arrangement two p and two d 
orbitals are available for double bond formation, 
while in the angular arrangement only two strong 
m« bonds are possible, one of which must be 
through a d orbital, the other of which may be 
formed by either a d or a p orbital. The other two 
possible x bonds are weakened by the fact that 
their orbitals belong to representations already 
used by the o bonds. 

The angular nature of the p? bonds in such 
molecules as HO and H.S is well known, and 
need not be discussed further. As examples of 
double bond formation in molecules having this 
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TABLE V. Resolution table for trigonal pyramid bonds. 


primary valence structure the nitroso compounds 
are typical. The prototype of these compounds 
is the nitrite ion NO.-. In this ion each oxygen 
atom is joined to the nitrogen by a o bond in- 
volving one of the p orbitals of the nitrogen. In 
addition to these bonds, however, there is a 7 
bond which belongs to the representation B,; and 
cannot be localized. Its component orbitals are 
the p, orbital of the nitrogen atom and the 
normalized sum of the p, orbitals of the two 
oxygen atoms. The remaining parts of the oxygen 
p orbitals are occupied by unshared pairs of 
electrons. Because of the distributed character 
of the 7 bond no single valence bond picture can 
be drawn for this molecule, but only the reso- 
nating pair 


O-—N 
\ 
O O-. 


In the true nitroso compounds RNO, and also 
in the nitrosyl halides CINO, etc., the situation 
is very similar. The lack of complete symmetry 
does not change any of the essential features of 
the structure of these molecules. In these cases, 
however, the difference in electro-negativity be- 
tween oxygen and the atom or group R may 
cause ionic structures to play an important role. 
Pauling has suggested for example that in CINO 
the important structures are 


CI-—N and CI-N 


\ 
O Or. 


The bond formation in both these structures is 
in accord with the present theory. 

Examples of the linear configuration are found 
in the ions of the Ag(NHs)2* type and those of 
the I;- type. In Ag(NHs)2* the configuration is 
certainly sp. In Is~ and its relatives it is not cer- 
tain whether the promotion of one of the 5p 
electrons of the central atom is to the 5d or 6s 
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TABLE VI. Resolution table for three unsymmetrical bonds 


in a plane. 


Coy Be 


orbital, but in either case the resulting configura- 
tion (ps or pd) should give the observed linear 
arrangement. 

The effect of double bond formation is shown 
clearly by COz and SOs. In CO: the valence 
configuration of the carbon is sp*. The primary 
structure of o bonds is sp, which requires the 
molecule to be linear. The other two orbitals 
form the two x bonds. In SOs, on the other hand, 
the primary bonds are certainly formed by two 
p orbitals, thus producing angular arrangement. 
The customary way of writing the structure of 
this molecule, 


(analogous to NO-~) is quite possible, since the 
a bond can be formed by the third p orbital of 
the sulfur, but the structure 


O=S 


N 
O 


is also possible and probably as important as the 
first. In this structure one of the d orbitals of the 
sulfur acts as an acceptor for one of the electron 
pairs of the oxygen. It should be noted that the 
formation of two double bonds by an atom does 
not require the linear arrangement. 


Coordination number 3 


For this number of o bonds three arrangements 
are important: the plane trigonal arrangement 
(group D3,) already discussed, the trigonal pyra- 
mid (group C3), and an unsymmetrical plane 
arrangement with two of the bond angles equal, 
but not equal to the third (group C2,). The reso- 
lution tables for the last two are given in Tables 
V and VI. As has already been shown, the plane 
arrangement is stable for the configurations sp’, 
sd*, dp’, and d’. (The pyramidal configuration is 


|| 
Cay Ai Az E 
1 0 0 
p 1 0 1 s 1 0 0 0 
d 1 0 2 p 4 0 1 1 
o 1 0 1 d 2 1 1 1 
1 1 2 o 0 1 0 
O=St+ 
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also possible, but less stable because of repulsive 
forces. Such cases will simply be omitted in the 
future.) The configurations p* and dp are there- 
fore the only ones leading to a pyramidal struc- 
ture.!° The configuration dsp leads to the un- 
symmetrical plane arrangement. 

The possibilities of double bond formation in 
the plane trigonal arrangement have already 
been discussed. In the pyramidal arrangement 
no strong double bonds are possible, although 
weak 7 bonds are possible with d orbitals. Thus 
the sulfite ion, for example, is restricted to the 
structure 


It is interesting to contrast this structure with 
that of SO3. The most easily removed pair of 
electrons in SO; is the pair of 3s electrons in the 
sulfur atom. If these are removed, however, one 
of the 3p pairs forming a o bond falls to the 3s 
level, thus making the valence configuration sp? 
and the molecule plane. The vacated 3p orbital 
is then filled by an unshared pair from the 
oxygen atoms to give a a bond. Two more 
unshared pairs of electrons from the oxygen 
atoms may interact with the 3d orbitals thus 
forming two more z bonds, which are somewhat 
weaker than the first. The SO; molecule has 


TABLE VII. Reduction table for tetrahedral bonds. 


T2 


TABLE VIII. Reduction table for tetragonal plane bonds. 


Aig Bu By Bou Eu 


therefore three comparatively strong double 
bonds, while the SO; ion has only weak double 
bonds, if any at all. It is this lack of stabilizing 
double bonds which gives SO;~ its relative 
instability. 

The remaining configuration, dsp, should give 
rise to three bonds in a plane forming two right 
angles. No examples of molecules of this con- 
figuration are known. This is hardly surprising, 
however, in view of the instability of this ar- 
rangement compared to the others. If molecules 
of this arrangement are to be found at all they 
would be complex ions of the transition elements. 
For example, if the ion Ni(NHs;)3++ existed it 
would be of this structure. 


Coordination number 4 


With a coordination number of four there are 
three arrangements which need consideration. 
The reduction tables for these are given in 
Tables VII, VIII, IX and X. Table VII is for 
the regular tetrahedral arrangement of the bonds. 
It is easily seen that the configurations sp* and 
d’s lead to this arrangement. For double bond 
formation there remain two d orbitals which can 
form strong 7 bonds and two other d orbitals 
(if the o bonds are sp*) or two p orbitals (if the « 
bonds are d*s) which can form weaker 7 bonds. 
These possibilities of double bond formation are 
in accord with Pauling’s suggestion that the 
structure of ions of the type XO, should be 
written with double bonds. Thus for the sulfate 
ion the primary o structure arises from the con- 
figuration sp*, but two of the unshared pairs of 
the O- ions can be donated to the vacant d 


TABLE IX. Reduction table for tetragonal pyramidal bonds. 


Ao Bz 


."° The bond angles of 90° in this structure are not pre- 
dicted by the reduction table alone, but can be easily 
found by the process of forming the actual orbitals as 
shown in the second section. 
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orbitals of the sulfur atom, thus giving the 
structure 


O- O 
SN 

O- O 

It must, however, be doubtful that two further 

m bonds are formed to give the structure 


O O 


o” 


The situation in such molecules as SiCl, is 
similar. 

From Table VIII it is seen that tetragonal 
_ planar (square) bonds are possible with con- 
figurations dsp? and d*p*. Beside the primary 
bonds, four strong z bonds are possible, using one 
p and three d orbitals. In Ni(CN).4-, for example, 
the primary valence configuration of the Ni is 
dsp”, so that the structure 


N 


Cc 

II 

N 
is one possible structure for this ion. By a shift of 
one of the 7 pairs forming the C—N triple bond 
to the nitrogen, however, the carbon atom is left 
with an empty orbital, which can accept an elec- 
tron pair donated by the nickel atom. In this 
way three 7 bonds can be formed between the 
nickel and carbon. At the same time one of the 
triple bond 7 pairs can shift in the opposite 
direction and be donated by the carbon to the 
empty » orbital of the nickel. These shifts would 
lead to the structure 
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This structure affords an interesting contrast 
with that of Ni(CO),; in which one more pair of 
electrons must be accommodated. This pair 
occupies the d orbital of the nickel, thus pushing 
the valence configuration up to sp* and making 
the molecule tetrahedral. There are then only 
two 7 bonds possible, and the structure is 


The configuration d?p* is found in such ions as 
ICl,-, which are known to be plane. 

From Table IX it appears that the configura- 
tions dp* and d*p can lead to a tetragonal pyra- 
mid structure, while from Table X it appears 
that the “irregular tetrahedron” is also possible 
for these same configurations and also for the 
configuration d’sp. By ‘‘irregular tetrahedron” is 
meant a structure in which three of the bonds are 
directed to the corners of an equilateral triangle 
and the fourth along the line perpendicular to the 
triangle at its center. The central atom is not 
necessarily located at the center of the triangle, 
but it is probably somewhat above it in the direc- 
tion of the fourth bond. The choice between 
these two structures in the cases of the con- 
figurations dp* and d*p is very close, with the 
advantage somewhat on the side of the irregular 
tetrahedron. The configuration d* can only have 
the pyramidal structure. 

Examples of these configurations are rare. In 
the tetrahalides of the sulfur family we have the 
configuration p*d produced by the promotion of 
one of the p electrons of the central atom to a d 
orbital. Unfortunately the spatial arrangements 
of these molecules are unknown. The other two 
configurations seem to be unstable. Such ions as 
Fe(NH3),** and Fe(CN),.= would have the con- 
figuration d’sp, but do not seem to exist, which 
is perhaps some indication of this instability. 


Coordination number 5 


For this coordination number four bond ar- 
rangements are possible. In Table XI is given the 
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TABLE XI. Reduction table for trigonal bipyramidal bonds. 
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reduction table for a trigonal bipyramid, which is 
seen to be stable for the configurations dsp* and 
dsp. In Table XII the bonds are directed from 
the center to the corners of a square pyramid. 
This is the arrangement expected for the con- 
figurations d*sp?, d*p®, and d‘p. Table XIII 
is for the case of five bonds in a plane, directed 
toward the corners of a regular pentagon, the 
arrangement for the configuration d*p?, and Table 
XIV is for the case of five bonds directed along 
the slant edges of a pentagonal pyramid, the 
arrangement for the configuration d°. 

In PCI; and other molecules of this type, the 
pentavalent state is formed by the promotion of 
an s electron to the vacant d shell. The valence 
configuration is therefore dsp* and the bipyra- 
midal structure is to be expected, in agreement 
with the results of electron diffraction studies. 
Table XI shows that two of the remaining d 
orbitals are capable of forming strong bonds; 
the other two, two weak 7 bonds by accepting 
pairs of electrons. Neglecting weak bonds, the 
structure of PCI; is best written 


cr Cr 


In the molecule Fe(CO); the valence configura- 
tion is again dsp’, and the bipyramidal structure 
is to be expected. Here 7 bonds can be formed by 
donation of pairs of d electrons from the iron 
through the carbon to the oxygens, as in Ni(CO), 
giving the structure 


Ot O 

\ 

4 | 
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TABLE XII. Reduction _ Ef five tetragonal pyramidal 
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TABLE XIII. Reduction table for pentagonal plane bonds. 


Dsh Ay’ Ay” A A” Ey E;” 
s 1 0 0 60 6 606 60 
rp 1 0 oO 1 6 0 0 
d i 0 0 0 0 1 1 oO 
1 0 0 0 i 
r 0 1 i o 1 1 1 1 
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ooo 


In IF;, however, the valence configuration is 
pd’, and the molecule should have the square 
pyramid structure. Here again two strong 7 
bonds can be formed by donations of pairs of 
electrons from the F to the I, so that the struc- 


ture may be 
F 
Ft 


F 

F 
but the high electro-negativity of fluorine makes 
the existence of the double bonds doubtful. The 
spatial configuration of this molecule has not 
yet been determined experimentally. 

The pentagonal configurations crowd the 

atoms so much that they must be unstable. 


Coordination number 6 


Six bonds may be arranged symmetrically in 
space in three ways: to the corners of a regular 
octahedron, to those of a trigonal prism, or to 
those of a trigonal antiprism (an octahedron 
stretched or compressed along one of the tri- 
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TABLE XV. Reduction table for octahedral bonds. 


On Ary Am Eg Eu Tig To Tru 


s i 0 0 0 0 0 0 0 0 90 
rp 0 0 0 0 0 0 0 1 0 0 
d 0 0 0 0 1 0 0 0 1 90 
1 0 0 0 1 0 0 1 0 90 
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TABLE XVII. Reduction table for trigonal antiprismatic 
bonds. 
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Oe 


gonal axes). The reduction tables for these possi- 
bilities are given in Tables XV, XVI, and XVII. 

Table XV shows that octahedral bonds are 
formed by the configuration d*sp* and no other. 
The commonness of this arrangement is therefore 
due to the fact that the configuration d?sp* is the 
usual configuration of six valence electrons, 
rather than any particular virtue of the octa- 
hedral arrangement. This configuration arises 
both in the 6-coordinated ions of the transition 
elements, and in the molecules of the SF. type 
in which one s and one ? electron are promoted 
to the next higher d level. Table XV also shows 
that all three remaining d orbitals are capable of 
forming strong 7 bonds. In the ferrocyanide ion 
Fe(CN),~*, these d orbitals can donate their 
pairs of electrons to the nitrogen atoms, giving 
Pauling’s structure 


al 
| 
N-=C=Fe-—C=N 
/ || 
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If SF, the remaining d orbitals are empty, and 
can accept pairs of electron from the fluorine 
atoms, giving the structure 


F+=S*—F 
F Ft 


It is this structure which accounts for the great 
resistance of SF, to hydrolysis. It has been sug- 
gested by Sidgwick" that the hydrolysis of 
halides takes place either by their accepting a 
pair of electrons from the oxygen of a water 
molecule, with subsequent decomposition, or by 
donation of an unshared pair on the central atom 
to a hydrogen atom, with loss of hypohalous 
acid. With the usual single bonded structure of 
SF, the first mechanism is possible, and one 
would expect the hydrolysis to take place easily. 
With the double bonded structure, however, all 
of the electrons in the outer shell of the sulfur 
are taking part in bond formation, so that one 
should expect the same inertness as that dis- 
played by CCl. It is interesting to note that 
SeF; is like SFs, but TeFs is easily hydrolyzed. 
This must be due to the possibility of the Te 
atom’s accepting a pair of oxygen electrons in the 
vacant 4f orbital. 

In MoSe and WS: in the crystalline form the 
valence configuration is d‘sp and the arrangement 
is prismatic, as Hultgren‘ has pointed out. Table 
XVI indicates that the two empty # orbitals of 
the metal atom cannot form strong 7 bonds to 
the sulfur atoms. The configuration d°p should 
also lead to the prismatic arrangement, but no 
examples are known. 

The cases of the ions SeBre~? and are 
interesting in that the configuration in each of 
these should be #*d*. This configuration should 
lead to the antiprismatic arrangement, i.e., the 
octahedral symmetry should not be perfect. If, 
however, the unshared pair of s electrons is 
promoted to a d orbital, and one of the valence 
pairs slips into its place, the configuration will 
be the octahedral d?sp*. The deciding factor here 
may be the possibility of double bond formation 
offered by the octahedral but not by the anti- 
prismatic arrangement. 


 Sidgwick, The Electronic Theory of Valency (Oxford, 
1927), p. 157. 
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TABLE XVI. Reduction table for trigonal prismatic bonds. 
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The configurations d’sp?, d's and d‘p? do not 
permit any arrangement in which all of the bonds 
are equivalent. While the configuration d*sp? 
might form bonds of a mixed type, e.g., the tetra- 
hedral d*s+the angular p*, the bonds will be 
weak, and the configuration unstable. Molecules 
which might be expected to have this arrange- 
ment, e.g., Fe(CN),-*, rearrange the electrons to 
obtain the more stable d*sp* configuration. In 
Fe(CN),-* for example, the odd electron, which 
normally should occupy a 4 orbital, instead goes 
to a 3d orbital making the third 4 orbital 
available for the formation of octahedral bonds. 


Coordination number 7 


The coordination number 7 is extremely rare. 
The fact that it appears only in the heavier 
atoms, such as Zr, Cb, I, and Ta, leads one to 
suspect that for stability f electrons are neces- 
sary, although it is possible that the determining 
factor is ion size. Two arrangements of seven 
bonds have been observed :” the ZrF;~ structure, 
which may be obtained from the octahedron by 
adding an atom at the center of one face; and 
the TaF,~ structure, which may be obtained 
from the trigonal prism by adding an atom at the 
center of one of the square faces. The reduction 
tables for these two arrangements are given in 
Tables XVIII and XIX. 

From Table XVIII we see that if no f orbitals 


TABLE XVIII. Reduction table for ZrF;-* type bonds. 
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Ss 1 0 0 0 
p 1 0 1 1 
d 2 1 1 1 

f 2 1 2 2 
o 3 1 2 1 
T 3 3 4 4 


® Hampson and Pauling, J. Am. Chem. Soc. 60, 2702 
(1938); Hoard, J. Am. Chem. Soc. 61, 1252 (1939). 
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TABLE XX. Reduction table for cubic bonds. 
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TABLE XXII. Reduction table for dodecahedral bonds. 


Va Al Ae Bi Bz E 
s 1 0 0 0 0 
p 0 0 0 1 0 
d 1 0 1 1 1 
¢ 2 0 0 2 2 
. 2 2 2 2 4 


TABLE XXIII. Reduction Fas for face-centered prismatic 


— 


are involved, the ZrF;-* type bonds can arise 
from the configurations d*sp* and d°sp. The con- 
figuration expected for ZrF;-* is d®sp, so that it 
is not necessary to appeal to the f orbitals to 
form these bonds. This, however, does not pre- 
clude the possibility that f orbitals are used to 
strengthen the bonds. 

Table XIX shows that the TaF;~ structure is 
possible with the configurations d’sp', d‘sp?, d°sp, 
d‘p* and d*°p?. The ion TaF;~ itself is isoelectronic 
with ZrF;-*, so that it is difficult to understand 
why it should prefer its structure to that of 


‘ZrF;-. It is quite possible that the bonds in all 


the 7-coordinated molecules are so ionic that the 


197 

TaBLE XXI. Reduction table for tetragonal antiprismatic 
‘bonds. 
Cry Ai Bz 
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TABLE XIX. Reduction table for TaF 7 type bonds. 
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TABLE XXIV. Summary of stable bond arrangements and 
multiple bond possibilities. 


COORDINATION 
CONFIGURATION 


ORBITALS 


ARRANGEMENT 


linear 

linear 

angular 

angular 

angular 

trigonal plane 
trigonal plane 
trigonal plane 
trigonal plane 
unsymmetrical plane 
trigonal pyramid 
trigonal pyramid 
tetrahedral 
tetrahedral 
tetragonal plane 
tetragonal plane 
irregular tetrahedron 
irregular tetrahedron 
irregular tetrahedron 
tetragonal pyramid 
bipyramid 
bipyramid 
tetragonal pyramid 
tetragonal pyramid 
tetragonal pyramid 
tetragonal pyramid 
pentagonal plane 
pentagonal pyramid 
octahedron 

trigonal prism 
trigonal prism 
trigonal antiprism 
mixed 

mixed 

mixed 

ZrF 73 

ZrF 7-3 

TaF7? 

TaFy? 

TaF7? 
dodecahedron 
antiprism 
face-centered prism 


| | SE | 


v's 


— 
> 


d? 
d 
d 
d 
d 
pd? 
a3 


directed nature of covalent bonds has little im- 
portance in determining the atomic arrangement. 
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Coordination number 8 


Until very recently the spatial arrangement of 
no molecule involving a coordinate number of 8 
had been determined experimentally. It had been 
supposed that the arrangements which were most 
probable were the cubic and the tetragonal anti- 
prismatic. Hoard and Nordsieck" have now de- 
termined the arrangement of the ion Mo(CN) s~* 
and have found that the structure is neither cubic 
nor antiprismatic, but instead that of a dode- 
cahedron with triangular faces and symmetry 
V.(D2*). The reduction tables for these three 
cases are given in Tables XX, XXI, and XXII. 
Table XXIII is the reduction table for an ar- 
rangement which is derivable from a trigonal 
prism by placing two atoms at the centers of two 
of the rectangular faces. It is seen that the cubic 
requires the use of f orbitals, and even then arises 
only from the configurations d*fsp* and d°f's. The 
cubic arrangement should therefore not be or- 
dinarily found. The antiprismatic arrangement is 
stable for the configurations d‘sp* and d*p*, and 
the dodecahedral for only d‘sp*. The observed 
arrangement for Mo(CN) 5, which has the con- 
figuration d‘sp*, indicates that the dodecahedral 
arrangement has greater stability than the anti- 
prismatic. The structure of such ions as TaFs~* 
and the molecule OsFs which should have the 
configuration d°sp? should be the face-centered 
prismatic structure of Table XXIII. 

Table XXIV summarizes these results. In it 
are given the most stable bond arrangements for 
each configuration of two to eight electrons in 
s, p or d orbitals. The orbitals available for 
strong z bonds and weak z bonds are also given. 
In those cases where there is a choice of two or 
more orbitals when only one can be chosen, the 
orbitals are enclosed in parentheses. 


13 Hoard and Nordsieck, J. Am. Chem. Soc. 61, 2853 
(1939). 
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The theory of the first paper of this series is extended to include reactions at electrode 
surfaces, The resulting expression for the rate of an electrode reaction includes both activation 


controlled and diffusion controlled reactions. The problem of the effect of stirring rate on 
reaction rate is discussed. The results are applied to the theory of reactions at ‘‘floating”’ 
electrodes, e.g. the solution of metals in acids. It is shown that the reaction rate expressions at 
such electrodes are of a very complex form, but can be approximated in special cases. It is ~ 
found that the rate of such a reaction will in general depend not only on the concentration of 


the principal reactants, but also on the concentration of every ion in the solution. 


NE of the commonest types of chemical 
reaction takes place between a metal and a 
solution containing ions. As examples there. are 
the processes of electrolytic deposition and solu- 
tion, the solution of metals in acids, the displace- 
ment of one metal by another, and corrosion of 
metals. Each of these processes can be regarded 
as consisting of processes of formation of ions 
or discharge of ions at a metal surface. It is 
therefore appropriate to include them all in the 
category of electrode reactions. 

As a metal atom leaves the electrode surface 
it passes through three regions: first, the surface 
of the metal, on which there may or may not be 
a layer of adsorbed atoms; second, a layer of 
adsorbed ions; and third, a diffusion layer in 
which the ionic concentration gradually changes 
from that near the layer of adsorbed ions to 
that of the main body of the solution. These 
layers are shown schematically in Fig. 1. 

The rate at which metal atoms go into solution 
as ions, or vice versa, is therefore determined by 
two factors: the rate of formation of ions in the 
adsorbed layer and the rate of passage of these 
ions through the diffusion layer. Both of these 
processes have been discussed at great length in 
the literature. 

The principal difference between these reac- 
tions and other heterogeneous reactions is caused 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

‘For reviews see Glasstone, Annual Reports of the 
London Chemical Society 34, 107 (1937); Bowden and 
Agar, ibid. 35, 90 (1938). 
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by the electrical effect of the layer of adsorbed 
ions and the electrical image of this layer in the 
electrode. This electrical double layer leads to a 
difference in potential between the solution and 
the electrode. It is, of course, impossible to 
measure absolutely the value of this potential 
difference? but nevertheless it is possible to 
measure changes in it. The effect of this potential 
on the reaction rate can be seen from Fig. 2. In 
this figure is given the potential energy of an 
atom as a function of its distance from the 
electrode surface. The curve ABCD is the curve 
for the potential energy of the neutral atom, 
while the curve ECFG is that for the ion. At the 
crossing point C quantum-mechanical resonance 
rounds.off the crossing so that the two curves no 
longer intersect. The actual charging of the ion 
will therefore take place at this point. The energy 
difference between the point C and the point B is" 
essentially the activation energy of the charging 


Electrode 
Layer of adsorbed atoms 
Layer of adsorbed ions 
Diffusion Layer 
Solution 


Fic. 1. Diagram of layers at an electrode surface. 


Guggenheim, Modern Thermodynamics (Methuen, 


1933), Chapter X. 
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Fic. 2. Potential energy curves of adsorbed ions and atoms 


process. The energy at C and at B, however, is a 
function of the potential difference between the 
solution and the electrode. The activation energy 
of the charging process is therefore also a 
variable quantity and will change with the 
potential difference. 

According to the theory of absolute reaction 
rates* the rate of formation of ions in the ad- 
sorbed layer from the metal surface will be 
given by 


h 


where AF;* is the free energy difference between 
an atom at B and one at C when the potential 
drop is zero, Vc and Vz are the electrical poten- 
tials at C and B, [A ] is the activity of the atoms 
A in the’surface or adsorbed on it, and « is the 
charge of the ion. Letting 


_k 


this rate may be written in the form 
Ri[A Jem 
The total current from the electrode due to this 
process is then given by 
Iq! = eki[A 
In a similar way the current (still measured as 


flowing from the electrode) caused by the dis- 
charge of ions is 


* Eyring, J. Chem. Phys. 3, 107 (1935); Eyring and 
Wynne-Jones, ibid. 3, 492 (1935); Kimball, ibid. 6, 447 
{f939)° Eyring, Glasstone and Laidler, ibid. 7, 1053 
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where [At ]aas is the activity of the ions in the 
adsorbed ion layer and k2 is a constant. The net 
current J,, flowing from the electrode, is there- 
fore given by 


I,= A J 
— A+ (1) 


The rate at which ions migrate through the 
diffusion layer will follow the Nernst‘ law. 


D 


where D is the diffusion constant, / is the thick- 
ness of the diffusion layer, and [At] is the 
activity of the ions in the main body of the solu- 
tion. The corresponding current, J,’’’, is given by 


eD 
I (2) 


Except for an extremely short period at the 
beginning of an experiment, during which the 
adsorbed ion layer is established, these two 
currents must be equal, for otherwise the ad- 
sorbed jon layer would not come to equilibrium. 
The activity [A*]nas is therefore determined by 
the equation J, =J,’’’. Setting these two equal and 
solving for [A+ ],as immediately gives 


kil A 4. (D/) [At] 
(Vo-V PRT 4 /] 


[A*]oas= 


Substitution of this equation in (2) gives for the 
current when the adsorbed ion layer is at 
equilibrium 


ko 


where 


= 
D)/] 


If we adopt the arbitrary convention that the 
potential of the solution is zero, and assume 
furthermore that Vg=V, the potential of the 
electrode, Vr =0 (neglect of the diffusion poten- 
tial), and put Vc=aV, Eqs. (3) and (4) can be 


4 Nernst, Zeits. f. physik. Chemie 47, 52 (1904). 
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put in the form 
k 
(5) 
2 


(6) 


We may simplify these equations still further if 
we note that at the equilibrium potential V4, 
I,=0 if [A]=[A+*]. This condition gives im- 
mediately 


alk? 


and 


The equation for \ can be put in the form 
Va) /kT 
4. 0) 
where and Q=D/l, or 
kT 44 /Q, (8) 


This equation may be simplified in the two 
limiting cases 


O> Pee and Qc 


In the first case the reaction is activation con- 
trolled, and approximately 


A= (8a) 
I,= eP 4 | 
(7a) 


In this case the kinetics are essentially those of 
Eyring, Glasstone and Laidler.* In the second 
case 


A=Q, (8b) 


Ta= A J—[A*]}. (7b) 


It is interesting to note in this case that the rate 
is independent of the value of a. 

The transition from activation control t 
diffusion control comes when 


Q= (9) 


If the electrode process involves the charge or 
discharge of a positive ion the reaction will be’ 
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diffusion controlled for V>V, and activation 
controlled for V&V,, while for negative ions the 
reverse is true. It is quite possible that when 
several electrode processes are taking place at 
one electrode some may be diffusion controlled 
while others are simultaneously activation con- 
trolled. 

The effect of the stirring rate on these reactions 
is essentially the effect of the change of the thick- 
ness of the diffusion layer. If we assume simply 
that the thickness of this layer is inversely pro- 
portional to the stirring rate, ¢, then Q will be 
of the form 


Q=0Qo, (10) 


where Qo is a constant. The form assumed by 
Eq. (7) is then 


(11) 


where K, and Kye are independent of oc. The 
form of this equation is shown in Fig. 3. This 
equation is of the same form as the Langmuir 
adsorption isotherm, and has the same relation- 
ship to the usual form of stirring curve 


I=Ac" 


as the Langmuir isotherm does to the Freundlich 
isotherm. Thus the suggestion of Nernst and 
Merriam’ and of Brunner*® that the rate is pro- 
portional to the two-thirds power of the stirring 
rate must be regarded as open to doubt, for their 
data may be fitted about as well by Eq. (11). 
With the exception of overvoltage measure- 
ments, which will not be discussed here, most of 
the experimental data on rates of electrode 
reactions have been obtained from experiments 
with what might be called a ‘‘floating”’ electrode, 
i.e., a piece of metal simply immersed in a 


I 


Fic. 3. The effect of stirring rate on solution rate. 


5 Nernst and Merriam, Zeits. f. physik. Chemie 53, 235 


* Brunner, Zeits. f. physik. Chemie 47, 56 (1904). 


|| 
e 
ot 
| 
he 
by 
(2) 
the 
the 
— 
| by 
|_| 


202 


solution with no control over its electrical poten- 
tial. In such a situation there are always at least 
three reactions which may take place: formation 
of metal ions in the solution from the metal, 
discharge of positive ions on the metal, and 
discharge of negative ions on the metal. It is a 
characteristic of electrode reactions that each of 
these processes tends to produce a potential in 
the metal of such a sign as to slow down the 
reaction producing it. The first and last of these 
reactions tend to produce a negative potential 
on the metal, while the second tends to produce 
a positive potential. After a very short time, 
therefore, the potential of the electrode will 
reach such a value that the positive current 
from the metal caused by reactions of the first 
and third types is exactly equal to the negative 
current caused by reactions of the second type. 
When this potential is once attained, the reac- 
tions proceed with no net current, so that this 
potential is maintained. 

To illustrate this process, let us consider a case 
in which the discharge of negative ions is 
negligible, one in which a metal A is dissolving 
to form singly charged positive ions A* and 
simultaneously singly charged positive ions Bt 
are discharging to form a second metal B. 
Since reactions of this type are usually diffusion 
controlled, we shall also assume that the currents 
due to these reactions are given by equation 
(7b). We shall also assume that the back re- 
actions (discharge of A+ and solution of B) can 
be neglected. With these simplifying assumptions 
the currents due to the two reactions are 


Ta= (12) 

— Q,[ (13) 

The steady-state potential is then determined by 
A] =0,[ Bt], 


so that 
(14) 
Qa [A] 
If this is substituted into (12) we find 
I,= «Q,[ Bt]. (15) 


In other words the potential of the electrode 
adjusts itself to just such a value that the rate 
of solution of A is exactly equal to the rate at 
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which B+ ions can pass through the diffusion 
layer. 

It is obvious that this result is oversimplified, 
for it gives a rate of solution which is independent 
of the metal dissolving, the same, for example, 
for sodium as for platinum. The difficulty is 
obviously in the neglect of the back reactions, in 
particular the solution of B, for the activity of 
the discharged B atoms rises very rapidly. If the 
volume of the solution is large, however, the 
discharge of A+ can be safely neglected, so that 
only Eq. (13) needs to be modified. If we include 
the back reaction in this equation it becomes 


Th= B]—[Bt]}. (16) 
The equation determining V is now 
A B]—0,[ B+] =0 
or 
Bt] 


Now if B is more noble than A, V, will be larger 
than V, and unless Q, is much smaller than Q,, 
the term 


17) 


will be small compared to 
[A 


Then we will have approximately 


Q. [A] 
I,= Bt ] (18) 


as before. On the other hand, if B is less noble 
than A, the term 


will be large compared to 
[Bt] 
[B] 


that is, the potential becomes the equilibrium 
potential of the B electrode and the solution of 


eValkT 


and 


and 


— 


=o 
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Qn, 


(18) 
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potential: 


Since in this case Vs Va, the solution proceeds 
very slowly. The net effect is therefore a slow 
corrosion of the electrode. 

To illustrate the effect of negative ions let us 
consider an electrode consisting of a metal A in 
contact with a solution containing positive ions 
Bt and negative ions C~. The three currents 
from the electrode are then 


B+]}, 
—[C-]}. 

The steady-state potential is therefore given by 


ax+b—c/x=0, (20) 
where 


c= 
The solution of (20) is 
—b+(b?+4ac)! 


x= ’ 
2a 


where, since x must necessarily be positive, only 
the positive sign before the ea is significant. 
Hence 


(b?-+4ac)! 
x= ; (21) 
2a 


We shall distinguish three cases. First suppose 
b>(4ac)!. This corresponds to a reaction in 
which the fastest reaction is the discharge of the 
negative ions. In this case (Case I) we may 
approximate Eq. (21) by 


—b+b(1+ (2ac/b*)) 
2a 


x= 


=c/b. 
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A proceeds at the rate determined by this 


To a first approximation 


b=Q.[C J, 
so that 
[C] 
c= eV clkT 
and 
V=V.+— In 
[c-] 


In other words the potential is the equilibrium 
potential for the discharge and charge of C-. 
Then, if J,’ is the solution current of A, 


MI 


(22) 
[e-] 
the condition for this case is 
X A B 5/47} 


which will usually hold if V.>V., Vi>V., and 
Q.>Qa+Q>. Because of the negative value of 


e =€ 


' the exponent in Eq. (22) reactions of this type 


are necessarily slow. 
Case II is that in which |b|<«(4ac)', in which 
case Eq. (21) becomes approximately 


x=(2c/a)'. 


This case in turn may be subdivided into two 
subcases, depending on whether A is more or 
less noble than B. If A is less noble than B 
(Case Ila), a is given approximately by 


a=Q,[A 
and x and J,’ become 
I,’ =€(20.0.[A te" (23) 


On the other hand, if A is more noble than B 
(Case IIb), 


Q.[B]. 
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The condition for this case is 


which may be expected to hold if V.>V. or 
V.> Eq. (23) will therefore hold if 
or and Eq. (24) when V.>V.>V, 
or Va>V> Vi. 

Case III arises when b<&-—(4ac)!, when 
approximately 

x=—b/a 

and 


Ta! = A 


If Va<Vz (Case IIIa) this is approximately 
Tq’ = (26) 
If Vk<Va (Case IIIb) Eq. (25) approximates to 


X (27) 


The conditions for this case are V.>V., Vi 
QatQ>Q.. If Var>Vi>V. either Eq. (22) or 
Eq. (27) will hold, depending on whether 


Q.>>Q.+Q or whether If Vi>V.> V- 
Eq. (22) or Eq. (26) will hold. 

It should be noted that the conditions given 
apply only when [C] is appreciably large. If [C] 
is extremely small, Case II becomes impossible, 
and the kinetics will belong to Case I or Case III 
according to the relative magnitudes of Qu, Qo 
and Q,. Case I and Case IIIb will give very slow 
rates under these conditions, which include most 
corrosion problems. The reactions of metals with 
acids and similar reactions belong principally to 
Case IIIb. 

The complexity of these rate expressions indi- 
cates very strongly that the experimental study 
of electrode reactions by means of experiments 
with “floating” electrodes is apt to be quite 
unreliable. The rates obtained for the solutions 
of metals in acids, for example, show very com- 
plicated behavior when the negative ion is 
changed and when salts are added to the solution. 
In part these effects are due to changes in 
activity coefficients, but the results derived above 
show that there is always a direct effect due to 
discharge of the negative ions. 

It would seem best, therefore, to conduct 
experiments with controlled potential between 
the electrode and the solution. If the above 
theory is correct, the rate of solution of a metal, 
whatever the solution used, should depend only 
on the potential difference between the electrode 
and the solution. Experiments are now being 
carried out in this laboratory to investigate the 
truth of this assumption and to obtain measure- 
ments of the constants P, Q and a. 
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A theory of phase changes due to hindered molecular rotation in crystals is developed on 


the basis of classical statistical mechanics. It is demonstrated that local hindrance of relative 
rotation of neighboring molecules in the lattice can produce a non-uniform distribution in 
orientation below a certain critical temperature. The theory provides a semi-quantitative 
description of the behavior of the hydrogen halides in certain of their lambda-transitions, 
with a relative torque potential of the form (¢/2) cos y for a pair of neighboring molecules, 
the axes of which are inclined at an angle + to each other. The barrier height ¢ is estimated to be 


between 200 and 300 cal./mole for hydrogen iodide and hydrogen bromide. 


NOMALIES are frequently observed in 
the heat capacity-temperature curves of 


crystals containing molecular or ionic groups 


capable of rotation in the lattice. The relation- 
ship of such anomalies to hindered molecular 
rotation has been established by Pauling.' Since 
they generally occur not as smooth maxima but 
in the form of lambda-points with a discon- 
tinuity on the high temperature side, it must be 
concluded that they are associated with phase 
changes of the second kind. In order to account 
for phase changes either .of the first or second 
kinds on the hypothesis of hindered molecular 
rotation, it is necessary to suppose that the 
hindering torque acting on a molecule in the 
lattice depends upon the state of rotation of its 
neighbors. Processes of this type have been called 
cooperational phenomena by Fowler. Fowler? has 
developed a cooperational theory of the freezing 
of molecular rotation in crystals which provides 
a qualitative description of the observed phase 
changes. However, since Fowler’s theory con- 
tains certain arbitrary features, it seems de- 
sirable to formulate a theory based more directly 
on the principles of statistical mechanics and in 
which the detailed role of molecular interaction 
is made more clear. It is the purpose of the 
present article to describe such a theory, carried 
to a stage of approximation equivalent to that 
of the Bragg-Williams* theory of super-lattices in 
alloys. In place of the Bragg and Williams long 


' Pauling, Phys. Rev. 36, 430 (1930). 
*Fowler, Statistical Mechanics (Cambridge, 1936), 
Second edition, Chap. XXI 


*Bragg and Williams, Proc. Roy. Soc. 99 
(1934) and Williams, Proc. Roy. Soc. A145, 6 
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range order parameter, a set of functions of 
molecular orientation appear, each of which 
specifies the distribution in orientation of the 
molecules located on one of the interpenetrating 
lattices of the crystal. By straightforward 
methods, the conditions for a transition from a 
state of uniform to non-uniform distribution in 
molecular orientation are formulated, and ex- 
pressions for the thermodynamic functions of the 
crystal in terms of the distribution functions are 
presented. From this point of view, the transition 
manifests itself as a segregation of the molecules 
from a random distribution in orientation to 
regions of preferred orientation, rather than as a 
change from rotating to librating states as in the 
Fowler theory. Local hindrance of the relative 
rotation of neighboring molecules in the lattice 
is assumed permanently present, and the exist- 
ence of a transition temperature is demonstrated, 
below which local hindrance becomes effective in 
producing a non-uniform distribution in orienta- 
tion of the crystal as a whole. The applicability 
of the general method of the present theory is 
not limited to phase changes associated with 
molecular rotation. Work is now in progress on a 
theory of fusion in which the method is applied 
to the lattice vibrational motion of a crystal. 
We consider a crystal composed of m inter- 
penetrating simple lattices, each site of which is 
occupied by one of a set of N identical molecules 
possessing rotational degrees of freedom. The 
rotational configuration of the entire crystal 
may be partially characterized by a set of 
distribution functions gi(w), -gn(w), so defined 
that (N/n)gi(w) specifies the density of molecules 


206 


on lattice k in orientation w. To an arbitrary 
specification of the functions g,(w), subject to 
normalization, 


f (1) 


there corresponds a subset of the entire set of 
rotational states accessible to the crystal. If the 
fluctuation in energy among the states of the 
subset is neglected, the entropy S corresponding 
to the specified distribution functions g;(w) may 
be calculated as k log 2 where Q is the total 
number of states in the subset. Assigning the 
quasi-classical statistical weights dpwidw;/h" in 
the rotational momentum-configuration space of 
each molecule 7 of 7 rotational degrees of freedom, 
we calculate 


S/Nk=—— | gelog gudw+S'/k, (2) 


nN k=1 


where S’ is the nonconfigurational rotational 
entropy of a free molecule. We shall henceforth 
suppress .S’ since it plays no role in the transitions 
with which we are concerned. The average 
potential energy E, calculated with a similar 
neglect of fluctuations in energy in the subset of 
configurations corresponding to the specified 


gx(w), is 


1 n 
E/N=— 


Qn k=1 l=1 


J J VO (co, 
(3) 


Nin 

0") =D 0"), 
b=1 
where V.,“" is the mutual potential energy of 
two molecules, respectively situated on site a 
of lattice k in orientation w and site b of lattice / 
in orientation w’. In the same approximation the 
configurational free energy F may be written 


1 n 
=— fe log gidw 
nN k=1 


(4) 


1 ni on 
2nkT 


Extremalization of the free energy, by means of 
the calculus of variations, with respect to the 
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functions gi(w), subject to the normalization 
restraint, leads to the following set of integral 
equations 


(5) 


i=1 
Kurlw, w’) = — Vw, w’) /RT. 


The form of the solutions of these equations is 
determined by the kernels Ky, which in turn 
depend upon the form of the V“(w,w’) as func- 
tions of molecular orientation. As we shall 
demonstrate, interactions of certain types lead 
to phase changes of the second kind and lambda- 
points. Other types of interaction may lead to 
phase changes of the first kind, which are also 
observed in crystals with degrees of molecular 
rotational freedom. 

We shall now undertake the solution of Eqs. 
(5), in the case of cubic crystals composed of 
diatomic molecules, under certain assumptions 
concerning the interaction hindering relative 
molecular rotation. We shall suppose the poten- 
tial energy of a pair of nearest neighbors depends 
only on the angle y between the molecular axes 
of figure and that the interaction of molecules 
which are not nearest neighbors can be neg- 
lected. The potential V.,%” may be expanded in 
a Fourier series in y with an arbitrary constant 
term, conveniently assigned the value zero, and 
an initial term depending on molecular sym- 


metry, 


(6) 


€ 
rar (my) 


if the molecule has an m-fold axis perpendicular 
to its axis of figure. For homonuclear diatomic 
molecules m has the value two and for hetero- 
nuclear diatomic molecules such as the hydrogen 
halides m is unity. ¢€ is height of the potential 
barrier between adjacent minima. For cubic 
crystals, it is of course independent of the 
lattice indices k and / and on the neighbor 
indices a and 3, since all nearest neighbor sites 
of an arbitrary lattice site are equivalent. We 
shall consider here only the heteronuclear case, 


| 

7 
p 


n 
al 


5) 


m=1, for which Eqs. (3), (5) and (6) yield 


Kii=—acosy; k+l 
=(@); k=/j (7) 
a=ze/2kT, 


where z is the number of nearest neighbors of an 
arbitrary site on one of the interpenetrating 
lattices, belonging to each of the remaining n—1 
lattices. For a body-centered lattice is equal to 
2 and 2 is 8; for a face-centered lattice n is 4 and 
z is 4. If orientations are referred to an arbitrary 
polar coordinate system, substitution of Eqs. (7) 
into Eqs. (5) yields 


log = 


n 
-a > f f cos ygi(0’, ¢’) sin 


m= ff exp | ff cos 
0 


Xsin vivde'| sin ddddy, (8) 


cos y=cos cos 8’+sin sin 3’ cos (g—¢’). 


Due to the degeneracy of the kernels, cos y, the 
solutions of Eq. (8) may be constructed without 
difficulty. 


=— sinh (as;,); (9) 
as 


where L(x) is the Langevin function, coth x—1/x. 
The transcendental Eqs. (9) determining the 
parameters S, may be put in an alternative 


TABLE I. The order parameter s as a function of 
the temperature. 


T/T. s T/Te s 
1.0 0.00 0.5 0.79 
39 4 84 
8 54 a .89 
.65 .93 
6 ae ok .96 
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form 


We shall consider here only those solutions of 
Eqs. (10) for which a@ is positive (minimum 
interaction in antiparallel orientation) and oa» 
vanishes. In the body-centered case (n=2), we 
have 

k=1, 2, (11) 


where S is the positive real solution, minimizing 
the free energy, of the equation, 


s=L(as). (12) 


In the face-centered case (n= 4), we have several 
alternatives, all of which are equivalent, two of 
the s; being equal to the positive real root s of 
Eq. (12) and two equal to the negative real root, 
—s, for example, 


Se=+(—1)*s; R=1, 2, 3,4. (13) 


In both cases, there exists a critical value of a, 
equal to 3, below which Eq. (12) has no real 
solutions different from zero, and above which a 
single nonzero positive real root s, of minimum 
free energy, exists. This situation leads to a 
phase change of the second kind with a transition 
temperature, 7., given by 


(14) 


The parameter s is tabulated as a function of 
T/T, in Table I.* Above the transition tempera- 
ture, the stable distribution functions g, have the 
uniform value 1/47, corresponding to a random 


* If Eq. (12) is to be solved graphically, the change of 
variable, y= as, is convenient. 


y/a=L(y). 


The intersections of a straight line through the origin, with 
slope 1/a, and the graph of L(y) provide the desired solu- 
tions of the above equation. When 1/a is greater than }, 
the limiting slope of L(y), the line and the curve intersect 
only at the origin, and s is constrained to the value zero. 
When 1/a is less than }, there are three intersections at 
zero and tas, where s is the equilibrium value of the 
order parameter. For actual computation it is more con- 
venient to assign a set of values to y and to calculate a and 
s by the parametric equations, 


a=y/L(y) 
s=L(y), 


from which a tabulation of s as a function of a may be 
obtained by analytical methods. 
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distribution in orientation on each lattice. Below 
the critical temperature, the distribution be- 
comes non-uniform, the molecules of half of 
the interpenetrating lattices partially oriented 
parallel to and the molecules of the remaining 
lattices antiparallel to the polar axis of the 
reference system. The orientation of this axis 
relative to the lattice vectors is of course arbi- 
trary, as is to be expected from the fact that the 
interaction between neighbors depends only on 
their relative orientation. In Fowler’s theory® of 
hindered rotation, distribution functions of the 
form exp (as cos #) were introduced in an ad hoc 
manner. Fowler’s surmise that s was the fraction 
of molecules librating turns out to be unjustified. 
Eqs. (8), (9) and (10) show that s is the mean 
value of cos # for molecules of the lattice under 
consideration. The transcendental equation de- 
termining s is therefore quite different from that 
of Fowler. 

The thermodynamic functions of the crystal, 
when expressed in terms of the order parameter s, 
have the form 


4a 
S/R=log = sinh as| —as*, 


as 


E/RT=—as*/2, (15) 


F/RT = —log sinh as |-+ast/2. 


as 


From Eqs. (12) and (15), the configurational 
contribution to the heat capacity C, may be 
calculated as 
a’s*L'(as) dL 
C,/R=————__,_ L(x) = —. (16) 
1—aL’(as) dx 
Eqs. (12) and (16) yield the following value for 
the heat capacity anomaly, AC,, at the transition 
temperature, 


AC,/R=5/2. (17) 


The heat capacity anomaly AC, of 5 cal./mole is 
considerably smaller than anomalies in the con- 
stant pressure heat capacity AC, observed at 
lambda-points which are attributed to hindered 
rotation. A theory based on two accessible 
relative molecular orientations rather than a 
continuous set only increases AC, to 6 cal./mole. 
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Moreover, the neglect of fluctuations in energy 
cannot be expected to introduce an error of 
more than about 20 percent in AC,, if it is of the 
same magnitude as in the superlattice case. It 
seems necessary to attribute the large observed 
ratios of AC,/AC, to a dependence of the barrier 
height on the expansion of the lattice.‘ 
A general formula relating AC, and AC, in 
phase changes of the second kind has been given 
by Bethe and the writer.® 


AC, (1+08,T.)? v de 
o=—-—, (18) 


where 8, and «x, are the coefficients of expansion 
and compressibility of the crystal just above the 
transition temperature and v is the molar volume. 
If o is sufficiently large, the ratio, AC,/AC,, may 
become very large. If the denominator of Eq. (18) 
becomes zero or negative, the transition of the 
second kind must be replaced by one of the first 
kind occurring at a higher temperature than 7.. 
If €is supposed to be proportional to a power —m 
of the lattice parameter, o is equal to m/3. 
Supposing m to be of the order of magnitude 10, 
we may estimate rather crudely that o is of the 
order of magnitude unity. However, since theo- 
retical calculation of ¢« and a is not yet feasible, 
it is necessary to leave them as parameters to be 
determined from experiment. 

If the molecules of the crystal are polar, the 
influence of rotational transitions on dielectric 
behavior is of considerable interest. It is first 
necessary to consider the influence of dipole- 
dipole interaction in the present theory and the 
possibility of electric Curie points. Thus far we 
have considered only short range interaction 
between the molecules, effectively limited to 
nearest neighbors in the crystal. Electric Curie 
points can be produced by short range interac- 
tions if « and a, Eq. (7), are negative. It seems 
likely that the few observed examples of electric 
Curie points are in fact due to short range inter- 
action, although for their discussion a more 
complicated potential than (€/2) cos y is doubt- 
less necessary. In cubic crystals, nearest neighbor 
dipole-dipole interaction makes no contribution 
to the kernels K;., of Eq. (7). However, due to 


4 Eisenschitz, Proc. Roy. Soc. A168, 546 (1938). 
6 Bethe and Kirkwood, J. Chem. Phys. 7, 578 (1939). 
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the long range character of dipole-dipole forces, 
they are found to make a contribution when 
summed over the entire crystal. If this is done, 
Eq. (7) is modified as follows, 


Ki=(B—a) cosy; k+l 
=B cos 7; k=l (19) 
B= fu?/AkT, 


where A is the volume of the unit cell, u the 
molecular dipole moment, and f is a factor 
depending on the shape of the crystal, ranging 
from zero for a spherical specimen to 87/3 for a 
slab bounded by two parallel planes of infinite 
extent. We shall not discuss in detail the solutions 
of Eqs. (5) with the kernels (18), but remark 
that electric Curie points (o9+0) are excluded if 
e>4fu?/zA. Eqs. (14) and (15) remain valid, the 
transition temperature and thermodynamic func- 
tions being determined by the short range param- 
eter a alone. Otherwise electric Curie points 
determined by both a@ and £6 are predicted. 
Unfortunately, Curie points are predicted by the 
above criterion for crystals in which they are 
not observed, for example the hydrogen halides. 
We must therefore conclude that the degree of 
approximation of the present theory is insufficient 
for the treatment of long range forces due to 
dipole-dipole interaction. However, the short 
range force theory based on Eq. (7) remains 
formally valid for a spherical specimen or for a 
specimen of arbitrary shape immersed in an 
infinite medium of its own dielectric constant. 
Since there is no evidence that rotational lambda- 
points are affected by the shape of the crystal, 
the difficulty appears to lié entirely in Eq. (19) 
and the manner in which dipole-dipole inter- 
action is introduced there. Were fluctuations in 
energy among states of specified rotational dis- 
tribution taken into account, we surmise that the 
long range effects of dipole-dipole interaction 
would be nullified leaving a residual short range 
effect which would appear in a higher approxima- 
tion of the present theory. . 
Direct application of the present theory to the 
polarization of a crystal in an external electric 
field is prevented by similar difficulties with 
dipole-dipole interaction. In fact, if the short 
range interaction is zero, the present theory 
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leads to results identical with the Lorentz theory 
of dielectric polarization. Fortunately, a general 
theory of dielectric polarization can be formu- 
lated in which the long-range effects of dipole- 
dipole interaction are taken properly into account 
at the beginning.*:7 This theory leads to the 
following formula for the mean dielectric con- 
stant ¢ of a polar crystal,* where v is the molar 
volume ap the optical polarizability of a molecule 
and 


(e-1)v  3eN 


€ 


is the average dipole moment of a macroscopic 
spherical specimen of the dielectric, immersed in 
an infinite medium of its own dielectric constant, 
produced by fixing one of its molecules in a 
specified orientation, and (u-t), is the mean 
value of the scalar product of wu and the dipole 
moment wu of the fixed molecule averaged over all 
orientations of the latter. In the original formula- 
tion of the theory for liquids, the final average 
was unnecessary since in that case u-p is inde- 
pendent of the orientation of the fixed molecule. 
Although in liquids u reduces to the sum of the 
moment of a molecule and that which it induces 
in its immediate environment, this is not in 
general true of crystals. Long range dipole-dipole 
interaction having been implicitly taken into 
account in Eq. (20), it is possible to obtain a 
moderately good approximation to (u-f)» with 
the rotational distribution functions of Eq. (9). 
For (u-@)a we may write, where (cos 7i;)a is the 


(21) 


mean value of the cosine of the angle y;; between 
a pair of molecules i and 7, averaged over the 
orientations of both, the sum extending over all 
molecules of the spherical specimen except the 
specified molecule i. The condition that the 
specimen be immersed in an infinite medium of 
its own dielectric constant, involved in the 
definition of f, seems at first irrelevant if the 

®L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 

7J. G. Kirkwood, J. Chem. Phys. 7, 911 (1939). 

* The mean dielectric constant is defined as one-third 
of the diagonal sum of the dielectric constant tensor. 


It is of course the observed dielectric constant of the 
crystal powder. 
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distribution functions of Eq. (9) are employed, 
until it is remembered that these functions are in 
fact only valid under these conditions and 
otherwise would require correction due to a 
reaction field produced by the fixed molecule i 
and the exterior boundary of the specimen. 
Making use of Eqs. (9) and (11) or (13), we 
calculate 

(cos (—1)**!s?, (22) 


where k and / are the indices of the lattices on 
which the molecular pair i and j are situated 
and the sum of Eq. (21) becomes 


(—1)! 
(23) 


which, with » equal to two or four, is exactly 
equal to —s’, so that (u-w)s for both body- 
centered and face-centered cubic crystals is 
equal to w?(1—s?) and Eq. (20) becomes 


(e—1) 3eN 


9 


From Eq. (24), we observe that the dielectric 
constant is a continuous function of the tem- 
perature but has a discontinuous temperature 
derivative since the parameter s of Eq. (12) hasa 
discontinuous temperature derivative. Above the 
transition temperature Eq. (23) reduces to the 
Onsager equation, which is a fairly good approxi- 
mation for such substances as the hydrogen 
halides. Below the transition temperature, the 
dipole contribution to the dielectric constant 
diminishes to zero as s increases to its saturation 
value, unity. 

Phase changes of the second kind possibly 
attributable to hindered molecular rotation were 
found for the solid hydrogen halides by Giauque 
and Wiebe.* Two lambda-points at 70°K and 
125°K were found in the heat capacity curve of 
hydrogen iodide, and three at 89°K, 113°K, and 
117°K in the heat capacity curve of hydrogen 
bromide. None were found for hydrogen chloride, 
but instead one transition of the first kind. 
Considering only the transitions of the second 
kind, we may ask which of them are directly 


’Giauque and Wiebe, J. Am. Chem. Soc. 50, 2193 
(1928); 51, 1441 (1929). 


aot (24) 


attributable to hindered rotation. The dielectric 
constant measurements of Smyth and Hitchcock® 
provide a fairly conclusive answer to this ques- 
tion. The two upper transitions in hydrogen 
bromide and the upper transition in hydrogen 
iodide produce no marked effect on the dielectric 
constant. Below the 89°K transition in hydrogen 
bromide, the dielectric constant decreases rapidly 
with diminishing temperature. Although the 
measurements on solid hydrogen iodide were not 
carried through the transition at 70°K, their 
trend suggests a behavior similar to that of 
hydrogen bromide. The transitions of interest in 
the light of the present theory therefore appear 
to be the transition at 89°K in hydrogen bromide 
and the transition at 70° in hydrogen iodide. 
However, before discussing them, a few remarks 
concerning the other transitions should be made. 
The special solutions of Eq. (5) for cubic crystals 
predict only one rotational lambda-point. For 
rhombic crystals more than one lambda-point 
could result. But, since the hydrogen halide 
structures are face-centered cubic or face-centered 
tetragonal with an axial ratio differing only 
slightly from unity, and since the other transi- 
tions should have a marked effect on the di- 
electric constant if they were purely rotational, 
it seems necessary to attribute them to other 
causes, possibly to a coupling between molecular 
rotation and the lattice vibrations. The occur- 
rence of a transition of the first kind in hydrogen 
chloride, evidently to be associated with hindered 
rotation from the dielectric constant evidence, is 
difficult to understand unless the dependence of « 
on lattice expansion is admitted. On the basis 
of the present theory it could be attributed to a 
torque potential of different form than (€/2) cos y. 
It is interesting to remark the solutions of Eq. (5) 
with a two-minima potential, (€/2) cos 2y, lead 
to a phase change of the first kind. Such a 
potential should be appropriate to homonuclear 
diatomic molecules such as oxygen or nitrogen, 
but is difficult to see why hydrogen chloride 
should conform to it. It seems more likely that 
the parameter o of Eq. (18) is sufficiently large 
to replace a transition of the second kind by 
one of the first kind in the case of hydrogen 
chloride. 


9Smyth and Hitchcock, J. Am. Chem. Soc. 155, 1830 
(1933). 
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Fic. 1. The dielectric constant of solid hydrogen bromide 
as a function of temperature. Curve A, measurements of 
Smyth "7 Hitchcock at 20 kc. Curve B, theoretical, 
Eq. (24). 


Identifying them as the rotational lambda- 
points, we may use the present theory to in- 
terpret the transition at 70°K in hydrogen iodide 
and the transition at 89°K in hydrogen bromide. 
We shall treat the lattices of both as face-centered 
cubic. The hydrogen iodide’: !! lattice is face- 
centered tetragonal with an axial ratio of 1.08, 
thus very nearly cubic. The lattice of hydrogen 
bromide has not been determined with certainty, 
but it appears to be face-centered cubic above 
the transition, changing to one of lower sym- 


metry at lower temperatures. From Eq. (4) with 


z equal to 4, and the values of the transition 
temperatures, we calculate the barrier height € as 
210 cal./mole for hydrogen iodide and 267 
cal./mole for hydrogen bromide..In both cases 
the ratio AC,/AC, is large, about 4 for hydrogen 
iodide and 94 for hydrogen bromide, if the 
experimental values of AC, and the theoretical 
value of AC,, Eq. (17), are used. Unfortunately, 
no data on the coefficients of expansion and 
compressibility of the solid hydrogen halides 
seem to be available. In the case of hydrogen 
iodide, a value of approximately 1.6X10-* for 
the mean expansion coefficient between 82°K 
and 100°K may be calculated from the de- 
pendence of the lattice parameters on tempera- 
ture. The molar volumes of solid hydrogen iodide 
and hydrogen bromide at their melting points are 
42 cc and 33 cc, respectively. In the absence of 
the necessary data, a very rough estimate of the 
magnitude of « may be obtained from the ratios 

Ruhemann and Simon, Zeits. f. physik. Chemie 


(B) 15, 389 (1932). 
4 Natta, Mem. Accad. Italia Chim. 2, 15 (1931). 
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AC,/AC,and Eq. (18), with a value 10- atmos. 
for x, and a value 10-* for 8+ assigned to both 
of the hydrogen halides. For hydrogen iodide o 
turns out to be about 6 and for hydrogen 
bromide 9. Only the order of magnitude of these 
values and the fact that o is greater for hydrogen 
bromide than for hydrogen iodide are significant. 

It is of interest to compare the temperature 
dependence of the dielectric constant predicted 
by Eq. (24) with the measurements of Smyth 
and Hitchcock on the dielectric constant of 
hydrogen bromide. In Fig. 1, the experimental 
values of ¢ at a frequency of 20 kc are compared 
with the theoretical values at zero frequency 
calculated by Eq. (24). The dipole moment of 
the HBr molecule is assigned a value 1.17 Debye 
units to bring the observed and calculated values 
of « into agreement at 186°K. Although the 
moment of an HBr molecule in the crystal may 
differ somewhat from its value, 0.78 Debye, in 
the vapor, the rather large value 1.17 doubtless 
reflects a lack of accuracy in the theory, due to 
neglect of local ordering of the relative orienta- 
tions of neighboring molecules in the crystal. 
The agreement between the theoretical and 
experimental curves is far from quantitative, but 
it will be observed that the kink at the transition 
temperature is reproduced by the theory. The 
theory fails to reproduce the very sharp peak in 
the experimental curve at the transition tem- 
perature, which is probably due to local ordering 
of orientations. It should of course be remarked 
that the experimental curve at 20 kc is not 
strictly comparable with the theoretical curve 
which is valid only at zero frequency. However, 
hysteresis effects prevent an unambiguous ex- 
trapolation of the experimental dielectric con- 
stants to zero frequency. 

We seem justified in concluding that the 
present theory provides a semi-quantitative de- 
scription of phase changes of the second kind 
produced by the ordering of molecular orienta- 
tions in crystals. From its application to the 
solid hydrogen halides, we estimate a barrier 
height for the relative rotation of neighboring 
molecules in the lattice of the magnitude of 
several hundred calories per mole. From the 
ratio of AC,/AC, the rough estimates of o indi- 
cate a marked sensitivity of the barrier height to 
the lattice parameter of the crystal, greater for 
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hydrogen bromide than for hydrogen iodide. 
We surmise that in hydrogen chloride, the 
sensitivity of the barrier height to the lattice 
parameter is the cause of a transition of the first 
kind instead of one of the second kind. In the case 
of hydrogen bromide, the behavior of the di- 
electric constant in the transition is qualitatively 
described. The theory is subject to corrections of 
two types. The first of these would involve taking 
into account energy fluctuations in states of the 
same orientational distribution. The second re- 
lates to deviations from classical statistical 


AVRAMI 


mechanics. The present theory must be regarded 
as inadequate when the separation between 
neighboring rotational energy levels in the 
crystal becomes large relative to kT. Evidence 
for the necessity of a quantum-mechanical 
theory of cooperative phenomena in the case of 
molecules of low moment of inertia at low 
temperatures is provided by the differences in 
the heat capacity curves of methane and the 
deutero-methanes, which cannot reasonably be 
attributed to differences in intermolecular forces 
between neighbors in the lattice. 
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Transformation-Time Relations for Random Distribution of Nuclei* 
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Following upon the general theory in Part I, a consider- 
able simplification is here introduced in the treatment of 
the case where the grain centers of the new phase are 
randomly distributed. Also, the kinetics of the main types 
of crystalline growth, such as result in polyhedral, plate-like 
and lineal grains, are studied. A relation between the actual 
transformed volume V and a related extended volume Vj ex is 
derived upon statistical considerations. A rough approxi- 
mation to this relation is shown to lead, under the proper 
conditions, to the empirical formula of Austin and Rickett. 


A. REsuME 


N the previous paper! we have developed some 

of the general theory of phase change with 

the following assumptions strongly supported by 
experiment. 

(I) The new phase is nucleated by tiny ‘‘germ 
nuclei” which already exist in the old phase and 
whose effective number NW per unit nucleation 
region can be altered by temperature and dura- 
tion of superheating. These germ nuclei, which 
may be heterogeneities of any sort, usually con- 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
1J. Chem. Phys. 7, 1103 (1939). Referred to hereafter 


as K.P.C. I. 


(Received July 5, 1939) 


The exact relation is used to reduce the entire problem 
to the determination of Vex, in terms of which all other 
quantities are expressed. The approximate treatment of the 
beginning of transformation in the isokinetic range is shown 
to lead to the empirical formula of Krainer and to account 
quantitatively for certain relations observed in recrystal- 
lization phenomena. It is shown that the predicted shapes 
for isothermal transformation-time curves correspond well 
with the experimental data. . 


sist of small particles of the subcritical phase 
or tough films of the latter surrounding foreign 
inclusions. 

(II) The number of germ nuclei per unit 
region at time t, N=N/(t) decreases from N in 
two ways: (a) through becoming active growth 
nuclei, with number N’=N’'(t), as a result of free 
energy fluctuations with probability of occur- 
rence ” (a function of temperature, concentra- 
tion, etc.) per germ nucleus; (6) through being 
swallowed by grains of the new phase whose 
linear dimensions are growing at rate G (also 
a function of temperature, concentration, etc.). 
Letting V=V/(t) represent the volume of the 
new phase per unit volume of space we are led to 
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the relation 


between the density of germ nuclei and trans- 
formed volume. For the density of growth 
nuclei we find 


It is very useful to express these relations as well 
as the others in terms of the characteristic time 
scale of the substance and process, defined by 
changing the unit of time so that ndt=dr. This 
gives 


N(r) =Ne[1— V(z)], (5’) 
N\(r)=N f e[1—V(z) (5"”) 
0 


Beyond +=7, roughly corresponding to the time 
when the germ nuclei are exhausted, i.e., where 


N=1=Ne-*[1-V()] 
these expressions are to be replaced by* 


N=0, 


(III) Besides (5’), we found another relation 
in that V is a functional V|[N(r)]| of N, ice., 
at any time 7, it depends upon the values of 
N(r) throughout the interval 0 to 7. Assuming, 
as is usually true in crystalline growth, that 
where one grain impinges upon another, growth 
ceases, we found the form and significance of the 
terms Vm ex in the expansion of the functional 


V=Viex— Ve ex t+ 
(9) 


(IV) If now it is assumed, as is plausible from 
the underlying mechanisms and as the experi- 
ments verify, that 

G/n=a 


*If we wished to be meticulous we should add to the 


above the contribution of the last germ nucleus. This is 


most likely to nucleate at ordinary time ¢+1/n or charac- 
teristic time 7+1. 
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is approximately independent of temperature 
and concentration in a given “isokinetic range,” 
for a given substance and crystal habit, we have 


Vn ox= f f 
0 0 


x ++ +a(7—Zm) |dom. 


_ In particular 


0 


for pseudospherical or polyhedral grains, where 
o is a shape factor, equal to 47/3 for spherical 
grains. dm, Vmezr, and therefore the transformed 
volume V, in the characteristic time scale 
depend upon m and G only through the ratio 
a, which has been supposed constant. There- 
fore, the relative internal history of the trans- 
formation is independent of temperature, con- 
centration, or other variables (in the range 
in which u and G are proportional) ; or, to restate 
our conclusion: For a given substance and crystal 
habit, there is an ‘isokinetic range’’ of temperatures 
and concentrations throughout which the kinetics 
of phase change in the characteristic time scale 
remains unchanged. This observation enables us 
to translate results of isothermal reaction curve 
analysis to the description of phenomena under 
conditions of varying temperature and con- 
centration. 

The calculation of Vm ex, and the determina- 
tion thereby (with Eqs. (9) and (5’)) of the com- 
plete kinetics of the process, was shown to de- 
pend upon complex statistical considerations. 
Fortunately, in the important case of random 
distribution of the grain centers, we may use a 
general method which, expressing all quantities 
in terms of a suitably defined V; .z alone, enables 
us to avoid these detailed statistical calculations. 
We note that the assumption of random distri- 
bution of grain centers throughout the volume 
should be applicable with good approximation to 
cases where the germ nuclei tend to segregate at 
the boundaries of a fine grained structure in the 
old phase, as with many solid-solid phase 
changes, especially in metals. The results of this 
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TABLE II. 

E3(-—r) E3(-—71) T E3(—r) 
0.1 0.0000041}} 2.4 0.914 11 171.33 
0.2 | 0.0000641)} 2.5 1.057 12 227 

- 0.3 0.000318 2.6 1.223 13 293.67 
0.4 | 0.000987 2.7 1.4025 14 372.33 
0.5 | 0.00236 2.8 1.5995 15 464 
0.6 | 0.00481 2.9 1.815 16 569.67 
0.7 | 0.00874 3.0 2.050 17 690.33 
0.8 | 0.01462 3.1 2.305 18 827 
0.9 | 0.02307 3.2 2.582 19 980.67 
1.0 | 0.0346 3.3 2.881 20 1152 
t4 0.050 3.4 3.204 25 2315.67 
1.2 | 0.069 3.5 3.553 30 4079 
1.3 | 0.095 3.6 3.923 35 6567.33 
1.4 | 0.124 3.7 4.322 40 9905.67 
1.5 | 0.161 3.8 4.748 45 14,219 
1.6 | 0.205 3.9 5.202 50 19,632 
1.7 | 0.258 4.0 5.685 60 34,259 
1.8 | 0.317 5.0 12.34 70 54,786 
1.9 | 0.387 6.0 23.00 80 82,212 
2.0 | 0.468 7.0 38.67 90 | 117,539 
21 0.560 8.0 60.33 100 | 161,766 
2.2 | 0.666 9.0 89.00 4 
2.3 | 0.783 10.0 125.67 oe) 73/3! 


theory should then check approximately with 
those of a theory which explicitly takes into 
account the distribution along grain boundaries. 

Before proceeding, we should point out that 
the special formulae (12), (15), and (16) given in 
the preceding paper apply only to the case of 
pseudospherical or polyhedral grains, i.e., those 
growing at approximately the same rate in all 
directions. In the general case, where the growth 
is considerably unsymmetrical, we have to deal 
with an interesting and complex problem. We 
shall not study it here but merely note that 
besides the case, given above, of what we can 
call polyhedral growth, we may have two im- 
portant situations: First, where growth is 
primarily along two dimensions so that the 
grain is plate-like in shape; and second, where 
growth is primarily along one direction when we 
shall call it lineal. In these two cases, the ex- 
pressions corresponding to (16) are 


Vi =0'a* f (z)dz, (16’) 


Vi (16’’) 


respectively. Here a represents the ratio, in each 
case, of a suitably chosen G to n. Like (16), these 
formulae are applicable only if we are in a range 


AVRAMI 


approximately isokinetic, i.e., for G/n constant. 
In general for most substances, as crystallization 
velocities are made higher the type of growth 
changes from polyhedral to lineal. 


B. RELATION BETWEEN V AND 


1. Approximation. The Austin-Rickett formula as 
a special case 


In any region, selected arbitrarily, let us 
designate as “‘nonoverlapped’’ that part of the 
volume in which there is no transformed matter, 
i.e., no grains or parts of grains of the new phase. 
Then we have that, on the average, the ratio of 
nonoverlapped volume v' to extended volume 1; ez 
of a randomly selected region is equal to the density 
of untransformed matter 1—V at that time, i.e., 


—=1-V. (17) 

V1 ex 
If now we suppose that the grains of the new 
phase themselves have formed at random during 
the reaction, we might hope to apply (17) to 
these on the average. We cannot, however, do 
this directly because any actual grain is not 
randomly located in the entire unit volume. Its 
position is subject to the restriction that its 
nucleus lay somewhere outside of the transformed 
volume V(z) at the time z of its origin. We can 
circumvent this complication in one of two ways. 
In either case we get the same result. First we can 
restore complete randomness by including in our 
consideration all the ‘‘phantom”’ grains associ- 
ated with germ nuclei which would have been 
activated into growth if they had not been ab- 
sorbed by the growing new phase. At time z, the 
total rate of grain formation, including phantoms, 
is Ne-* (per unit volume), and it is this value of 
N(z) which should be used in the calculation of 
the various extended volumes V,, ez. In particular 


0 


where, as in (6) of the previous paper, 


v(7,2) =V1 ex(7,2) 


2R. F. Mehl and W. Johnson in Metals Technology, 
T.P. No. 1089 (1939), have used arguments similar to 
those in this section in the special case of spherical grains 
—— at a constant rate per unit untransformed 
volume. 


7 ( 
n 
— 
Ov 


iS 


logy, 
lar to 
zrains 
»rmed 


is the extended volume, at time 7, of a grain 
which started at z. The contributions of the 
phantom grains to the nonoverlapped and 
actual volumes, V;’ and V, respectively, ave zero, 
so that there is no physical difference introduced 
by their consideration. Strictly, we should dis- 
tinguish between V,, 2 as defined in this way 
and as used before, but for brevity’s sake we 
employ the same notation. From now on, unless 
otherwise specified, we shall work with the 
phantom-included population in calculating 
Vm ex This way of thinking is particularly useful 
in proving the theorem of the Appendix. Now 
we may apply (17) to the nonoverlapped volume 
of an average grain (over the phantom included 
population) in relation to its extended volume. 
Summing both numerator and denominator over 
all the grains (including phantoms) in unit 
volume we get 


Vi! 
l ex 
or 
(7,2)d (17') 
= e~*v(7,2)dz 
tae 


An alternative approach, without the use of 
“phantom’’ grains, which leads to the same 
result is as follows: The nucleus or center of an 
actual grain originating at time z was, with 
equal likelihood, anywhere in the volume 
1—V(z) (and only in this region). Then, at a 
later time 7, when the untransformed volume 
has diminished to 1— V(r), the ratio of non- 
overlapped to extended volume of the grain is, 
on the average, given by 


vy 1- V(r) V1 ex 
Viex 1—V(z) 1—V(r) 1—V(z) 


If we multiply both sides of the equation by the 
number of (actual) grains forming in time inter- 
val dz, and integrate from 0 to 7, we get 


V1’ 
=N f e~*v(1,2)dz 
1-—V 


which agrees entirely with (17’). 

We shall now consider an interesting special 
case: For transformations in which the new 
Phase is very fine grained, and therefore the 
overlapping is small, we may write as a rough 
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approximation 
V/=V 
and, combining with (17’), this gives 
V/A- V)= Vi ex. (17’’) 


Under isothermal conditions, and for phases 
which remain uniform in concentration, growth 
of linear dimensions is frequently at a constant 
rate G. n too is constant and as before, with 
a=G/n and polyhedral growth, we get 


Vie™ (r—2)%edz.* 
0 


In assuming the constancy of G under isothermal 
conditions, we have supposed that, once acti- 
vated, the germ nucleus soon passes the region 
of slow growth beyond which its rate becomes 
constant. This amounts to neglecting any ‘“‘in- 
cubation period”’ in the growth nucleus. Inte- 
grating, we get 
2! 3! 18’) 


ia 


n> 


where we have introduced the convenient ab- 
breviations 


B =60G°N (18’’) 
1 z 
En(—x) -—| (x 


xm 


m! 


TABLE III. 

E2(—7) 
0.1 0.0048375 0.0001625 
0.5 0.10653 0.01847 
1.0 0.36788 0.13212 
2 1.13534 0.86466 
3 2.04979 2.45021 
4 3.01832 5.98168 
5 4.00674 8.49326 
6 5.00248 12.99752 
7 6.00091 18.49909 
8 7.00034 24.99966 
9 8.00012 32.49988 

10 9.00005 40.99995 

20 19 181 

100 99 
72/2! 


* The generalization to the nonisothermal case as in 
quenching, and to that of varying concentrations (as in 
precipitation from noneutectic solutions), has the same 
form in the isokinetic range, but the unit of 7 in ordinary 
time is variable. 
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(18’) is valid only up to r=7, the time cor- 
responding to the exhaustion of the germ nuclei. 
Beyond this, the upper limit of the integral 
should be replaced by 7, and we may express the 
result of integration by 


—(r—7) J}. 


Formulae similar to (18’) and (18), with 
replaced by E2 and Fi, respectively, are obtained 
when we suppose the grain growth to be plate- 
like or lineal as in (16’) and (16’’). For brevity’s 
sake we shall use the designation E(7) with the 
understanding that 


(E3(—r) for polyhedral growth 
E(r)=4E2(—71) for plate-like growth 
| for lineal growth. 


E;(—71), Ex(—7r), and E,(—7) are tabulated as 
functions of 7 in Tables II and III. 
For n very large, (18’) goes over to 


Br? /3!=oG NB. 


For 7 small, i.e., ¢ small compared with 1/n, it 


goes into 


In either case, and for different types of crystal 
growth, it approximates the form B?', where Bisa 
function of temperature only and & is a constant. 
(17’’) gives then 


V(t)/(A— V(t) = Be* 


and this is exactly the formula for isothermal 
transformation of super-cooled austenite into 
bainite found empirically by J. B. Austin and 
R. L. Rickett.* The reasons for the approximate 
independence of temperature of k are now clear. 
Any observed deviations from (17’’) can be 
understood as due to the approximations in- 
volved in deriving it from (17’). 


(17’”) 


2. Exact relation 


To obtain the exact relation between V and 
Vi ex we may employ the following argument: 
_ Applying the same reasoning as in the previous 
section, not to the volumes of the single grains, 
but to the nonoverlapped and extended portions 
of the increments of these grains in an element of 


3J. B. Austin and R. L. Rickett, Metals Technology, 
September 1938, T.P., No. 964. 
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time, we get for the average grain (including 
phantoms) 

dv; /dv, cx = 1 V 
since the nonoverlapped increment of a grain is 


nothing other than the increment in transformed 
volume of that grain. For the unit volume 


dV 
=1-V. (19) 


Integrating, and rearranging, this gives 
V=1-—e-VYiex (20) 


as the desired fundamental relation, which is 
applicable in the case of random distribution of 
the grains. 

The derivation of (20) is entirely independent 
of any assumption about isothermal or isokinetic 
conditions. Thus even in the case where the 
distribution of grain centers is not uniformly 
random, as for instance in solidification in an 
ingot mold where there is a temperature gradient. 
as long as the distribution is fine grained and 
locally random, we may expect to be able to apply 
(20). Viez and V are then functions not only of 
the ordinary time ¢, but also of the space coor- 
dinates which, together with ¢, determine the 
thermal history and the local characteristic time 
at any point. ; 

Substituting (20) into (5) and (5’), we get 


N=Ne-(rt+"1 ex), (21) 


f dz. (22) 
0 


Thus the entire problem of determining the 
kinetics of reaction has been reduced to finding 
Vi ex from (18) in any particular case. The iso- 
kinetic case will be examined in detail in this and 
the next paper. 

Considerable doubt might have been raised 
as to the validity of (19), on the following grounds 
chiefly. The fraction 1—V measures the ratio 
of nonoverlapped to extended volume of any 
volume region selected at random. We cannot too 
easily assume that it also measures the corre- 
sponding ratio in what is effectively a selected 


_ surface region, i.e., the layers of grain increments 


in an element of time. The assumption requires 


| 
| 
q 
a oO 
4 in 
de 
t 


therefore more rigorous justification. It does 
turn out to be correct, primarily because of the 
complete randomness which we assume. The 
rigorous justification given in Appendix I brings 
out a number of very interesting relations. 
Expressions, which are valid in our case of 
random nucleation, for all Vx ex, Vx, and V;’ in 
terms of V; <x alone, are also derived there. 


C. TRANSFORMATION—TIME PHENOMENA 
IN THE ISOKINETIC RANGE 


1. The beginning of transformation, isothermal 
and nonisothermal. Derivation of Krainer’s 
formula. Recrystallization phenomena 


The kinetics of phase change, as concerning 
transformation-time and nucleation-time rela- 
tions are completely determined by (20) to (22), 
when is known. In the isokinetic case, ez 
is given by (18’) up to r=7, and we have 


V(r) =1—e7 82, (23) 
N(r) = 821, (24) 
N'(r) =" e dz (25) 

0 


7=7 is determined by 
#+BE(7)=lg N (26) 
and for r>7 we have, from (18’’’) 


N=0 (or, more exactly, 1 up to 7+1), (24’) 


N=N'=N f ds, (25’) 


In all cases when 7 is very small, i.e., at the 
beginning of transformation, and for polyhedral 
grain growth, (23) goes over into 


V = 1 — G3 /4 (27) 
or, even more approximately, 
Brt ca®N 


(27’) 


indicating that the V—¢ curve begins like a fourth 
degree parabola. The physical significance of 
this is that in the early stages of transformation 
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the number of growth nuclei present is propor- 
tional to the first power of the time and their 
respective volumes to the third power giving for 
the total volume transformed a dependence on 
the fourth power of the time. 

(27’) shows that V begins to depart appreci- 
ably from zero, i.e., the transformation “‘begins” 
when r‘/4! becomes appreciably different from 
zero. Now, for ordinary magnitudes of £8, this 
happens rather abruptly around r=1. Thus we 
would expect that, generally, the time of ‘‘begin- 
ning” ¢, of isothermal transformation would be 
approximately equal to 1/n at that temperature, 
i.e., 

ty(T) =1/n(T) 


or, if the circumstances under which transforma- 
tion is taking place are nonisothermal, the 
condition for the “beginning”’ is 


r= <1 


and this is precisely the result obtained empiri- 
cally by H. Krainer* in his experiments on the 
decomposition of austenite in steels. The root 
of the formula is seen to be the dependence of V 
upon the exponential of a high power of 7, a 
quantity which changes its value very con- 
siderably in the neighborhood of 7=1. 

For constant 8 (isokinetic case), not only does 
the “beginning” time vary inversely with m at 
different temperatures, but also the temperature 
—lg time curves for different degrees of trans- 
formation are all parallel. We have indicated 
the form of these curves in the previous paper. 

It is especially interesting to consider recrystal- 
lization phenomena in metals in the light of (27’). 
Here, even if a remains approximately constant 
(isokinetic range), we may vary N considerably 
through the degree of deformation which the 
old phase is made to undergo. Now there is 
strong experimental evidence’ for believing that 
the “degree of deformation” in a polycrystalline 
material is measured by Nadai’s ‘‘octahedral 
shear.’’® This measure has the advantage of giving 
results which are independent of the mode of 


‘H. Krainer, Archiv f. Eisenhiittenwesen 9, 619 (1936). 
a M 39} F. Kaiser and H. F. Taylor, Trans. A. S. M. 27, 227 


5M. Nadai, J. App. Phys. 8, 205 (1937). 
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deformation, whether, for instance, tension, 
torsion or rolling. Nadai has given formulae 
whereby the octahedral shear may be computed, 
in any case, from the principal specific extensions 
or contractions. Thus, for instance, in deforma- 
tion by elongation in tension, the octahedral 
shear y, is related to the specific reduction q 
through the formula 


Y¥n=V2 Ig (1+). 


For small plastic deformations y, will evidently 
be proportional to g. For larger deformations, we 
must distinguish between the two. 

Now, an inspection of the data on the number 
of nuclei which appear as a function of the 
deformation, e.g., that of Kornfeld and Pawlow,’ 


? Kornfeld and Pawlow, Physik. Zeits. der Sowjetunion 
6, 537 (1934). 
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suggests a parabolic dependence of N upon q, 
Or Yn, 1.€. 

Since, by (27’), if we keep temperature and other 
things constant, ¢, varies inversely as the 4th 
root of N, we would expect to find 


1 
or —8lgtz,« lg «lg lg 
Yn 


which for small deformations can be written 
—8 lg tp «lg q. 


If the recrystallizing grains are growing lineally, 
instead of polyhedrally, ¢, varies inversely as the 
2nd root of N and we would expect to find 


—4 1g t «lg lg(1+q)~lg 


In the absence of data on the “time of 
beginning”’ we use, as a rough indication, the 
data on the so called “‘time of recrystallization.” 
This is generally taken to mean the time ¢o.5 re- 
quired for approximately 50 percent transforma- 
tion. As will be seen, this is not, by any means, 
always simply related to ty). In many cases, how- 
ever, it is found that the transformation, once 
started, proceeds rapidly, i.e., the rise of the 
S-shaped transformation-time curve is quite 
steep and fo.; does not differ very much from ¢,. 
Here we may expect the data to fit an inverse 
fourth or eighth power relation depending on 
whether the growth is lineal or polyhedral. 

In the experiments of Karnop-Sachs* on the 
recrystallization of stretched copper wires, the 
deformations were small enough to permit us to 
replace the octahedral shear by the specific 
reduction. In Fig. 4 we have plotted the loga- 
rithm of ‘“‘time of recrystallization”’ against the 
logarithm of percentage reduction in area as 
obtained by Karnop and Sachs. We note that 
the data, over the entire range, all plot on straight 
lines with remarkable accuracy. The slope of the 
line passing through the points taken (from 
several wires) at 370°C is exactly equal to —4, 
and the slope of the line through the points 
taken at 284° is equal to —8. This tends to 
verify our assumption that the number of germ 
nuclei is proportional to the square,root of the de- 


8 R. Karnop and G. Sachs, Zeits. f. Physik 60, 464 (1930). 
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gree of deformation, and suggests that the growth 
is lineal at the upper temperature and polyhedral 
at the lower temperature. 

If we suppose a to be approximately constant 
and fix a definite time of annealing, (27) indi- 
cates that the N required to produce appreciable 
recrystallization varies inversely as n‘ or, taking 
logarithms, 


lg N« —4 lg n=4(Q+A)/RT 


where, as shown in the first paper of this series, 
Q+A is the total energy (per mole) required for 
a germ nucleus to become a growth nucleus. 
We shall suppose that we are in the temperature 
range where this has become constant. Now if we 
assume, as above, that 


N« y,'« [Ig (1+g) } 
we may write, for a given time of anneal f,, 
Ig Ig (1+q) «1/T. 


This proportionality is well satisfied over a 
wide range of deformation as may be seen from 
the straightness of the lines in Fig. 5, where we 
have plotted 1/T against lglg (1+g) for the 
data of Oberhoffer-Oertel’ on electrolytic iron 
and of Hanemann" on hot forged mild steel. 
This result also constitutes evidence in favor 
of the assumption of proportionality of N to 
some power of the degree of deformation 


2. Transformation-time curves in the isothermal 
range 

Let us consider the entire V—¢ curve, for 
polyhedral grain growth as given by (23) and 
(23’). We first examine the case where N is very 
large, i.e., exhaustion of the germ nuclei does not 
occur until the end of transformation. Under 
such conditions we may consider two limiting 
cases. When n is very small (7 small), (27) will 
be an adequate representation over the entire 
transformation, and this corresponds to the 
special case considered by Mehl-Johnson.? On 
the other hand, for m very large, 7 will be large 
during most of the reaction, and (23) is repre- 


® Oberhoffer and Oertel, Stahl und Eisen 44, 560-561 
(1924); ibid. 39, 1061 (1919). 
” H. Hanemann and F. Lucke, Stahl und Eisen 47, 1117 


tisas) H. Hanemann, Zeitschrift fiir Metallkunde 17, 316 
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sented by 
(28) 


For intermediate values of m the dependence of V 
upon ¢ will lie between (27) and (28). 

If N is small so that exhaustion occurs near 
the beginning of transformation, i.e., for 


N 


which (26) becomes in this case, most of the 
region of transformation is represented by 


(28’) 


which is practically the same as (28). The occur- 
rence of the factor N—1 instead of N is due to 
the neglect of the last germ nucleus in the 
derivation of (23’). Had we included this we 
should have obtained a result identical with (28). 

The physical interpretation of the two limiting 
cases is readily given. (27) was seen to be inter- 
pretable in terms of a linear increase in number 
and a third power increase in extended volume 
of the growth nuclei with time. (28) evidently 
corresponds to the case where the great majority 
of the growth nuclei start near the beginning of 
transformation, and the growth in total ex- 
tended volume is due merely to their continued 
growth. A little analysis will show that all other 
cases, of nucleus exhaustion at intermediate 
stages of the reaction, give a dependence of V; ., 
upon ¢ somewhere in between # and f*. 

For plate-like and lineal growth an analysis 
similar to the above gives for limiting cases 
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Fic. 6. 
Plate-like, 2 small : 
Plate-like, 1 large: 
(30) 
Lineal, ” small 


Lineal, large: 
(32) 


all of which may be written, by analogy with the 
Austin Rickett formula, 


(20’) 


It is interesting to note that all of the cases 
(29) to (32) (except the last which has the shape 
of a first-order reaction curve) give the typical 
S-shape for V against ¢, starting slowly, rising 
rapidly, and then leveling off again. Since, in 
each, 68 and 7 appear only in the combination 
Br* or Bt* (k=1, 2, 3, 4), on a V—lg# plot all 
isothermal reaction curves for a particular value 
of k will have the same shape, differing only in 
position. Thus it will be convenient to designate 
as the shape of a V—t curve the form which it 
takes in such a plot. We have plotted the four 
type curves in this fashion in Fig. 6. We note 
the increasing flatness of the reaction curves as 
the growth changes from polyhedral to lineal. 
This is in conformity with the available experi- 
mental data, which shows that when the habit 
changes from a polyhedral to an “acicular” 
needle-like, lineal type (e.g., for lowered tem- 
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perature), the transformation-time curves become | 


flatter. The effect is quite apparent in the data 
of Bain and Davenport on various alloys." 

For purposes of rough comparison with experi- 
ment, the difference in shape between the various 
limiting cases may be indicated by taking the 
ratio of times for two fixed degrees of trans- 


formation. A convenient representative index is _ 


the ratio of times for 75 percent and 25 percent 
transformation, respectively. This is equal in 
the various limiting cases to the fourth, third, 
second, and first roots, respectively, of 


lg (1—0.75) 
lg (1—0.25) 


Thus we should expect to find observed values 
of the ratio lying between the corresponding 
limits for the various types of growth, i.e., 


to.7s 
1.48¢—= £1.69 for polyhedral growth. 
to.2s T 0.25 
to.75 
1.69< <2.2 for plate-like growth. 
to.25 
to.7s 
2.2 <— £4.82 for lineal growth. 
to.25 
TABLE IV. 
TEMPER- 
ATURE 10.751IN #076 TYPE OF 
STEEL °C SECONDS SECONDS  o.25 GROWTH 


A 340 500 365 1.37 


275 2000 1350 1.48  Polyhedral, small 
250 4900 3200 = 1.53 
200 25500 17500 1.46 
B 340 110 60 1.83  Plate-like 
320 130 70 1.86 
300 170 96 1.89 
275 290 155 1.86 
260 420 180 2.33 Lineal, small 
250 540 220 2.45 
C 340 33 16 2.08 
340 69 41 1.68 Polyhedral, large 
E 340 245 145 1.69 ane 
275 680 410 1.66 
250 880 390 2.26 Lineal, » small 
F 340 650 480 1.35 
275 4500 3300 = 11.36 
250 10100 7400 1.37 


1 Bain and Davenport, Trans. A.I.M.E. 90, 117 (1930); 
100, 13 (1932). 
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To test these conditions, we have examined the 
ratios as observed in a set of different steels 
studied very completely by Bain and Daven- 
port.” The results are tabulated in Table IV. 
We notice two things. First, the ratios for a given 
steel, at different temperatures within an appre- 
ciable range, are approximately the same. This 
is evidence that the phenomena in this range are 
isokinetic. Second, except for A at 340°C and 
for F, all of the ratios lie within the theoretically 
predicted limits. What is more, the values for A 
lie approximately at the polyhedral—n small 
limit, i.e., where 

Vi ex © 


those for D and E above 250°C at the polyhedral 
——n large limit where 


Vi 


and for C, B and E below 260° at the lineal—n 
small limit where 


Vi ex 


When we consider that the time dependences of 
Vi ex leading to the indicated limits were derived 
completely theoretically, these agreements are 
strong confirmation of the suggested mechanisms. 

The values for B above 260°C are, of course, 
within the allowed range. They correspond to the 
5/2th root of 4.82, and we have suggested that 
the growth is plate-like in this region. It may be, 
however, that plate-like growth is not common, 
and that this intermediate value is due to some 
other cause. The whole question will bear further 
experimental and theoretical investigation. It is 
interesting that the outstanding discrepancy, 
that of A at 340°C and F, corresponds almost 
exactly to the 5th root of 4.82. Whether this is 
connected with some time variation of ” which 
leads to a 5th power dependence of V; ex upon ¢ 
we cannot at present say. 

We have plotted the experimental data, asso- 
ciated by Table IV, with each of the limits, upon 
the corresponding limit curve in Fig. 6. We note 
the falling off of the experimental values from the 


2 The soqgeiiene and descriptions of the steels are 


given in the Trans. A.I.M.E. 90, 117 (1930). I wish here to 
acknowledge my indebtedness, not | to Messrs. Bain 
and Davenport for making available their original data, 
but also to Dr. J. B. Austin of the U. S. Steel Research 
werk vid his help in reading the values off the charts. 
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curve at the end of reaction in each case. A still 
more clearcut comparison of the data with the 
limit curves may be made on the basis of (20’), 
which implies that a plot of lg lg 1/1— V against 
Ig t should give a straight line with slope k. In 
fact, it is found that the experimental data plot 
even more accurately on these straight lines than 
on those given by Austin-Rickett. This is 
especially true at the beginning and end points. 
There is, however, the same falling off of the 
experimental points at the end of reaction as we 
observe in Fig. 6. We can trace this deviation to 
two factors both operating in the same direction: 
First, in writing 
Vi ox = Bt* 


with k& constant, we have ignored the fact that 
the dependence of V;.x upon ¢, as given by (18), 
occurs with a diminishing power of ¢ as r in- 
creases in E(r). Second, as apparent in (18’”), 
there is an additional dropping off of Vi ex from 
the limit curve after the exhaustion of the germ 
nuclei. To obtain an exact fit for the experi- 
mental data we should use (23) and (23’) with 
properly chosen values of 8 and 7. This, we have 
not attempted here. The deviation of experi- 
mental data from (17’”) or (20’), always in the 
same direction (falling off), is important in 
affording strong additional confirmation of the 
whole germ nucleation mechanism, for if we 
assumed only a constant rate of spontaneous 
nucleation per unit untransformed volume we 
should obtain only the limit curves (27) to (31). 


APPENDIX I 


Expression of Vz, Vz ., and V;’ in Terms of V; ., 


First we shall show how (20), and therefore 
(19), may be derived by taking into account all 
of the higher order terms in V. To do this we 
proceed as follows. Let vmx represent the ex- 
tended volume of an average region in which m 
grains overlap, and let v,,’ represent the volume 
of that part of the region in which only m grains 
overlap. Then, assuming complete randomness, 
we have by the same reasoning as that which 
leads to (17) 

Um’ /Um ex =1—V. 


Multiplying both numerator and denominator 
by the total number of regions of this sort in 
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unit volume, we get 
Van V, (I, 1) 
where V;,,’ and Vm ex have the significance dis- 


cussed. in the previous paper. If we put this 
expression for V;’ into Eq. (10’) of that paper 


(—1)?-"Cn?V > ex) (10’) 
we get 
(1-V) a= (—1)?-*"C,.?V > ex 
p=m 


which, upon rewriting and putting down to- 


_ gether with (9), gives us an infinite system of 


linear equations for the determination of the 
infinite set of V, ex in terms of V. 


(I, 2) 


0=VV next Vp ex. 


p=m+1 

The determinants involved in the solution of this 
system do not seem to have been studied in the 
literature. They are connected with the rather 
remarkable type of continued fraction to be 
discussed below. Since the matter is largely of 
pure mathematical interest we shall confine 
ourselves here to the consideration only of the 
determinants which appear in the solution for 
Viw 


1-1 --- (-—1) 
0 3 «++ 
V =—3 (—1)?-?C2? 
V —k «++ (—1)?-*C,? 


V =—3 (—1)??°C,? 
«-- (—1)?-*C,? 


Upon multiplying alternate rows and columns 
by —1 all the minus signs become plus. Ex- 
panding both numerator and denominator ac- 
cording to the elements of the first column, and 
dividing through, we get 

Vie= -fV), (I, 3) 


where f is the fraction 


36-- C2” 6 Co” 
V 4 +--+ V Cg? 


The numerator and denominator of this fraction 
separately approach infinity, but the ratio is 
finite for values of V in the range of physical 
interest between 0 and 1. For V=0 f is evidently 
equal to 3; for V=1 it goes over into a ratio of 
determinants of a type studied by Zeipel.” 
Terminating the rows and columns after any 
fixed number, e.g., —1, the numerator and 
denominator have the values p—1 and p+2, 
respectively. For p> the ratio approaches 1. 
f then varies between } and 1 as V goes from 0 
to 1, and we see how the Austin Rickett formula 
with f=1 can give an approximate representa- 
tion of (I, 3). We shall now proceed to show how 
Vi ex, and therefore f, may be obtained as a 
function of V. 

If, by convention, we set Vo ex, “the extended 
volume of no transformation,’ equal to 1, 
(I, 2) may be written in the compact form 


0= VV, at Vin exe 
m=k+1 
k=0,1---. (I, 2’) 


Now let us arbitrarily cut off each of the equations 

after the term containing V, ex. We may later 

allow p to become infinite. Our system becomes 
Pp 


m=k+1 
k=0,1---p—1. (I, 4) 


Let us add to each equation the last equation 
multiplied by 


(—1)?-# Ci? 


p—l 


Our system becomes 


0=VV eet ox 


m=k+1 


V 
+ ( 1) -—] 


k=0,1---p—2, (I, 4’) 


2V. v Zeipel, Lunds Universitets Arsskrift, 1 (1865); 
VIII (1871). 
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of which the last (p—2) equation has the form 


V 


If now we add to each of the equations in (I, 4’) 
the last equation multiplied by 


1———— 


(- 1) p—k-iC,P- 1 
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1— 
p—k-j+2 


our system becomes 


p—3 
0= VVi ext (- Vin ex 


k+1 


Vp-2 ex» 


k=0,1,---p—3. (I, 4”) 


Continuing in similar fashion, we find that 
the last remaining term in the kth equation 
after the jth operation of this type is 


V 
V p—k—j+3 


Vi-V/3 


in which there appears a remarkable novel type 
of generalized ascending-descending continued 
fraction. It is characterized by the property 
that from each and every branch except the 
terminating ones there are two branches of 
higher order. We may perhaps describe it as a 
“tree fraction,’’ and define its order as the order 
number of branchings, e.g., j in (I, 5). Upon 
letting the number of operations j equal p—1, 
the sole remaining equation for k=0 is 


Vi-V/5 

Vi-V/3 

2 Vi-V/4 

Vo 31-V/2 

V1i-Vv/3 
| 


and this determines V; .x in terms of V. 


In general a tree fraction, like that in the 
brackets, defines a function which could be 
evaluated approximately for any value of the 
parameter by cutting off the “‘tree’”’ after a 
certain number of branchings. In our case the 
function defined is a familiar elementary one 
as we shall now show. Let us designate the 
expression in the brackets by y. Then 


V=yVi cx. (I, 6’) 


Since, for our infinite tree fraction, each suc- 
cessive lower branch is always exactly the same 


in structure as the entire tree, and is therefore 
equal to y, we have the relation 


/ 


2y 3y 4y 


which may be rewritten 
V y2 y3 
y=1-—+—-~ 


2!y Sly? 


from which 
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2 1-V/2 
(I. 5) 
p—k—j+1 
V 3 1-—V/2 
2 V1-V/3 
= 
1865); 
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and 
y=—V/Ig (1—V). 


Upon substituting into (I, 6’) this gives 


—Ig (1-V)=Vi x (I, 7) 
which, upon rearrangement, gives 
V=1—e-"1 ex, (20) 


thus establishing the desired result. 
Comparing (I,7) with (I, 3), we find f in 


terms of V, 
Ig (1— 


an interesting result for the ratio of the two 
infinite determinants given above. 
Upon expanding (20) we get 


Vi ex® Vi 
V= Vi ex 
2! 3! 
Vi 
m! 


If we compare this with (9) the suggestion 
crops up that, in our case of complete random- 
ness, the relations 


Vi a> Vi ex*/k! 


hold. That this is true can be seen as follows. 
From the previous paper 


Vi= 
m=k 
from which 
Vil = Vi— Views 


and, combining with (I, 1) and (20), this gives 
(I, 8) 


Now we must compare this with another relation, 
the generalization of (19), which we shall now 
obtain. Applying the same reasoning as above. to 
increments in extended volume dv; ¢x and volume 


AVRAMI 


dv; of the average kth-order region, we find 


dv;, 
=1—V=e-VYiex 
dp, ex 
For unit volume 
dv, 
=e-Viex, (I, 9) 
dV, ex 


Differentiating (I, 8) and comparing with equa- 
tions of the form (I, 9) we get 


dV Vi ex d Vi ex 


which is evidently equivalent to the desired 


result 
] ex 


Vi 


(I, 10) 


To express the V;, also in terms of Vi cx we use 
the recursion formula (I, 8) successively to get 


= 1—e-"1 ex(1+ V; Vi. ox) 


and, substituting (I, 10) and rearranging, this 
becomes | 


Veer ex), (I, 11) 


where E;, has the significance defined in (18”). 7 
Also 


V; m 
= e~Viex, (I, 12) 
! 

which is evidently identical with (10’) in our 
case, since 

¢ 
bed ex” t 
p=m p! 
p—m t 
Viex™ @ Viex c 
= (-1)-™ —. ‘ 
m!  p=m (p—m)! 

This completes the process of expressing all the 
other geometrical quantities in terms of Vi ex tl 
alone. 
1 
1¢ 
Al 
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The Magnetic Susceptibility of Silver Sub-Fluoride 


Smon FREED, NATHAN SUGARMAN, AND RICHARD P. METCALF 
G. H. Jones Laboratory, University of Chicago, Chicago, Illinois 


(Received December 4, 1939) 


The metallic character of the crystals of silver sub-fluoride was studied in terms of the mag- 
netic susceptibility. It was found that the conduction electrons contribute a considerable 
diamagnetism which was related especially to the anisotropy of the crystals. 


RYSTALS of silver sub-fluoride (AgsF) were 
prepared electrolytically from a solution of 


‘silver fluoride according to the procedure of A. 


Hettich! as modified by Scholder and Traulsen.? 
At the beginning of the electrolysis, the plates 
depositing on the cathode resembled highly 
burnished brass in their luster, but later the 


‘reflections became predominantly dark green. 


The different colors were re lected from different 
faces of the same crystal.* 

There is no longer any question concerning the 
stoichiometrical relation between silver and fluo- 
rine. According to x-ray measurements, the 
silver atoms are in equivalent positions in a 
cadmium iodide lattice and at the same distance 
from each other as prevails between the silver 
atoms in metallic silver. The distance between 
the silver and fluorine is the same as the inter- 
ionic distance in silver fluoride (AgF). 

The magnetic measurements were undertaken 
to obtain information concerning the conduction 
electrons. The crystals conduct the electric cur- 
rent about as well as graphite.‘ In addition they 
rectify an alternating current.’ Such properties 
corresponding to partially bound electrons might 
be expected to lead to a net paramagnetism 
somewhat higher than furnished by free elec- 
trons, not only because of a higher paramagnetic 
contribution by partially bound electrons but 
also because of the diminished diamagnetic sus- 
ceptibility.® 

The magnetic susceptibility was measured by 
the standard Gouy method at room temperature 


1A. Hettich, Zeits. f. anorg. Chemie 167, 67 (1927). 

*R. Scholder and K. Traulsen, Zeits. f. anorg. Chemie 
197, 57 (1931). 

*H. Steinmetz and A. Hettich, Zeits. f. anorg. Chemie 
167, 75 (1927). 

*H. Terry and H. Diamond, J. Chem. Soc. 2820 (1928). 

°N. F. Mott and H. Jones, Properties of Metals and 
Alloys (Oxford, 1936), p. 209. 


and at 77°K by an adaptation for low tempera- 
tures resembling one already described. The 
curve shown in Fig. 1 represents the apparent 
change in weight as the current was varied in 
the electromagnet. The calibration was derived 
from the change in weight observed when water 
filled the tube in place of the crystals. Correc- 
tions for the empty tube were determined at 
77°K since the method of substitution by water 
was no longer applicable. The molal suscepti- 
bility of silver sub-fluoride (Ag2F) was found to 
be —64.3+1.2 e.m.u.X10~-* both at room tem- 
perature and at 77°K. 

The molal susceptibility of silver fluoride 
(AgF) as computed from the measurements of 
Kido’ are —26.2 (Agt)+(—12.2)(F-) = —38.4 
(AgF). On the other hand, Brindley and Hoare® 
obtained —23.6+0.06 for the molal suscepti- 
bility potassium fluoride (KF) whereas Kido’s 
value was —25.7 by direct measurement. We 
shall accept —37.0 as the average value for the 
molal susceptibility of silver fluoride. The sus- 
ceptibility of metallic silver has been measured 
by de Haas and van Alphen® and found to be 
—20.5+0.2 per gram atom, independent of the 
temperature from 289°K to 20.4°K. According 
to Wiedemann’s law, the molal susceptibility 
of the sub-fluoride would be —37.0+(—20.5) 
= —57.5. However, the ideal independence of the 
constituent magnetic susceptibilities would also 
allow for an increase in the paramagnetism of 
the conduction electrons because of the increase 
in molal volume—the difference between the 
molal volume of silver sub-fluoride and of silver. 


*S. Freed, J. Am. Chem. Soc. 52, 2702 (1930). 

7K. Kido, Sc. Repts. Tohoku Imp. Univ. (1) 21, 149 
(1932); Sc. Repts. Tohoku a3 Univ. (1) 21, 288 (1932). 

8G, W. Brindley and F. E. Hoare, Proc. Roy. Soc., 
London 49, 342 (1937). 

® W. J. de Haas and P. M. van Alphen, Proc. Acad. Sci., 
Amsterdam 36, 263 (1933). 
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molecule=1, » is the Bohr unit of magnetic 
moment and AV! is the increase in the } power 
of the molal volume. The factor 3 before the 
bracket includes the subtraction effected by the 
diamagnetism of free electrons the magnitude of 
which is one-third of the spin paramagnetism.” 
The additivity on an “‘ideal’’ basis leads to the 
susceptibility —57.5+5.6= —51.9, a value which 
differs from — 64.3 the observed susceptibility by 
considerably more than any experimental un- 
certainties. 

The plate-like form of the crystals suggests 
that the enhanced diamagnetism is to be related 
to their anisotropic character, that is, to the 
“bindings” of the electrons which differ in 
magnitude with the directions in the crystal. 
If the surface of constant energy can be repre- 
sented by a spherical surface, the ‘‘binding”’ or . 
the increase in the “effective mass” of the 
electrons would, as we have mentioned in the 
introduction, have increased the paramagnetism 
and decreased the diamagnetism. However, 
anisotropic binding may lead to large diamag- 
netic contributions." Instead of the Landau 
formula for the diamagnetism of free electrons 


3h? 


the diamagnetic susceptibility becomes 


where a, a2, a3 represent the bindings in different 
directions, a3 corresponding to the direction of 
the magnetic field. This formula was derived on 
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If the electrons be regarded as free, the addi- 
tional net paramagnetic susceptibility is 


(=) av'|- 2(1.88A V2) 


= 3[1.88(27.7!— 10.3!) ]=5.6, 


where m is the mass of the electron, / is Planck’s 
constant, ” is the number of free electrons per 


the supposition that the surfaces of constant 
energy could be approximated by ellipsoidal 
surfaces, instead of spherical surfaces. It is clear 
that the diamagnetism becomes very large when 
the ratio in the brackets is large. In the case of 
silver sub-fluoride, a rather weak anisotropy can 
account for the diamagnetic susceptibility which 
was obtained experimentally. 


Landau, Zeits. f. Physik (1930); N. F. 


Mott and H. Jones, reference 5, gh Fy 
11H. Jones, Proc. Roy. Soc. A147, 506 (1934); N. F. 
Mott and H. Jones, reference 5. 
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Lattice Vibrations of Crystals and the Corresponding Vibrations 
of Their Solutions 


S. I. WEISSMAN AND SIMON FREED 
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(Received December 4, 1939) 


Crystals of hydrated europium chloride possess structures in their absorption spectrum which 
are faint but accurate repetitions of the intense patterns arising from electronic transitions in 
in the europium ions. These repetitions occur on both sides of the intense patterns at intervals 
which are ascribed to oscillations of the lattice, thermally as well as optically excited oscillations. 
The spectra of aqueous solutions of europium chloride and also of europium nitrate closely 
parallel the structures in the spectrum of the crystals. In the region where the crystals absorb, 
the solutions absorb also. The structures in the spectra of the solutions are somewhat more 
diffuse than the corresponding structures in the spectra of the crystals. Especially is this true 
for the structures at the intervals corresponding to lattice vibrations. 


HE spectra of crystals of the rare earths 

consist of rather widely spaced groups of 
lines. The separate groups owe their origin, it is 
generally accepted, to the decomposition of the 
electronic levels of the rare earth ions by the 
electric fields of their environment, somewhat as 
a line in the spectrum of a gas separates into a 
group of lines in the presence of an external 
electric field. Recently H. Ewald! has taken the 
absorption spectra through thick layers of 
crystals at the temperature of liquid hydrogen. 
He examined neodymium chloride NdCl3-6H.O, 
the nitrate Nd(NOs;)3-6H2O the mixed mag- 
nesium neodymium nitrate 3Mg(NOs)e-2Nd- 
(NOs3)3-24H2O, the corresponding zinc salt 
and the sulfate 
Nd2(CO,)3:8H2O. He found most of the intense 
groups accurately repeated at intervals toward 
shorter wave-lengths. The intervals between the 
intense pattern and the repetitions which were 
less intense conformed in most instances to well- 
known vibrational frequencies. For example the 
patterns from the crystal of the hydrated nitrate 
were repeated at intervals which accurately 
matched the frequencies of the internal vibra- 
tions of the nitrate ion as derived from Raman 
spectra. The intervals associated with the vibra- 
tions of nitrate ion were the same in the mixed 
as in the simple salt. In addition, repetitions of 
the patterns existed at intervals equivalent to 
the vibrations of water molecules. And when 
heavy water was substituted for light water the 


'H. Ewald, Ann. d. Physik 34, 209 (1939). 


faint patterns were moved toward the main 
pattern, that is, the intervals now accurately 
represented the lower frequencies characteristic 
of heavy water. The intervals we have been 
discussing range from about 800 cm~ to 1600 
cm, There were, however, smaller separations 
from the main patterns, from 100 cm~ to about 
500 cm. At these intervals the faint structures 
were usually rather diffuse. When heavy water 
replaced the light water, these lines also moved 
toward the main pattern and the lines were 
nearer the main group from the mixed zinc 
nitrate than they were from the isomorphous 
mixed magnesium nitrate. That is, an increase 
in the mass of the constituents of the lattice 
brought about a lower frequency and for this 
reason, especially, Ewald identified these in- 
tervals with lattice vibrations. To summarize, 
it may be stated that the electronic transitions 
in these crystals couple with the internal vibra- 
tional frequencies of the molecules as well as 
with various modes of oscillation of the lattice. 
In neodymium chloride, only the intervals of 
the water and lattice vibrations were observed in 
the spectrum. Hellwege? confirmed this interpre- 
tation. Reststrahlen revealed frequencies in these 
crystals which agreed with the intervals in the 
visible spectrum. 

The crystals of europium homologous to those 
of neodymium possess such sharp absorption 
spectra at room temperature that we examined 
them for the phenomena Ewald described. The 


2K. H,. Hellwege, Zeits. f. Physik 113, 192 (1939), 
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spectra of hydrated europium chloride at room 
temperature showed faint repetitions of the main 
pattern (especially at 4650A) not only on the 
side at shorter wave-lengths but also on the side 
toward longer wave-lengths. The intervals in 
both directions were the same. The lattice fre- 
quencies 115 cm~ and 127 cm“ were especially 
prominent. When the crystals of the hydrated 
chloride were placed in liquid nitrogen the 
spectrum no longer contained the lines on the 
long wave-length side. These represent then 
thermally activated lattice vibrations. 

Rather concentrated aqueous solutions of 
europium chloride (0.15M/, 0.75M, 1.7M and the 
corresponding light paths 90 cm, 24 cm, 12 cm) 
possessed definite regions of absorption in the 
spectrum on both sides of the main electronic 
pattern at 4650A and at intervals which paralleled 
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closely the lattice vibrations in the crystals. 
The lines were diffuse and the resulting patches 
extended with especial intensity over a range 
about the width of the main pattern and were 
removed from the latter on both sides by the 
same interval, about 125 cm. The solutions of 
europium nitrate possessed practically the same 
features in absorption. From more dilute solutions 
the structures were more diffuse and appeared 
to extend over a greater range of frequencies. 

In as much as reststrahlen confirmed the 
presence of the lattice oscillations in the crystals, 
there is reason to expect that infra-red or Raman 
frequencies would register these quasi-lattice 
vibrations in the liquid. 

This investigation was made possible with the 
aid of the Penrose Fund of the American 
Philosophical Society. 
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